Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,251,437	0	C	T	31.0%	72.8 / 47.3	58	E211K (GAA→AAA)	yebC	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/0);  new base T (18/0);  total (58/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTT  >  minE/1251407‑1251478                               |                                          cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:709963/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:322198/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:304585/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:294717/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:375910/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:25404/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:393721/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:85839/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:190288/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:179726/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:447480/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:452127/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:509963/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:530132/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:81344/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:1057311/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACggg               >  1:1112854/1‑59 (MQ=255) cGATCAGACGCATCAGTTTCGGTACGGTTCCTGCATCCATATCAGCTTTGGTAGACggg               >  1:360942/1‑59 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATTGAAACtt  >  1:1095147/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:511603/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:73536/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:730660/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:656230/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:623215/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:116164/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:360803/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:13277/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:81824/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:850840/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:930706/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:406369/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:339826/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACt   >  1:951044/1‑61 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACt   >  1:88889/1‑61 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACt   >  1:73319/1‑61 (MQ=255)           cATCAGTTTCGGTGCGGTTTTTGCATACATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:262566/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGAc                  >  1:624084/1‑46 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:938065/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:967670/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:779052/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:922708/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:749540/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:744214/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:970960/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:681100/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:673595/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:542919/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:540411/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:449560/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:419134/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:396516/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:243184/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:214040/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:174621/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACtt  >  1:120253/1‑62 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACt   >  1:502582/1‑61 (MQ=255)           cATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACt   >  1:164761/1‑61 (MQ=255)           cATCAGTATCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACg                 >  1:1126785/1‑47 (MQ=255)                               |                                          CGATCAGACGCATCAGTTTCGGTGCGGTTTCTGCATCCATATCAGCTTTGGTAGACGGGATCATGGAAACTT  >  minE/1251407‑1251478

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