Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,596	0	G	A	48.7%	26.1 / 26.1	39	F13F (TTC→TTT)	ymgD	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (15/5);  new base A (0/19);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 7.71e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCACC  >  NC_000913/1222448‑1222745                                                                                                                                                         |                                                                                                                                                       cAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGaa                                                                                                                                                      <  2:66200/151‑1 (MQ=255)        ggCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAaaacaaac                                                                                                                                               >  1:409831/1‑151 (MQ=255)        ggCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAaaacaa                                                                                                                                                 <  1:13627/149‑1 (MQ=255)        ggCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAaaacaa                                                                                                                                                 <  2:927014/149‑1 (MQ=255)         gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAaaacaaac                                                                                                                                               >  1:357029/1‑150 (MQ=255)             tacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑c                                                                                                                                          <  2:681911/150‑1 (MQ=255)                 cattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGcc                                                                                                                                     >  2:127316/1‑151 (MQ=255)                      ccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTa                                                                                                                                 >  1:392188/1‑150 (MQ=255)                             attttCTGCACCCTGTGCGTTTTAAGTTTTTGCCT‑TTTCTTTGTTTTCTGCCCAAGTAATATGGTTTTCTGCCGGCAAC‑CCTTTTTTAATTTTG‑TATTCCTCCCCCCAGGTCATTG‑GTGCAAAAACAAACGGG‑CCCGCCAGTAAAGCAAAt                                                                                                                         <  2:678799/149‑1 (MQ=255)                                 tCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAAttttt                                                                                                                     >  2:760471/1‑151 (MQ=255)                                    ccacccggtgCGCGCTAAGTTTTAGCCT‑TTTGTTTTCGTTCTCCGAAAGTAATAGAGTTTTC‑GGCCGCAAA‑ACCTTTTTAAAATTGTTGGGCCGGCCCCCAGTTTA‑TG‑GGGGGAAAACAAACCGGCCCCGCCAGTAAAGCAAAttttt                                                                                                                     <  1:980062/141‑1 (MQ=255)                                        ccTTTGCGTTGTCA‑TTTTTGCCATTTTTGTGCGTTTCTTCAAAAATAATATAGGTTTCGGACAGCAAC‑CCTTTTTTAATATGG‑TGGCCCGTCCCCCAGATCTTTT‑GTGCAAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTtcatca                                                                                                               <  1:139133/150‑1 (MQ=255)                                           gtgCGCTTTCATTTTTTGCCT‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCTTTTTTAATATTG‑TAGCCCTGCGCCCAGGTCATTG‑GTGAAAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTtca                                                                                                                  <  2:293724/144‑1 (MQ=255)                                            tgCGCTAGCAGTTTTTGCCT‑TTTCGTTGGTTTCTGCGAAAGTAATATAGGTTTGTGA‑GGCAAACCCTTTTTTAATATTG‑TAGTCCGGCCCCAAGTTCATTG‑GTGAAAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTtcatcatca                                                                                                            <  2:387598/149‑1 (MQ=255)                                             gcgcTGTCATTTTTTGCCT‑TTTCGTTGGTTTCTGCGAAAGTAATATAGTTTTCTGACGGCAAC‑CCGTTTTAAATATTG‑TAGTCCTGCCCCAAGGTCATTG‑GTGAAAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTtcatcat                                                                                                              <  2:531653/146‑1 (MQ=255)                                               gcTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                       >  1:420455/1‑151 (MQ=255)                                                   tgggTTTTTGCTT‑TTTCTTTGTCTTCTACCCAATTATGATTGTTTTCAAACGGAAAT‑CTGTTTTTATTTGTT‑TGGCCCGGCTCCCGGTTCTTGG‑GTGCAAAAAAAAACAGG‑GCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTa                                                                                                     <  1:862757/146‑1 (MQ=255)                                                       tttttGCCT‑TTTCTTTGGTTTCTGCGCAATTAATATAGGTTTCTGACGGCAAC‑CCTTTTTTAATATTG‑TAGTCCTGCCCCCAGGTCTTGG‑GTGAAAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTa                                                                                                     <  2:597138/145‑1 (MQ=255)                                                       ttttgGCTA‑TTTCTTTGTTTTCTGCGCAATAAATATTGGTTTCGGACGACAAC‑CCTTTTTTAATTTTG‑TAGCCTTGCGCCCGTGTTTTGG‑GTGCAAAAACAACCAGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCa                                                                                                <  1:556591/150‑1 (MQ=255)                                                          ttGCCA‑TTTCGTTGGTTTCTGCGCAAGTAATATAGTTTTCGGACGGCAAC‑CCTTTTTAAATATTG‑TAGCCCTGCCCCCAGTTTTTTG‑GTGCAAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAttt                                                                                             <  1:122164/150‑1 (MQ=255)                                                            gCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAc                                                                                          >  2:546831/1‑151 (MQ=255)                                                             ccA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                         >  2:546153/1‑151 (MQ=255)                                                                      ttGGTTTCTGCACAGAAAATAGAGTTTGCGGACGCAAAA‑CCTTTTTTAATATGA‑TAGTCCTGC‑CCCCGGTTTTTGTGTGCAAAAACAAACAGC‑ACCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGaa                                                                                    <  1:526243/148‑1 (MQ=255)                                                                        ggCTTCTGCAAAAGTAATATAGCTTTCTGACGACAAA‑CCGTTTTAAATATTG‑GAGCCCTGCCCCAAGGTCATTG‑GTGCAAAAACAACCAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                                 <  2:171935/149‑1 (MQ=255)                                                                                               tttCTGACGCAAAC‑CCTTTTTTAATTTTG‑TAGCCCTGCCCCCAGTTTTTTG‑GTGAAAAAAAAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACgg                                                         <  2:885001/150‑1 (MQ=255)                                                                                                  cTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCAcc                                                     >  2:599751/1‑151 (MQ=255)                                                                                                         cAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATc                                              >  2:13996/1‑151 (MQ=255)                                                                                                           aa‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGc                                            >  2:388811/1‑151 (MQ=255)                                                                                                                   tttAATTTTG‑TAGCCCTGCCCCCAGGTCTTTG‑GTGCAAAAAAAAACAGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                      <  2:899350/150‑1 (MQ=255)                                                                                                                   tttAATATTG‑TACCCCTGCCCCCAGGTCTTGG‑GTGAAAAAAAAAACAGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                      <  1:652798/150‑1 (MQ=255)                                                                                                                          ttG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCa                               >  2:746945/1‑150 (MQ=255)                                                                                                                            g‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTg                            >  1:917467/1‑151 (MQ=255)                                                                                                                               aGTCCTCCCCCCAGTTAATTG‑TTGCAAAAAAAAACAGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGcc                                <  2:103095/145‑1 (MQ=255)                                                                                                                                  cctccccccAGGTAATTG‑GTGCAAAAAAAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGc                           <  2:825532/141‑1 (MQ=255)                                                                                                                                  ccTGCGCCCAGTTCATTG‑GTGCAAAAACAAACAGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCg                         <  1:536119/149‑1 (MQ=255)                                                                                                                                   cTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTACAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACCGCACCTTTAATCGCCCCTTCGAGCCCCTTCCCTGt                       >  1:57839/1‑150 (MQ=255)                                                                                                                                        cccAGATCTTTG‑GGGCAAAAACAACCGGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCa                     <  2:20202/147‑1 (MQ=255)                                                                                                                                         ccAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAAc                 >  2:722413/1‑150 (MQ=255)                                                                                                                                                   g‑gTGCAAAAAAAAACAGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAg            <  2:610970/145‑1 (MQ=255)                                                                                                                                                         gAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCAcc  >  1:585654/1‑150 (MQ=255)                                                                                                                                                         gAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCAcc  >  2:923770/1‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                       CAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTG‑TAGTCCTGCGCCCAGGTCATTG‑GTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCACC  >  NC_000913/1222448‑1222745

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑