Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	188,326	T→A	100%	L294Q (CTG→CAG)	accA →	acetylCoA carboxylase, carboxytransferase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	188,326	0	T	A	75.0%	57.2 / 19.0	32	L294Q (CTG→CAG)	accA	acetylCoA carboxylase, carboxytransferase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/8);  new base A (24/0);  total (24/8)
Fisher's exact test for biased strand distribution p-value = 9.51e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAG  >  minE/188275‑188347                                                    |                      gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGCCGTGtt     >  1:1604022/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACgtgt      >  1:2014246/1‑69 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:910504/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:853912/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:821282/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:814660/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:543843/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:2769340/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:2460124/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:2033273/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:2015478/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1993462/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:183906/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1718241/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1727662/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1636358/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1616123/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1612825/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1311064/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1325290/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1330911/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1591384/1‑70 (MQ=255) gTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:1564129/1‑70 (MQ=255) gTAACCAGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCAGGCCGATCTCGACGTGtt     >  1:372492/1‑70 (MQ=255)    aCCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAg  <  1:1107755/70‑1 (MQ=255)    aCCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAg  <  1:3994/70‑1 (MQ=255)    aCCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAg  <  1:2387522/70‑1 (MQ=255)    aCCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAg  <  1:2323244/70‑1 (MQ=255)    aCCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAg  <  1:2128565/70‑1 (MQ=255)    aCCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAg  <  1:2048731/70‑1 (MQ=255)    aCCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAg  <  1:1646034/70‑1 (MQ=255)    aCCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGAAGTGTTAAg  <  1:1488870/70‑1 (MQ=255)                                                    |                      GTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAG  >  minE/188275‑188347

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