Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,049,235	0	C	A	69.8%	27.8 / 45.5	43	A94E (GCG→GAG)	ydhJ	undecaprenyl pyrophosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/13);  new base A (30/0);  total (30/13)
Fisher's exact test for biased strand distribution p-value = 2.73e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGC  >  minE/1049180‑1049250                                                        |                gatAACACGTCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1029369/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:2420202/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:941836/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:2616040/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:2840067/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:2942223/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:320786/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:343479/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:3550434/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:3631085/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:433016/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:540372/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:733047/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:742706/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:919929/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:937930/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1023719/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1891103/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1106439/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:117332/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1272781/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1398350/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1401771/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1450903/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1532349/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1603029/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1609088/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:1723076/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCacac           >  1:190801/1‑59 (MQ=255) gatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGAGCaca            >  1:176684/1‑59 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:668197/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:1074834/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:1100576/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:926090/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:894304/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:2673635/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:2681229/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:3336308/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:270175/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:3174304/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:3527735/62‑1 (MQ=255)          cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:3373086/62‑1 (MQ=255)          cccTTTCAGATCCCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGc  <  1:858560/62‑1 (MQ=255)                                                        |                GATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGC  >  minE/1049180‑1049250

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