Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	704,575	0	T	G	52.9%	6.7 / 13.5	17	V211G (GTT→GGT)	nagE	N‑acetylglucosamine‑specific PTS enzyme IIABC component
Reads supporting (aligned to +/- strand):  ref base T (2/6);  new base G (8/1);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 1.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCAT‑GTCCGGCTTCTTCCCGATCATGATGTTCGGTCTGCCGGGTGCGGCGCTGGCGATGTACTTCGCAGCACCGAAAGAGCGTCGTCCGATGGTTGGCGGTATGCTGCTTTCTGTTGCTGTTACTGCGTTCCTGACCGGTGTGACTGAGCCGCTGGAATTCCTGTTCATGTTCCTTG  >  NC_000913/704329‑704805                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        gcAAACGCTTGGTGCCGGTTGCCTCCGAGTTCTGCTGCCTGATGGTGGCCACGATTTTTGGTGACGGCTGGCCGCCGGTACAGCACGCTTTCCATGCAGGCGGCGAGTGGATCTTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCATCCGTCTGCTTATCCCACTCGGTCTGCATCAGGTACTGAACACCATCGCCTTGTTTCATATTAGTGAATTCTCCAACACCGCGGGTACGGTTTTc                                                                                                                                                                                                                                     <  2:144909/251‑1 (MQ=255)    aaCGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTGGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCAc                                                                                                                                                                                                                                  <  2:178666/251‑1 (MQ=255)                                        tGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGGTTTCCACGGTGGCATTAACCGCTTCTTTGCCGGTGACGGcccc                                                                                                                                                                                              >  1:138292/1‑248 (MQ=255)                                                                    cTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTTCACGGTGACATTAACCGCTTCTATTCCCGGGGCGGCCCCGCGGGGGTGGTCAT‑GTCCGGCTTCTTcc                                                                                                                                                                 >  2:38225/1‑251 (MQ=255)                                                                       gccgccGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGAAGGGTTTCCACGGTGACATTAAACGCTTCTTTGCCGGGGCCGGCACCGCGGGGGTGGTCAT‑GTCCGGCCTCTTcccc                                                                                                                                                               >  1:59794/1‑249 (MQ=255)                                                                        ccgccgGTACAGCACGCTAGCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTTGTGAATTCACCAACGCGGCGGGTACGGGTTTCCGCGGTGACATTAACCGCTTCTATTCCGGTGGCGGCACCGCGGGGATGTTCAT‑GTCCGGCTTCTTCCCGAt                                                                                                                                                             <  1:113036/251‑1 (MQ=255)                                                                                     aCGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCAT‑GTCCGGCTTCGTCCCGATCATGATGTTCGGt                                                                                                                                                <  2:278782/251‑1 (MQ=255)                                                                                         tATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTCCACGGTGACATTAACCGCCTCTAAGCCGGGGACGGCACCGCGGGGATGTTCAT‑GTTCGGCTTCTTCCCCATCATGATGTTCGGTCTGc                                                                                                                                            >  1:243108/1‑251 (MQ=255)                                                                                                       gcgAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTCCACGGTGACATTAACCCCTTCTATGCCCGGGACGGCACCGCGGGGGTGTTTAT‑GTCCGGCTTCTTCCCGCGCATGATGTTCGGGCTGGCGGGGGCGGCGCTg                                                                                                                              >  2:2079/1‑251 (MQ=255)                                                                                                          aGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGCGAATTCACCAACGCGGCGGGTACGGGTTTCCACGGTGACAATAAACGCTTCTTTGCCCGGGACGGGCCCGCGGGGAGGTTCATTGTC‑GGGTTCTTCCCGATCCTTATGCTCGGTCTTCCGGGTGCGGCGCTGGCg                                                                                                                           >  1:131607/1‑251 (MQ=255)                                                                                                                                gcTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGGGACGGCACCGCGGGGATGTTCAT‑GTCCGGCTTCTTCCCCATCATGATGTTTGGGCTGCCGGGTGCGGCGCTGGCGATGTACTTCGCAGCACCggaa                                                                                                      >  1:27536/1‑247 (MQ=255)                                                                                                                                         cGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCAT‑GTCCGGCTTCTTCCCGATCATGATGTTCGGTCTGCCGGGGGCGGCGCTGGCGATGTACTTCGCAGCACCGAAAGAGCgtcgtc                                                                                            >  1:207606/1‑251 (MQ=255)                                                                                                                                                                     ctgATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGGTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCCCCGCGGGGATGTTCAT‑GTCCGGCTTCTTCCCGATCATGATGTTCGGTCTGCCGGGGGCGGCGCTGGCGATGTACTTCGCAGCACCGAAAGAGCGTCGTCCGATGGTTGGCGGGATGCTGCTTTCTGt                                                                >  1:40251/1‑251 (MQ=255)                                                                                                                                                                         tCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTCCGGTTTTCCCCGGTGACCTTAACACATTCTTTGCCCGCGACGGCACCGCGGGGGTGTTCAT‑GTCCGGCTTTTTTCCGATCATGATGTTCGGTCTGCCGGGTGCGGTGCTTGTGATGTACTTCGCCGCCCCGGAAGAGGGTTGTCCGATGGTTGGCGGGCTGCTGCTTTCTGTTgtt                                                            >  2:124760/1‑249 (MQ=255)                                                                                                                                                                          cccAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCAT‑GTCCGGCTTCTTCCCGATCATGATGTTCGGTCTGCCGGGTGCGGCGCTGGCGATGTACTTCGCAGCACCGAAAGAGCGTCGTCCGATGGTTGGCGGTATGCTGCTTTCTGTTGCTg                                                           <  1:90444/251‑1 (MQ=255)                                                                                                                                                                                                                    ccAGATTGGCGAATTCACCAACGCGGCGGGTACGGTTTTCCGCGGTGACATTAACCGCGTCTATGCCGGTGACGGCACCGCGGGGATGTTCAT‑GTCCGGCTTCTTCCCGATCATGATGTTCGGTCTGCCGGGTGCGGCGCTGGCGATGTACTTCGCAGCACCGAAAGAGCGTCGTCCGATGGTTGGCGGTATGCTGCTTTCTGTTGCTGTTACTGCGTTCCTGACCGGTGTGACTGAGCCGCTGGAATTcc                 <  2:131607/251‑1 (MQ=255)                                                                                                                                                                                                                                   cACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCAT‑GTCCGGCTTCTTCCCGATCATGATGTTCGGTCTGCCGGGTGCGGCGCTGGCGATGTACTTCGCAGCACCGAAAGAGCGTCGTCCGATGGTTGGCGGTATGCTGCTTTCTGTTGCTGTTACTGCGTTCCTGACCGGTGTGACTGAGCCGCTGGAATTCCTGTTCATGTTCCTTg  <  1:38225/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        GCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCAT‑GTCCGGCTTCTTCCCGATCATGATGTTCGGTCTGCCGGGTGCGGCGCTGGCGATGTACTTCGCAGCACCGAAAGAGCGTCGTCCGATGGTTGGCGGTATGCTGCTTTCTGTTGCTGTTACTGCGTTCCTGACCGGTGTGACTGAGCCGCTGGAATTCCTGTTCATGTTCCTTG  >  NC_000913/704329‑704805

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