Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,453,828	0	G	C	56.9%	21.8 / 64.2	51	V56V (GTG→GTC)	atoB	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (22/0);  new base C (0/29);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 6.41e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGT  >  minE/1453802‑1453871                           |                                            gTTGATGAAGTGATTATGGGTAACGTGTTACAAGc                                     >  1:724308/1‑35 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTg   >  1:1971291/1‑69 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:35298/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:2799819/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:2769086/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:2827492/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:2658999/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:2510818/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:2222118/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:2165169/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:2143420/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:367637/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:393313/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:444169/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:1032236/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:1101073/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:1120710/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:1197604/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:746859/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:1598563/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:1603329/1‑70 (MQ=255) gTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGt  >  1:1748257/1‑70 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:536251/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:555828/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:3298530/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:3275880/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:3226268/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:3158846/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:3088923/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:3017476/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:2843839/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:999209/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:78516/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:183187/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:1081979/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:1085918/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:1214682/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:1280623/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:1287640/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:1739451/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:1756104/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:2761040/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:1880182/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:2023473/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:2247341/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:232308/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:2360669/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:2573153/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:2575620/53‑2 (MQ=255)         aGTGATTATGGGCAATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:2609917/53‑2 (MQ=255)         aGTGATTATGGG‑AATGTCTTACAAGCAGGCCAAGGACAAAACCCTGCACGta           <  1:1279851/52‑2 (MQ=255)                           |                                            GTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGT  >  minE/1453802‑1453871

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: ‑