Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,639,640	A→G	60.1%	A133A (GCT→GCC)	ypfH ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,639,640	0	A	G	60.1%	31.0 / 46.8	70	A133A (GCT→GCC)	ypfH	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (7/21);  new base G (22/20);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 2.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCC  >  minE/1639576‑1639702                                                                 |                                                               gAATCAAGTGACTTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:3104448/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGTCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:1616058/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:98706/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:3208103/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:707741/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:2554656/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:2278969/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:218441/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:2101012/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:2095870/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:717826/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:87440/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:163753/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:1559821/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:1052038/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:1089191/69‑1 (MQ=255) gAATCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:985144/69‑1 (MQ=255)  aaTCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:770605/68‑1 (MQ=255)  aaTCAAGTGAATTGTTGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:1930508/68‑1 (MQ=255)          gAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGc                                                              <  1:619585/58‑1 (MQ=255)          gAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGc                                                              <  1:754930/58‑1 (MQ=255)          gAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGc                                                              <  1:1887950/58‑1 (MQ=255)          gAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGc                                                              <  1:629067/58‑1 (MQ=255)          gAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGc                                                              <  1:1631884/58‑1 (MQ=255)          gAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGc                                                              <  1:2941823/58‑1 (MQ=255)          gAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGc                                                              <  1:1218332/58‑1 (MQ=255)           aaTTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGc                                                              <  1:1001633/57‑1 (MQ=255)              tGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGc                                               <  1:2847773/69‑1 (MQ=255)                     gCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGccc                                        >  1:284923/1‑69 (MQ=255)                           gACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAGGCaa                                                            <  1:1498108/43‑1 (MQ=255)                                    tCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATg                         <  1:3271286/69‑1 (MQ=255)                                    tCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATg                         >  1:2616148/1‑69 (MQ=255)                                    tCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATg                         <  1:627674/69‑1 (MQ=255)                                    tCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATg                         <  1:1943890/69‑1 (MQ=255)                                                   tAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa          <  1:3270919/69‑1 (MQ=255)                                                   tAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa          <  1:2896188/69‑1 (MQ=255)                                                   tAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa          <  1:2895963/69‑1 (MQ=255)                                                   tAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa          <  1:346286/69‑1 (MQ=255)                                                   tAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa          <  1:2760691/69‑1 (MQ=255)                                                   tAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa          <  1:229065/69‑1 (MQ=255)                                                   tAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa          <  1:11085/69‑1 (MQ=255)                                                   tAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa          <  1:742006/69‑1 (MQ=255)                                                   tAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa          <  1:819227/69‑1 (MQ=255)                                                     gCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATaa         >  1:3203415/1‑68 (MQ=255)                                                     gCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAAt        >  1:263708/1‑69 (MQ=255)                                                     gCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAAt        >  1:3039347/1‑69 (MQ=255)                                                     gCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAAt        >  1:3315591/1‑69 (MQ=255)                                                     gCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAAt        >  1:526054/1‑69 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCATCGCCATAATCGCGcc  >  1:1047946/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCgcgc   >  1:46990/1‑67 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:649611/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:2827321/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:708863/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:950482/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:90058/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:2040516/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:1033027/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:1422416/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:1662318/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:1707283/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:1737069/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:1745458/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:491806/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:2074712/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:2119961/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:2571231/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:2825548/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:3305928/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:3324512/1‑68 (MQ=255)                                                            ttAAAGGCAATGACGCGTGAGGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCGCCATAATCGCGcc  >  1:3400380/1‑68 (MQ=255)                                                              aaaaGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGcc  >  1:2783087/1‑66 (MQ=255)                                                              aaaaGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGcc  >  1:834926/1‑66 (MQ=255)                                                              aaaaGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGcc  >  1:1517284/1‑66 (MQ=255)                                                              aaaaGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGcc  >  1:1301349/1‑66 (MQ=255)                                                                 |                                                               GAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCC  >  minE/1639576‑1639702

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