Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,926,875	1	.	C	50.0%	17.0 / 18.7	12	coding (818/936 nt)	lysR	DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base . (0/6);  new base C (6/0);  total (6/6)
Fisher's exact test for biased strand distribution p-value = 2.16e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCGTCAGCCTGATCCG‑CCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAA  >  minE/1926860‑1926930                 |                                                        ccGTCAGCCTGATCCGCCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGt              >  1:168467/1‑60 (MQ=255) ccGTCAGCCTGATCCGCCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTg             >  1:1046416/1‑61 (MQ=255) ccGTCAGCCTGATCCGCCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTg             >  1:576571/1‑61 (MQ=255) ccGTCAGCCTGATCCGCCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTg             >  1:744294/1‑61 (MQ=255) ccGTCAGCCTGATCCGCCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTg             >  1:812057/1‑61 (MQ=255) ccGTCAGCCTGATCCGCCCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTg             >  1:901186/1‑61 (MQ=255)          tGATCCG‑CCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACaa  <  1:1003630/62‑1 (MQ=255)          tGATCCG‑CCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACaa  <  1:1156233/62‑1 (MQ=255)          tGATCCG‑CCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACaa  <  1:341002/62‑1 (MQ=255)          tGATCCG‑CCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACaa  <  1:572998/62‑1 (MQ=255)          tGATCCG‑CCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACaa  <  1:592415/62‑1 (MQ=255)          tGATCCG‑CCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACaa  <  1:872477/62‑1 (MQ=255)                 |                                                        CCGTCAGCCTGATCCG‑CCCCCTGCACCGCCCGTCATCAGCGCTGGTTCAGGCGTTTAGTGGGCATTTACAA  >  minE/1926860‑1926930

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