Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,967,207	0	G	A	45.0%	17.8 / 24.2	20	P475S (CCG→TCG)	tktA	transketolase 1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base G (11/0);  new base A (0/9);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 5.95e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGT  >  minE/1967182‑1967254                          |                                                cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGt                                >  1:925733/1‑43 (MQ=255) cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                   >  1:603906/1‑56 (MQ=255) cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                   >  1:970879/1‑56 (MQ=255) cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                   >  1:908490/1‑56 (MQ=255) cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                   >  1:841229/1‑56 (MQ=255) cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                   >  1:503030/1‑56 (MQ=255) cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                   >  1:461836/1‑56 (MQ=255) cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                   >  1:433551/1‑56 (MQ=255) cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                   >  1:428912/1‑56 (MQ=255) cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                   >  1:369472/1‑56 (MQ=255) cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                   >  1:189761/1‑56 (MQ=255)            gACGTGCTCGATCGACTGGTGGGTCGGGCCGTCTTCGCCCAGGCCGATCGAATCGTGGgtga  <  1:480844/62‑2 (MQ=12)            gACGTGCTCGATCGACTGGTGGGTCGGGCCGTCTTCGCCCAGGCCGATCGAATCGTGGgtga  <  1:1003303/62‑2 (MQ=12)            gACGTGCTCGATCGACTGGTGGGTCGGGCCGTCTTCGCCCAGGCCGATCGAATCGTGGgtga  <  1:620481/62‑2 (MQ=12)            gACGTGCTCGATCGACTGGTGGGTCGGGCCGTCTTCGCCCAGGCCGATCGAATCGTGGgtga  <  1:661560/62‑2 (MQ=12)            gACGTGCTCGATCGACTGGTGGGTCGGGCCGTCTTCGCCCAGGCCGATCGAATCGTGGgtga  <  1:894375/62‑2 (MQ=12)            gACGTGCTCGATCGACTGGTGGGTCGGGCCGTCTTCGCCCAGGCCGATCGAATCGTGGgtga  <  1:33402/62‑2 (MQ=12)            gACGTGCTCGATCGACTGGTGGGTCGGGCCGTCTTCGCCCAGGCCGATCGAATCGTGGgtga  <  1:924550/62‑2 (MQ=12)            gACGTGCTCGATCGACTGGTGGGTCGGGCCGTCTTCGCCCAGGCCGATCGAATCGTGGgtga  <  1:1030947/62‑2 (MQ=12)            gACGTGCTCGATCGACTGGTGGGTCGGGCCGTCTTCGCACAGGCCGATCGAATCGTGGgtga  <  1:948588/62‑2 (MQ=255)                          |                                                CGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGT  >  minE/1967182‑1967254

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