Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,081,210	(T)8→7	coding (885/1260 nt)	yccE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,081,203	0	T	.	84.0%	159.2 / 35.1	50	coding (878/1260 nt)	yccE	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (5/3);  new base . (24/18);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACATTATTTGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTCTTTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTACCCGAGTTTATA  >  REL606/1081128‑1081276                                                                            |                                                                          aGACATTATTTGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGctttct                                                                           >  2:1488935/1‑76 (MQ=255)  gACATTATTTGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGctttctt                                                                          <  1:820493/76‑1 (MQ=255)  gACATTATTTGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGctttctt                                                                          <  2:283382/76‑1 (MQ=255)   aCATTATTTGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGctttcttt                                                                         >  1:928246/1‑76 (MQ=255)      ttattTGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTCtttttt                                                                      <  1:818140/76‑1 (MQ=255)       tattTGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTCttttttt                                                                     <  1:1204990/76‑1 (MQ=255)        attTCGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATCTCATGGAGCTTTCttttttta                                                                    <  2:779147/76‑2 (MQ=255)          ttGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTCTTTTTTTaaa                                                                  <  2:666734/76‑1 (MQ=255)          ttGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTCTTTTTTTaaa                                                                  >  1:941357/1‑76 (MQ=255)          ttGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTCTTTTTTTaaa                                                                  >  1:1287721/1‑76 (MQ=255)            ggCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTCTTTTTTTAaaaa                                                                <  2:1385354/76‑1 (MQ=255)                 tGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAAt                                                          >  1:179024/1‑76 (MQ=255)                  gATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTCTTTTTTTTAAAAAACAAt                                                          >  2:810592/1‑76 (MQ=255)                     cTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCAt                                                      >  2:1118243/1‑76 (MQ=255)                     cTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCAt                                                      >  1:870434/1‑76 (MQ=255)                      tGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCAtt                                                     <  1:683699/76‑1 (MQ=255)                       gCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTg                                                    >  1:689381/1‑76 (MQ=255)                       gCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTg                                                    <  1:8103/76‑1 (MQ=255)                       gCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTg                                                    >  2:1096312/1‑76 (MQ=255)                       gCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTg                                                    >  2:1096326/1‑76 (MQ=255)                       gCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTg                                                    <  2:313573/76‑1 (MQ=255)                          aTCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATa                                                 >  1:584810/1‑76 (MQ=255)                              gACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTg                                             <  1:792355/76‑1 (MQ=255)                              gACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTg                                             <  2:786624/76‑1 (MQ=255)                              gACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTg                                             <  2:867084/76‑1 (MQ=255)                               acaacaAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGc                                            <  1:1415127/76‑1 (MQ=255)                                     aaaTTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAAt                                      <  2:661617/76‑1 (MQ=255)                                      aaTTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATa                                     <  2:1309060/76‑1 (MQ=255)                                      aaTTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATa                                     <  2:1310836/76‑1 (MQ=255)                                       aTTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATaa                                    <  2:2599/76‑1 (MQ=255)                                        ttCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATaaa                                   <  2:1047871/76‑1 (MQ=255)                                        ttCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATaaa                                   <  1:271617/76‑1 (MQ=255)                                          cAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGa                                 >  1:113689/1‑76 (MQ=255)                                              gtgtGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGaaaaa                             >  2:798988/1‑76 (MQ=255)                                              gtgtGGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGaaaaa                             >  2:893871/1‑76 (MQ=255)                                               tgtgGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAAt                            <  1:67200/76‑1 (MQ=255)                                               tgtgGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAAt                            >  1:1199569/1‑76 (MQ=255)                                                gtgGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATg                           >  1:734421/1‑76 (MQ=255)                                                 tgGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGt                          >  2:217679/1‑76 (MQ=255)                                                 tgGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGt                          >  1:1459609/1‑76 (MQ=255)                                                 tgGGTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGt                          >  1:1276535/1‑76 (MQ=255)                                                  tggTAGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTa                         >  2:702383/2‑76 (MQ=255)                                                      aGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCt                     >  1:106328/1‑76 (MQ=255)                                                      aGTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCt                     >  1:735289/1‑76 (MQ=255)                                                       gTCCCATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTg                    <  2:677712/76‑1 (MQ=255)                                                         cccATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCa                  <  1:576860/76‑1 (MQ=255)                                                          ccATGTCATGGAGCTTTCTTTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCa                  <  2:1302448/76‑1 (MQ=255)                                                          ccATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCAt                 <  2:786898/76‑1 (MQ=255)                                                          ccATGTCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCAt                 <  2:786909/76‑1 (MQ=255)                                                               tCATGGAGCTTTCTTTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTAc             >  2:1444837/1‑76 (MQ=255)                                                               tCATGGAGCTTTCTTTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTAc             >  2:1444831/1‑76 (MQ=255)                                                               tCATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGCTGCTTGAATAAAGAAAAATGTATCCTGCATTTAcc            >  1:544775/1‑76 (MQ=255)                                                                cATGGAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTAccc           >  1:171595/1‑76 (MQ=255)                                                                    gAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTACCCGAGt       >  2:335084/1‑76 (MQ=255)                                                                    gAGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTACCCGAGt       >  2:948145/1‑76 (MQ=255)                                                                     aGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTACCCGAGtt      >  2:487637/1‑76 (MQ=255)                                                                     aGCTTTC‑TTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTACCCGAGtt      >  2:990401/1‑76 (MQ=255)                                                                       ctttc‑tttTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTACCCGAGTTta    >  2:1433347/1‑76 (MQ=255)                                                                        tttctttTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTACCCGAGTTta    >  2:967118/1‑76 (MQ=255)                                                                          tctttTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTACCCGAGTTtata  >  2:729897/1‑76 (MQ=255)                                                                            |                                                                          AGACATTATTTGGCGTTGATCTGCAATCTGACAACAAAATTCAGGGTGTGGGTAGTCCCATGTCATGGAGCTTTCTTTTTTTTAAAAAACAATTCATTGATAGTTGCTTGAATAAAGAAAAATGTATCCTGCATTTACCCGAGTTTATA  >  REL606/1081128‑1081276

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑