Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	98,996	C→A	45.5%	S271S (TCC→TCA)	guaC →	GMP reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	98,996	0	C	A	45.5%	35.3 / 98.9	97	S271S (TCC→TCA)	guaC	GMP reductase
Reads supporting (aligned to +/- strand):  ref base C (22/30);  new base A (10/34);  total (32/65)
Fisher's exact test for biased strand distribution p-value = 5.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

GGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACC  >  minE/98937‑99059                                                            |                                                                ggcggTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTc                                                            >  1:245159/1‑65 (MQ=255) ggcggTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGc                                                         >  1:1085548/1‑68 (MQ=255) ggcggTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGc                                                         >  1:2956266/1‑68 (MQ=255) ggcggTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCg                                                        >  1:2292923/1‑69 (MQ=255) ggcggTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCg                                                        >  1:1912499/1‑69 (MQ=255) ggcggTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCg                                                        >  1:2511851/1‑69 (MQ=255) ggcggTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCAGTTCTACGGCATGAGCTCCGAGTCTGc                                                         >  1:1570441/1‑68 (MQ=255)     gTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGa                                                    <  1:2656463/69‑1 (MQ=255)        gCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAAc                                                 >  1:903385/1‑69 (MQ=255)          aTCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGt                                               <  1:1250995/69‑1 (MQ=255)          aTCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGt                                               <  1:568481/69‑1 (MQ=255)          aTCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGt                                               <  1:1353499/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCTGAGTCTGCGATGAAACGTCACGTTGtcggcg                                  <  1:2073168/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2459481/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:997210/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2271561/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:90261/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2349462/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:73193/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2398394/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2413251/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:675148/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:305568/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:572248/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2561471/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:3214698/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2674761/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:3086416/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2743196/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2747356/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:3066029/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1192073/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1057168/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1630979/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1485230/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1642164/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1887167/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1801506/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1174041/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1398083/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1856119/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATGAGCGCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2451256/69‑1 (MQ=255)                       aCGGCGAGAAATTTATGCTGTTCTACGGCATAAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:867934/69‑1 (MQ=255)                       aCGGAGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:2812770/69‑1 (MQ=255)                         ggCGAGAAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:404911/67‑1 (MQ=255)                            gagaAATTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1209945/64‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTTCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:3409367/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGGGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:3273709/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:3145164/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:3077610/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:1388107/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:3250005/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:2886064/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:2843321/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:3328618/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:15687/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:3343772/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:2212791/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:192173/69‑1 (MQ=255)                             agaAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:1869691/69‑1 (MQ=255)                               aaaTTTATGCTGTTCTACGGCATGAGCTCAGAGTCTGCGATGAAACGTCACGTTGgcggcg                                  <  1:1914137/61‑1 (MQ=255)                                aaTTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAat                         <  1:3049926/69‑1 (MQ=255)                                               aCGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGaa           <  1:408819/68‑1 (MQ=255)                                               aCGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGaa           <  1:3158827/68‑1 (MQ=255)                                               aCGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGaa           <  1:3147571/68‑1 (MQ=255)                                               aCGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGaa           <  1:188291/68‑1 (MQ=255)                                               aCGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGaa           <  1:2151421/68‑1 (MQ=255)                                               aCGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGaa           <  1:719844/68‑1 (MQ=255)                                               aCGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGaa           <  1:3311151/68‑1 (MQ=255)                                               aCGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGaa           <  1:897908/68‑1 (MQ=255)                                               aCGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGaa           <  1:336082/68‑1 (MQ=255)                                                cGGCATGAGCTCCGAGTCTGCGATGAAACGTCAcgtt‑gcggcgttgcgg                            <  1:2530646/49‑1 (MQ=255)                                                     tGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCgg                            <  1:2969533/45‑1 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:3294708/1‑69 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:415200/1‑69 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:714112/1‑69 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:2984974/1‑69 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:1562010/1‑69 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:2724303/1‑69 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:2711461/1‑69 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:2379255/1‑69 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:218912/1‑69 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAc   >  1:880142/1‑68 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAc   >  1:1464792/1‑68 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAc   >  1:544949/1‑68 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAc   >  1:2594447/1‑68 (MQ=255)                                                       aGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGAAAAAcc  >  1:1422411/1‑69 (MQ=255)                                                       aGCTCAGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:1345902/1‑69 (MQ=255)                                                       aGCTCAGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:1250439/1‑69 (MQ=255)                                                       aGCTCAGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:3253848/1‑69 (MQ=255)                                                       aGCTCAGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:3047401/1‑69 (MQ=255)                                                       aGCTCAGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:2871659/1‑69 (MQ=255)                                                       aGCTCAGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:2842205/1‑69 (MQ=255)                                                       aGCTCAGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:1640362/1‑69 (MQ=255)                                                       aGCTCAGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:2417870/1‑69 (MQ=255)                                                       aGCTCAGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:2070757/1‑69 (MQ=255)                                                       aGCTCAGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAAcc  >  1:1026509/1‑69 (MQ=255)                                                            |                                                                GGCGGTCGCATCGTTGAGGAGAACGGCGAGAAATTTATGCTGTTCTACGGCATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACC  >  minE/98937‑99059

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