Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	252,565	C→T	60.2%	L528L (CTA→TTA) ‡	malZ →	maltodextrin glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	252,565	0	C	T	60.2%	41.7 / 69.0	98	L528L (CTA→TTA) ‡	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base C (30/9);  new base T (19/40);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 2.67e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

CAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGTGCTGA  >  minE/252502‑252631                                                                |                                                                   cAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGt                                                               <  1:244703/69‑1 (MQ=255) cAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGt                                                               <  1:46354/69‑1 (MQ=255) cAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGt                                                               <  1:3199187/69‑1 (MQ=255)      tACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATgg                                                          >  1:1562220/1‑69 (MQ=255)      tACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATgg                                                          >  1:2576878/1‑69 (MQ=255)      tACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATgg                                                          >  1:2490856/1‑69 (MQ=255)      tACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATgg                                                          >  1:1398658/1‑69 (MQ=255)        cGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                                                        >  1:3154953/1‑69 (MQ=255)        cGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                                                        <  1:299191/69‑1 (MQ=255)        cGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                                                        >  1:3017812/1‑69 (MQ=255)        cGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                                                        >  1:2485216/1‑69 (MQ=255)          gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                                                      >  1:1711662/1‑69 (MQ=255)          gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                                                      >  1:832853/1‑69 (MQ=255)          gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                                                      >  1:87045/1‑69 (MQ=255)          gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                                                      >  1:1367298/1‑69 (MQ=255)          gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                                                      >  1:464509/1‑69 (MQ=255)          gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                                                      >  1:2010498/1‑69 (MQ=255)          gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                                                      >  1:3117341/1‑69 (MQ=255)          gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                                                      >  1:1125018/1‑69 (MQ=255)          gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcgg                                                       >  1:3041790/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTggcggc                                                      >  1:3249120/1‑67 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTgg                                                          >  1:263061/1‑63 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:2442256/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:2726916/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:331678/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:649614/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:2282454/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:2394072/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:1940176/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:3171573/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:2804174/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:282264/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:1212644/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:1529389/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:2954152/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:1428314/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:1219385/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:1311457/1‑68 (MQ=255)            gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                                                     >  1:3085536/1‑68 (MQ=255)                  cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                                               <  1:2760026/68‑1 (MQ=255)                  cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                                               <  1:2165796/68‑1 (MQ=255)                  cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                                               <  1:136173/68‑1 (MQ=255)                  cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                                               <  1:3281756/68‑1 (MQ=255)                  cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                                               <  1:2697389/68‑1 (MQ=255)                  cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCATCGTGGCGGCTGTCAgg                                               <  1:1037326/68‑1 (MQ=255)                         tACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                                               <  1:911256/61‑1 (MQ=255)                                  aTGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGAt                              <  1:1350938/69‑1 (MQ=255)                                  aTGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGAt                              <  1:2449032/69‑1 (MQ=255)                                  aTGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGAt                              <  1:1136010/69‑1 (MQ=255)                                  aTGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGAt                              <  1:1239780/69‑1 (MQ=255)                                  aTGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGAt                              <  1:1928150/69‑1 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCt                                           >  1:3360673/1‑55 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTAt                                       >  1:622308/1‑59 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:2002061/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:3022539/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:136085/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:1381314/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:1437742/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:2878361/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:737284/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:638694/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:2125692/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:338559/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGATGTCAGGTGCTGTATGCGGAAGATa                             >  1:174774/1‑69 (MQ=255)                                   tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGAGGCTGTCAGGTGCTGTATGCGGAAGATa                             >  1:2068317/1‑69 (MQ=255)                                        gcgcTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                                               <  1:745542/46‑1 (MQ=37)                                                     aaggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATG‑GGAAGATAACGTGGTGGTATTTGTCcg          <  1:1508607/66‑1 (MQ=255)                                                      aggTCAGGGGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:1654220/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTTTCcg          <  1:2331030/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:915398/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:880735/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:700520/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:3318480/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:3371777/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:386415/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2271610/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:1540640/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:1823731/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:1978800/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2072933/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2094012/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:211663/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2125562/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:3078997/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2129619/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2410734/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2411073/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:247696/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2496216/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2553904/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2589292/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:1350901/67‑1 (MQ=255)                                                      aggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:1077937/67‑1 (MQ=255)                                                       ggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:506328/67‑1 (MQ=255)                                                       ggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:1434395/67‑1 (MQ=255)                                                       ggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:1176493/67‑1 (MQ=255)                                                       ggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:3299765/67‑1 (MQ=255)                                                       ggTCAGGCGTTGCGTCGTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCcg          <  1:2829291/67‑1 (MQ=255)                                                              cgcTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGTGCTGa  <  1:2995311/69‑1 (MQ=255)                                                               gcTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTAATTGTCCGCGTGCTGa  <  1:642084/68‑1 (MQ=255)                                                                |                                                                   CAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGTGCTGA  >  minE/252502‑252631

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑