Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,114,039	G→A	51.6%	I255I (ATC→ATT)	pps ←	phosphoenolpyruvate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,114,039	0	G	A	51.6%	‑4.4 / 24.8	31	I255I (ATC→ATT)	pps	phosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base G (11/4);  new base A (10/6);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

TCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACC  >  minE/1113983‑1114104                                                         |                                                                  tCCTGGGGCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGa                                                       <  1:1673022/69‑1 (MQ=255)  ccTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGa                                                       <  1:1295622/68‑1 (MQ=255)  ccTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGa                                                       <  1:3191502/68‑1 (MQ=255)  ccTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGa                                                       <  1:2029922/68‑1 (MQ=255)     gggTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTcg                                                   >  1:279784/1‑69 (MQ=255)      ggTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTcg                                                   >  1:3203186/1‑68 (MQ=255)                     cATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGttt                                   >  1:1133358/1‑69 (MQ=255)                      aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                  >  1:2718366/1‑69 (MQ=255)                      aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                  >  1:453238/1‑69 (MQ=255)                      aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                  >  1:3377156/1‑69 (MQ=255)                      aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                  >  1:2989575/1‑69 (MQ=255)                      aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                  >  1:2751304/1‑69 (MQ=255)                      aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                  >  1:2620023/1‑69 (MQ=255)                      aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                  >  1:2447483/1‑69 (MQ=255)                      aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                  >  1:1556679/1‑69 (MQ=255)                      aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTa                                  >  1:1667491/1‑69 (MQ=255)                      aTGCGGATTTTTTTAGACCCCATGGTGCGGCGCAAAATAGCCGGGCGATTCGCAGCAAGTGTCGGttt                                   >  1:2762440/1‑68 (MQ=255)                                             ggTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGtt           >  1:685895/1‑69 (MQ=255)                                             ggTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGtt           >  1:1340944/1‑69 (MQ=255)                                             ggTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGtt           >  1:1404951/1‑69 (MQ=255)                                             ggTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGtt           >  1:3060073/1‑69 (MQ=255)                                             ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCTGGtt           <  1:2503330/69‑1 (MQ=255)                                             ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCGGGtt           <  1:2499302/69‑1 (MQ=255)                                             ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCGGGtt           <  1:3202119/69‑1 (MQ=255)                                             ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCGGGtt           <  1:2886849/69‑1 (MQ=255)                                             ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCGGGtt           <  1:2171284/69‑1 (MQ=255)                                             ggTGCGGCGCACAATAGCCGGGCGATTCGCCGCAAGTGTCGGTTTATGCACATAAAACTCATCCGGGtt           <  1:275483/69‑1 (MQ=255)                                                      cACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCAcc  >  1:3121984/1‑69 (MQ=255)                                                      cACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCAcc  >  1:3033651/1‑69 (MQ=255)                                                      cACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCAcc  >  1:2008341/1‑69 (MQ=255)                                                      cACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCAcc  >  1:1120909/1‑69 (MQ=255)                                                         |                                                                  TCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACC  >  minE/1113983‑1114104

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