Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	511,079	A→G	59.2%	A374A (GCA→GCG)	ybhJ →	predicted hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	511,079	0	A	G	59.2%	16.1 / 27.7	44	A374A (GCA→GCG)	ybhJ	predicted hydratase
Reads supporting (aligned to +/- strand):  ref base A (16/2);  new base G (16/10);  total (32/12)
Fisher's exact test for biased strand distribution p-value = 8.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATTATCGCGGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAATGACACCTTCTCGCTGGCAGTTTA  >  minE/511021‑511126                                                           |                                                ggATTATCGCGGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTc                                       >  1:1825681/1‑69 (MQ=255)  gATTATCGCGGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCa                                      >  1:3325250/1‑69 (MQ=255)      aTCGCGGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATcc                                  >  1:1126020/1‑69 (MQ=255)        cgcgGGTTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCt                                 <  1:3173644/68‑1 (MQ=255)        cgcgGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCt                                 <  1:1140566/68‑1 (MQ=255)        cgcgGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCt                                 <  1:858794/68‑1 (MQ=255)        cgcgGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCt                                 <  1:657873/68‑1 (MQ=255)        cgcgGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCt                                 <  1:44967/68‑1 (MQ=255)        cgcgGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCt                                 <  1:3397393/68‑1 (MQ=255)        cgcgGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCt                                 <  1:2844764/68‑1 (MQ=255)        cgcgGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCt                                 <  1:2645629/68‑1 (MQ=255)        cgcgGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCt                                 <  1:1995694/68‑1 (MQ=255)        cgcgGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCt                                 <  1:1654052/68‑1 (MQ=255)                   cTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAAt                                    >  1:2992344/1‑54 (MQ=255)                   cTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAAt                                    >  1:759128/1‑54 (MQ=255)                   cTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAAt                                    >  1:1519448/1‑54 (MQ=255)                   cTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAAt                                    >  1:1897765/1‑54 (MQ=255)                   cTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAAt                         >  1:250507/1‑65 (MQ=255)                   cTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAAt                         >  1:2837845/1‑65 (MQ=255)                   cTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAAt                         >  1:2403893/1‑65 (MQ=255)                   cTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAAt                         >  1:3404110/1‑65 (MQ=255)                   cTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAATGaca                     >  1:1813642/1‑69 (MQ=255)                        ggTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAATGACACCtt                 <  1:2777284/68‑1 (MQ=255)                        ggTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAATGACACCtt                 <  1:1723267/68‑1 (MQ=255)                         gTAACTACGAAAACGTCATCGCTGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:960228/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCTGCGAATGCGCTGCGCGGTCAAACCTGCGGCAATGACACCTtct               >  1:2760552/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:92715/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:2871390/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:1982473/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:320880/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:1697939/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:3341571/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:1693959/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:2866449/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:3418469/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:356875/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:1750338/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:541914/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:627792/1‑69 (MQ=255)                         gTAACTACGAAAACGTCATCGCGGCGGCGAATGCGCTGCGCGGTCAATCCTGCGGCAATGACACCTtct               >  1:2844484/1‑69 (MQ=255)                           aaCTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAATGACACCTTCTCg             >  1:3182841/1‑69 (MQ=255)                           aaCTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAATGACACCTTCTCg             >  1:1465771/1‑69 (MQ=255)                                       gTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAATGACACCTTCTCGCTGGCAGTTTa  >  1:2877490/1‑68 (MQ=255)                                       gTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAATGACACCTTCTCGCTGGCAGTTTa  >  1:1987993/1‑68 (MQ=255)                                                           |                                                GGATTATCGCGGGCTGTTCTGGCGGTAACTACGAAAACGTCATCGCGGCGGCGAATGCACTGCGCGGTCAATCCTGTGGCAATGACACCTTCTCGCTGGCAGTTTA  >  minE/511021‑511126

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