Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	771,770	0	T	C	72.7%	33.5 / 11.5	22	S39S (TCA→TCG)	hflD	predicted lysogenization regulator
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base C (16/0);  total (16/6)
Fisher's exact test for biased strand distribution p-value = 1.34e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGC  >  minE/771706‑771800                                                                 |                               tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTtgagcga                              >  1:1577669/1‑67 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:216609/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:479632/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:467030/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:3110738/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:2999492/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:2918715/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:2603346/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:723876/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:2225793/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:1856418/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:1830521/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:1738236/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:1601538/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:1538809/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:803933/1‑70 (MQ=255) tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGCGacac                           >  1:990577/1‑70 (MQ=255)                      ccGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTTTAGCGCATCGGCATCACaa     <  1:1409955/71‑1 (MQ=255)                      ccGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACaa     <  1:89895/71‑1 (MQ=255)                      ccGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACaa     <  1:967300/71‑1 (MQ=255)                       cGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACaa     <  1:1380273/70‑1 (MQ=255)                          cAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGc  <  1:2319580/70‑1 (MQ=255)                                   gCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACaa     <  1:1490564/58‑1 (MQ=255)                                                                 |                               TGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGC  >  minE/771706‑771800

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