Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	766,716	0	G	T	40.9%	‑3.9 / 13.1	22	A245S (GCG→TCG)	mngA	2‑O‑alpha‑mannosyl‑D‑glycerate specific PTS enzyme II
Reads supporting (aligned to +/- strand):  ref base G (6/7);  new base T (9/0);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 1.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑ACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCA‑AGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCAC  >  NC_000913/766494‑766960                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      aaGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑AATCAATAAGGCCAACGCCCGGATTTTTTCGGCTGAAGACGCCAACA‑Aggggggggg                                                                                                                                                                                                                          >  1:458641/1‑244 (MQ=255)        gCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGGTCACAC‑CAATAGTGCAAACGCCCGGATTTTTTCGGCTGGAATCGCCATCA‑AGAagagggggggtttt                                                                                                                                                                                                                   >  1:578147/1‑241 (MQ=255)                             tttCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGGTCA‑AATCAAAAGAGCAACCGCCCGTATTTTTTTGGCGGAAGGCGCCCCCACAGG‑GAGTGAGCGGTTTTATGGCATTcaccccctttcc                                                                                                                                                                                              >  1:311789/1‑240 (MQ=255)                               tCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑ACTCAATAGTGCAACCCCCCGTATTTTTTCGGCTGAAATCGCCCTCA‑AGGGGAGTGAGCGGTTTTATGGCAATCCCGCGCttttcgg                                                                                                                                                                                            >  1:199878/1‑247 (MQ=255)                                       gcCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑AACCAATAAAGCAAACCCCCGTATTTTTTCGGCTGGAGGGGCCATCA‑AAGAGAGGGGGCGGTTTAATGGCCTTTCCGCGCtctgagtggaaaagg                                                                                                                                                                                    >  1:120391/1‑237 (MQ=255)                                                           cGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑ACTCAAAAGTGGAACCGCCCGTATTTTTTCGGCCGGAGTCGCCATCA‑AGGAGAGGGGGCGGTTTAATTGGATTTCCGCGCCTTTAGGGGCTGTTTCCGGGCCGCTTCggcctgcc                                                                                                                                                                >  1:257608/1‑244 (MQ=255)                                                                                                      atatCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑ACTCAATAGAGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCCTCC‑AGGAGAGTGAGCGGTTTAATGGCATTTCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGCTGTTCCAACAAGCGCTTACCCCCAAGCGTAGCgctg                                                                                                                     >  1:122637/1‑248 (MQ=255)                                                                                                                        aCGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑ACTCAAAAATGCAACCGCCCGGATTTTTGCGGCTGAAATCGCCCTCA‑AGGAGAGGGAGCGGTTTAATTGCCTTCCCCCGCTTTTAGTGCCTGTTGCCGAGCCCATTCGCCATGCAGAAGCGTTGGTCCAACAAGCGCCTACCCTCCAGCGTAGCGGTgaaaagccaaccctacagc                                                                                                   >  1:377355/1‑234 (MQ=255)                                                                                                                            aaaCTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGGTCA‑ACTCAATATTGCAACCCCCTGTATTTTTTCGGCTGAAGTCGCCATccaagggagt                                                                                                                                                                                                                             >  1:280781/1‑117 (MQ=255)                                                                                                                                                 gTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑ACTCAAAAGTGCAACCCCCTGCATTTTTTCGGCTGAAGTCGCCATCA‑AGGAGAGGGAGCGGTTTAATGGCATTCCCGCGCTTTTCGTGCCTGTTGCCGAGCCGATTCGCCATGCCGAAGAGTTGATCCACCAAGCGCTTACCCTCAAGCGTAGCGATGGGAAGCGTAACGTACAGCAAGATaagccacccgggcaaaaggg                                                                          >  1:285363/1‑232 (MQ=255)                                                                                                                                                      aaaaCAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑ACTCAATAGCGCAACCGCCTGCATTTTTTCGGCTGAAGGCGCCATCA‑AGGAGAGTGGGCGTTTTTATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGCTTCGCCATGCAGAAACGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATAAGCAACCGGGGGAAAGTGTCaaa                                                                      >  1:76589/1‑250 (MQ=255)                                                                                                                                                      aaaaCAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑AACCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCCTCA‑AGGAGAGGGAGCGTTTTAATGGCATTCCCGCGCTTTCTGTGCCCGTTGCCGAGCCGCTTCGCCATGCAGAAGCGTTGATCCCACAAGCGCCTACCCTCCAGCGGAGCGATGAGACGCGTACCGGACAGCAAGAAACGCAAACGGGGAAAAGTGTCaaa                                                                      >  1:209540/1‑250 (MQ=255)                                                                                                                                                                 cTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑ACTCAATAGGGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCA‑AGGAGAGTGAGCGGTTTAATGGCATTTCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCCTGCAGAAGCGTTGATCCCACAAGCGCTTACCCTCAAGCGGAGCGATGAGACGCGGACCCTACAGCAAGATACGCAACCCGGGAAAAGAGTCaaaaagggggtgaa                                                           >  1:378044/1‑240 (MQ=255)                                                                                                                                                                                          gggtacaataacaaatggggacAACCGCCTGTTTTTTTGCGGCTGAAGTCGCCAGCA‑AGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGTGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGGGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATctct                                  <  2:367944/230‑1 (MQ=255)                                                                                                                                                                                             ggaatcaacaaaagtgggaaaCCGCCTGTATTTTTGCGGCTGAAGTCGCCATCA‑AGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCtttt                               <  2:70878/232‑1 (MQ=255)                                                                                                                                                                                                acaaactaaagagacagcccGCCTGTATTTTTGCGGCTGAAGTCGCCATCA‑AGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCGGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGcc                            <  2:138229/233‑1 (MQ=255)                                                                                                                                                                                                       aaagagtacacccGCCTGTATTTTTGCGGCTGAAGTCGCCATCA‑AGGAGAGTGAGCGTTTTAAGGGCATTCCCGCGCTTTCAGTGCCGGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGt                     <  2:208278/240‑1 (MQ=255)                                                                                                                                                                                                              gCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCA‑AGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGt              >  1:242350/1‑251 (MQ=255)                                                                                                                                                                                                               aacccGCCTGTATTTTTGCGGCAGAAGTCGCCATCA‑AGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGt              <  1:221248/247‑1 (MQ=255)                                                                                                                                                                                                                aaCCGCCTGTATTGTTGCGGCTGTAGTCGCCATCC‑TGGAGAGTGAGCGGTTTAATGGCATTCCCGCGCTTTCAGTGCCAGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCAAACAAGCGCTTCCCCTCAACCGAAGCGCTGAGTCGCGGACCGTACAGCAAGATACGCATCCGGTGACAAGGGTCACAACGGAACTGAAACAGGCCCTGTTGAACGGAATCTCTTTTGCCGTACCGTTGATTGTcg            >  2:215924/1‑251 (MQ=255)                                                                                                                                                                                                                     ccTGTATTTTTGCGGCTGAGGGCGCCATCA‑AGGAGAGTGAGCGTTTTAAGGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGGAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCgggg       <  2:189698/251‑1 (MQ=255)                                                                                                                                                                                                                          aTTTTTGCGGCTGAAGTCGCCATCA‑AGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCAc  <  1:155205/251‑1 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      AAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCA‑ACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCA‑AGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCAC  >  NC_000913/766494‑766960

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