Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,830	0	T	A	52.4%	2.7 / 27.4	21	intergenic (‑385/‑4)	ylcI/nohD	DUF3950 domain‑containing protein YlcI/putative DNA‑packaging protein NohD
Reads supporting (aligned to +/- strand):  ref base T (3/7);  new base A (5/6);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGCCAGGGACTAT  >  NC_000913/580580‑581078                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         ttGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgt                                                                                                                                                                                                                                                          >  1:165667/1‑251 (MQ=255)                       attatCATTTTGTGGGTCCTTTCCTGCGATCCGACAGGTTTCGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACATGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAag                                                                                                                                                                                                                                     <  2:12612/250‑1 (MQ=255)                                             ccGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTg                                                                                                                                                                                                               <  2:132130/250‑1 (MQ=255)                                                          aGGTTACGGGGCGGCGACCACGCGGGTTTTCGCTCTTTATGAAAATTTTCCGGGTTAAGGCGTTTCCGTTCTTCTTAGTCATAACTTAATGTTTTTATTTAAACTACCCTCTGAAAAGAAAGGAAACGAAAGGTGCTGAAAGCGAGCTTTTTGGCCTATGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGAACAATGGAAGTCAACAACAAGCAGCTGGCTGACATTTTCGGTCAGAGTATCCGtccc                                                                                                                                                                                                 >  2:153248/1‑248 (MQ=255)                                                                                   gTTTTCGCTATTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGaa                                                                                                                                                                        <  1:127940/251‑1 (MQ=255)                                                                                                                 tAAGGTGTTTCCGTTCTTATTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAggga                                                                                                                                           >  1:65657/1‑249 (MQ=255)                                                                                                                 gcgcagctttCCGTTCTTCTTCGTAATAACTTAATGTTTTTATTTAAAATACCCTCTGATAAGAAAGGAAACGACAGGTGCTGTAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTATCAGTTTTTGTCCGTGGAATGAACAATGGAAGTCAACACAACGCAGCTGGCTGACATTTTCGGTGCGCGTATCCGTACCATTCAGAAATGTCAGGAACAGGGAATGCCCGTTCTGCGACGCGGTGGAAAggct                                                                                                                                           >  2:166713/8‑248 (MQ=255)                                                                                                                      cgtTTCCGTTCTTCTTCGTCATAACTTAATGTTTTTATTTAAAATACCCTCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGAACAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATg                                                                                                                                       >  1:202361/2‑249 (MQ=2)                                                                                                                                  tcttcGTCGTAACTTAATGTATTTATTTAAAAGTCCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTTGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTAGCCGTTCCATTCAGAACTGGCAGGAACAGGTAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTAt                                                                                                                            <  1:166202/248‑1 (MQ=255)                                                                                                                                     tcgtcgTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGAct                                                                                                                        <  1:262815/249‑1 (MQ=40)                                                                                                                                      cgtcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTg                                                                                                                     <  2:6815/251‑1 (MQ=40)                                                                                                                                                      atttatttAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAAGGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATg                                                                                                     <  2:165667/251‑1 (MQ=255)                                                                                                                                                        ttatttaAAATACCCTCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGAACAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATgg                                                                                                    <  2:144157/250‑1 (MQ=25)                                                                                                                                                           tttaAAATACCCCCTGAAAAGCAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCGGTCGTTTCCTGTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGAGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCAGTTCTGCAAGGCGGTGGCAAGGGACATGAGGTGCTTTATGAATCTGCCGCCGTCATAAACTGGTATg                                                                                                >  1:242278/1‑251 (MQ=14)                                                                                                                                                                 aTACCCCCTGAAAAGAAAGGAAACGCCAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCgca                                                                                           >  2:260931/1‑248 (MQ=255)                                                                                                                                                                                       aCGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGAACAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGa                                                                    <  1:71980/251‑1 (MQ=21)                                                                                                                                                                                          aCAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGAACAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGaaa                                                                  >  2:197725/1‑250 (MQ=25)                                                                                                                                                                                                  tGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGAACAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGcc                                                         >  2:193545/1‑251 (MQ=21)                                                                                                                                                                                                   gAAAGCGATCTTTTTGGCCTCTGATTTTTTCTTTTTCTGTTTTTGTCCGTAGAATGAACAATGGAAGTCAAAAAAAAGGAGCTGGCTGACATTTTCGGGGCGAGTATCCGTGCCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGTTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGGCATAAGATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTgcgc                                                          <  1:166713/249‑1 (MQ=11)                                                                                                                                                                                                                tgggCCTCTGTCGTTTCCTTTCTCTGTTTTGGTCCGTGGAATGAACAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATGGAGAACGAAAAGCTGCGCCGGGAGGGTgaaga                                            <  1:193545/248‑1 (MQ=255)                                                                                                                                                                                                                               tCCTTTCTCTGTTTTTGTCCGTGGAATTAACAATGGAAGTCAACAAAAAGTAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCGGCAGGCCAg                             <  1:209329/250‑1 (MQ=9)                                                                                                                                                                                                                                         gTTTTTGTCCGTGGAATGAACAATGGAAGTCAACAAAAAGGAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCATAACTGGCAGGAACAGGGAATGCCCGTTCTGCGATGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAAGTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCGGCAGGCCAGCGAGGCAGAt                   <  2:36540/250‑1 (MQ=11)                                                                                                                                                                                                                                                         ggtgCAAGGGAATTCAACAAAAAGCAGCTGGTTAACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGTTGTCAATGGTAATGAGGGGTTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGTGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGCCAGGGACTAt  <  1:153514/250‑1 (MQ=11)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         TTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAATGAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGCCAGGGACTAT  >  NC_000913/580580‑581078

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