Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,133,843	T→C	100%	intergenic (‑69/‑218)	ydiY ← / → pfkB	conserved hypothetical protein/6‑phosphofructokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,133,843	0	T	C	81.3%	111.8 / 20.0	48	intergenic (‑69/‑218)	ydiY/pfkB	conserved hypothetical protein/6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base C (39/0);  total (42/6)
Fisher's exact test for biased strand distribution p-value = 6.85e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACG  >  minE/1133776‑1133865                                                                    |                       tATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTa                      >  1:1673889/1‑70 (MQ=255) tATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAg                     >  1:2638716/1‑71 (MQ=255) tATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGCATCAAAGATTAg                     >  1:829772/1‑71 (MQ=255)            cGTACAACAAAAAAAGAGCCCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:2116199/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGCCCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:2523598/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCGTCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:988553/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTgg            >  1:2109336/1‑69 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:2004787/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:2066493/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1008890/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:2182409/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:24235/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:2616853/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:2635436/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:2763685/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:277056/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:2864112/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:2912400/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:459102/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:665172/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:672258/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:674556/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:741229/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1901593/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1024971/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:103521/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1256106/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1335764/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1403310/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1441305/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1459692/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1512759/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:155990/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1623422/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1633858/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1715780/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1716908/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1719006/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1863374/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1866577/1‑70 (MQ=255)            cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:1906818/1‑70 (MQ=255)            cGTAAAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGACTAGCGTCCCTGGa           >  1:124044/1‑70 (MQ=255)                    aaaaaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACg  <  1:1451465/71‑1 (MQ=255)                    aaaaaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACg  <  1:1416512/71‑1 (MQ=255)                    aaaaaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACg  <  1:21108/71‑1 (MQ=255)                    aaaaaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCATGGAAAGGTAACg  <  1:2414148/71‑1 (MQ=255)                     aaaaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACg  <  1:2719319/70‑1 (MQ=255)                     aaaaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACg  <  1:2738338/70‑1 (MQ=255)                                                                    |                       TATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACG  >  minE/1133776‑1133865

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