Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,323,554	0	C	T	72.6%	108.6 / 46.5	73	A124A (GCG→GCA)	glmM	phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base C (18/2);  new base T (53/0);  total (71/2)
Fisher's exact test for biased strand distribution p-value = 7.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGT  >  W3110S.gb/3323502‑3323571                                                     |                  aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCAt               >  1:145829/1‑57 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:1066499/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:897329/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:659799/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:486571/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:389532/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:2783222/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:2768755/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:2245932/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:2112225/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:2002392/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:1936857/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:1884710/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:1769037/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:1736420/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:168016/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc          >  1:1416040/1‑62 (MQ=255) aaTCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCAGc          >  1:2604195/1‑62 (MQ=255)        gCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTgg   <  1:2115818/62‑1 (MQ=255)        gCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTgg   <  1:1993008/62‑1 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCGTCTTCTACTGCATCCGGTAGTTTgg   >  1:1623160/1‑61 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTgg   >  1:1825403/1‑61 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTgg   >  1:2386227/1‑61 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:387622/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1057100/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2692514/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2743420/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1202046/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1169045/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2786095/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2869200/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2905733/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:34324/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:34376/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:34674/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:716185/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1612821/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1095237/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:427921/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:444734/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1078283/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:890940/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:545713/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:57572/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:624821/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:659604/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1066524/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2157208/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1450521/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1797301/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1821872/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1353188/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1287642/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1277425/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2069807/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1225214/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1642720/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2601719/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2198848/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:223965/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1217845/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2294741/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2339767/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2356833/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2405238/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2502405/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:2598402/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGGCAGTTTGGt  >  1:2629318/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTCCGCTTCGATGGCCTCTTCTACTGCATCCGATAGTTTGGt  >  1:138883/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCATTTACGCTTCGATGGCCTCTTCTACTGCATCCGGTAGTTTGGt  >  1:1625198/1‑62 (MQ=255)         cAGCTGATCTCTTTTTCCAGTTCCGCTTCGATGGCCTCTTCTACTGCATCCgg           >  1:793247/1‑53 (MQ=255)         cAGCTGATCTCTATTTCCATTTCCGCATCGAAGGACTCTTCTACTGCATCCGGGAGTTTGGt  >  1:2115599/1‑62 (MQ=255)         cAGCTGATCTCTATTTCCATTTCCGCATCGAAGGACTCTTATACTGCATCCGGTAGtt      >  1:2118881/1‑58 (MQ=255)                                                     |                  AATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGT  >  W3110S.gb/3323502‑3323571

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