Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	31,302	0	G	T	53.3%	‑3.6 / 21.2	15	C66F (TGC→TTC)	unknown	TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base G (4/3);  new base T (3/5);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCTTACGGCACAGGATTTGAGTTAGGCGGCGATTTAACATCGAATGACACATTGAAATTTTCCGCCATATTAAAACAAGACGTTCATCGTGAACATAATAGCAATGAACCAACCGCCGTCGATAAAGACCGCACTTATAGCTTCGGGTTGGAAAATACCTACCGTTTCAGCGAACGAACCAAATTGATTACCGGCATGAGTTTTGATCGCCGCGAAGCTAAACGGGCGGAAAATTATCAGGCACGCTGCCCGACCATGAAAGATCGCTATTCTTATCGTTTAAATCGGGCGGAACCCAATCCGTTCTTGGATCCTGAAATCGCCTACCACTACGAAGCCGGCTATACGCGCACCTTCGGCAATTGGTTAAGACTGGACGGCGCGTTGTTCTATTCCGAAGTGCGTGATGCCATTGAACTTGTGGCGCATCCGACAATTTCCAATAAAGAACAAAATCAAAACTACGGCAAAGAAGTGTTTAAGGGTGTGGA  >  1/31053‑31545                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    ccGCTTACGGCACAGGATTTGAGTTAGGCGGCGATTTAACATCGAATGACACATTGAAATTTTCCGCCATATTAAAACAAGACGTTCATCGTGAACATAATAGCAATGAACCAACCGCCGTCGATAAAGACCGCACTTATAGCTTCGGGTTGGAAAATACCTACCGTTTCAGCGAACGAACCAAATTGATTACCGGCATGAGTTTTGATCGCCGCGAAGCTAAACGGGCGGAAAATTATCAGGCACGCTg                                                                                                                                                                                                                                                     >  1:22315/1‑250 (MQ=255)   gCTTACGGCACAGGATTTGAGTTAGGCGGCGATTTAACATCGAATGACACATTGAAATTTTCCGCCATATTAAAACAAGACGTTCATCGTGAACATAATAGCAATGAACCAACCGCCGTCGATAAAGACCGCACTTATAGCTTCGGGTTGGAAAATACCTACCGTTTCAGCGAACGAACCAAATTGATTACCGGCATGAGTTTTGATCGCCGCGAAGCTAAACGGGCGGAAAATTATCAGGCACGCTGcc                                                                                                                                                                                                                                                   >  1:564382/1‑250 (MQ=255)      tACGGCACAGGATTTGAGTTAGGCGGCGATTTAACATCGAATGACACATTGAAATTTTCCGCCATATTAAAACAAGACGTTCATCGTGAACATAATAGCAATGAACCAACCGCCGTCGATAAAGACCGCACTTATAGCTTCGGGTTGGAAAATACCTACCGTTTCAGCGAACGAACCAAATTGATTACCGGCATGAGTTTTGATCGCCGCGAAGCTAAACGGGCGGAAAATTATCAGGCACGCTGCCCGAc                                                                                                                                                                                                                                               <  1:189357/251‑1 (MQ=255)                     gAGTTAGGCGGCGATTTAACATCGAATGACACATTGAAATTTTCCGCCATATTAAAACAAGACGTTCATCGTGAACATAATAGCAATGAACCAACCGCCGTCGATAAAGACCGCACTTATAGCTTCGGGTTGGAAAATACCTACCGTTTCAGCGAACGAACCAAATTGATTACCGGCATGAGTTTTGATCGCCGCGAAGCTAAACGGGCGGAAAATTATCAGGCACGCTGCCCGACCacaaacgttcgtt                                                                                                                                                                                                                                 >  1:650753/1‑238 (MQ=255)                             cggcgATTTAACATCGAATGACACATTGAAATTTTCCGCCATATTAAAACAAGACGTTCATCGTGAACATAATAGCAATGAACCAACCGCCGTCGATAAAGACCGCACTTATAGCTTCGGGTTGGAAAATACCTACCGTTTCAGCGAACGAACCAAATTGATTACCGGCATGAGTTTTGATCGCCGCGAAGCTAAACGGGCGGAAAATTATCAGGCACGCTGCCCGACCacaaacgttcgttcgcaaagca                                                                                                                                                                                                                        <  1:726145/251‑22 (MQ=255)                                cgATTTAACATCGAATGACACATTGAAATTTTCCGCCATATTAAAACAAGACGTTCATCGTGAACATAATAGCAATGAACCAACCGCCGTCGATAAAGACCGCACTTATAGCTTCGGGTTGGAAAATACCTACCGTTTCAGCGAACGAACCAAATTGATTACCGGCATGAGTTTTGATCGCCGCGAAGCTAAACGGGCGGAAAATTATCAGGCACGCTGCCCGACCacaaacgttcgttcgcaaagcatt                                                                                                                                                                                                                      <  1:359303/250‑24 (MQ=255)                                  aTTTAACATCGAATGACACATTGAAATTTTCCGCCATATTAAAACAAGACGTTCATCGTGAACATAATAGCAATGAACCAACCGCCGTCGATAAAGACCGCACTTATAGCTTCGGGTTGGAAAATACCTACCGTTTCAGCGAACGAACCAAATTGATTACCGGCATGAGTTTTGATCGCCGCGAAGCTAAACGGGCGGAAAATTATCAGGCACGCTGCCCGACCacaaacgttcgttcgcaaagcatttg                                                                                                                                                                                                                    >  1:883860/1‑225 (MQ=255)                                  aTTTAACATCGAATGACACATTGAAATTTTCCGCCATATTAAAACAAGACGTTCATCGTGAACATAATAGCAATGAACCAACCGCCGTCGATAAAGACCGCACTTATAGCTTCGGGTTGGAAAATACCTACCGTTTCAGCGAACGAACCAAATTGATTACCGGCATGAGTTTTGATCGCCGCGAAGCTAAACGGGCGGAAAATTATCAGGCACGCTGCCCGACCacaaacg                                                                                                                                                                                                                                       >  1:613315/1‑225 (MQ=255)                                                                                                                                                                                                                               gccttggttttgcgaagaaaacccGCTTCCCGACCATGAAAGATCGCTATTCTTATCGTTTAAATCGGGCGGAACCCAATCCGTTCTTGGATCCTGAAATCGCCTACCACTACGAAGCCGGCTATACGCGCACCTTCGGCAATTGGTTAAGACTGGACGGCGCGTTGTTCTATTCCGAAGTGCGTGATGCCATTGAACTTGTGGCGCATCCGACAATTTCCAATAAAGAACAAAATCAAAACTACGGCaaa                      >  1:584594/24‑251 (MQ=255)                                                                                                                                                                                                                                             aagaaaacccGCTTCCCGACCATGAAAGATCGCTATTCTTATCGTTTAAATCGGGCGGAACCCAATCCGTTCTTGGATCCTGAAATCGCCTACCACTACGAAGCCGGCTATACGCGCACCTTCGGCAATTGGTTAAGACTGGACGGCGCGTTGTTCTATTCCGAAGTGCGTGATGCCATTGAACTTGTGGCGCATCCGACAATTTCCAATAAAGAACAAAATCAAAACTACGGCAAAGAAGTGTTTAAggg        <  1:13689/242‑1 (MQ=255)                                                                                                                                                                                                                                             aagaaaacccGCTTCCCGACCATGAAAGATCGCTATTCTTATCGTTTAAATCGGGCGGAACCCAATCCGTTCTTGGATCCTGAAATCGCCTACCACTACGAAGCCGGCTATACGCGCACCTTCGGCAATTGGTTAAGACTGGACGGCGCGTTGTTCTATTCCGAAGTGCGTGATGCCATTGAACTTGTGGCGCATCCGACAATTTCCAATAAAGAACAAAATCAAAACTACGGCAAAGAAGTGTTTAAggg        <  1:149221/242‑1 (MQ=255)                                                                                                                                                                                                                                             aagaaaacccGCTTCCCGACCATGAAAGATCGCTATTCTTATCGTTTAAATCGGGCGGAACCCAATCCGTTCTTGGATCCTGAAATCGCCTACCACTACGAAGCCGGCTATACGCGCACCTTCGGCAATTGGTTAAGACTGGACGGCGCGTTGTTCTATTCCGAAGTGCGTGATGCCATTGAACTTGTGGCGCATCCGACAATTTCCAATAAAGAACAAAATCAAAACTACGGCAAAGAAGTGTTTAAggg        <  1:300953/242‑1 (MQ=255)                                                                                                                                                                                                                                             aagaaaacccGCTTCCCGACCATGAAAGATCGCTATTCTTATCGTTTAAATCGGGCGGAACCCAATCCGTTCTTGGATCCTGAAATCGCCTACCACTACGAAGCCGGCTATACGCGCACCTTCGGCAATTGGTTAAGACTGGACGGCGCGTTGTTCTATTCCGAAGTGCGTGATGCCATTGAACTTGTGGCGCATCCGACAATTTCCAATAAAGAACAAAATCAAAACTACGGCAAAGAAGTGTTTAAggg        <  1:329269/242‑1 (MQ=255)                                                                                                                                                                                                                                                  aacccGCTTCCCGACCATGAAAGATCGCTATTCTTATCGTTTAAATCGGGCGGAACCCAATCCGTTCTTGGATCCTGAAATCGCCTACCACTACGAAGCCGGCTATACGCGCACCTTCGGCAATTGGTTAAGACTGGACGGCGCGTTGTTCTATTCCGAAGTGCGTGATGCCATTGAACTTGTGGCGCATCCGACAATTTCCAATAAAGAACAAAATCAAAACTACGGCAAAGAAGTGTTTAAGGGTGTgg   <  1:391868/247‑1 (MQ=255)                                                                                                                                                                                                                                                   acccGCTTCCCGACCATGAAAGATCGCTATTCTTATCGTTTAAATCGGGCGGAACCCAATCCGTTCTTGGATCCTGAAATCGCCTACCACTACGAAGCCGGCTATACGCGCACCTTCGGCAATTGGTTAAGACTGGACGGCGCGTTGTTCTATTCCGAAGTGCGTGATGCCATTGAACTTGTGGCGCATCCGACAATTTCCAATAAAGAACAAAATCAAAACTACGGCAAAGAAGTGTTTAAGGGTGTGGa  >  1:344631/4‑251 (MQ=255)                                                                                                                                                                                                                                                    cccGCTTCCCGACCATGAAAGATCGCTATTCTTATCGTTTAAAt                                                                                                                                                                                                                >  1:71535/3‑44 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    CCGCTTACGGCACAGGATTTGAGTTAGGCGGCGATTTAACATCGAATGACACATTGAAATTTTCCGCCATATTAAAACAAGACGTTCATCGTGAACATAATAGCAATGAACCAACCGCCGTCGATAAAGACCGCACTTATAGCTTCGGGTTGGAAAATACCTACCGTTTCAGCGAACGAACCAAATTGATTACCGGCATGAGTTTTGATCGCCGCGAAGCTAAACGGGCGGAAAATTATCAGGCACGCTGCCCGACCATGAAAGATCGCTATTCTTATCGTTTAAATCGGGCGGAACCCAATCCGTTCTTGGATCCTGAAATCGCCTACCACTACGAAGCCGGCTATACGCGCACCTTCGGCAATTGGTTAAGACTGGACGGCGCGTTGTTCTATTCCGAAGTGCGTGATGCCATTGAACTTGTGGCGCATCCGACAATTTCCAATAAAGAACAAAATCAAAACTACGGCAAAGAAGTGTTTAAGGGTGTGGA  >  1/31053‑31545

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