Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	252,556	A→G	36.9%	S525G (AGT→GGT)	malZ →	maltodextrin glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	252,556	0	A	G	36.9%	39.5 / 44.5	73	S525G (AGT→GGT)	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base A (30/16);  new base G (19/8);  total (49/24)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATA  >  minE/252494‑252604                                                               |                                                 tGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAgg                                            <  1:959138/69‑1 (MQ=255) tGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAgg                                            <  1:886743/69‑1 (MQ=255) tGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAgg                                            <  1:3406420/69‑1 (MQ=255) tGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAgg                                            <  1:1340136/69‑1 (MQ=255) tGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAgg                                            <  1:179356/69‑1 (MQ=255) tGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAgg                                            <  1:1526477/69‑1 (MQ=255) tGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAgg                                            <  1:1460238/69‑1 (MQ=255)         cAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGt                                    <  1:46354/69‑1 (MQ=255)         cAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGt                                    <  1:3199187/69‑1 (MQ=255)         cAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGt                                    <  1:244703/69‑1 (MQ=255)              tACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATgg                               >  1:2490856/1‑69 (MQ=255)              tACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATgg                               >  1:2576878/1‑69 (MQ=255)              tACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATgg                               >  1:1562220/1‑69 (MQ=255)              tACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATgg                               >  1:1398658/1‑69 (MQ=255)                cGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                             >  1:2485216/1‑69 (MQ=255)                cGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                             <  1:299191/69‑1 (MQ=255)                cGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                             >  1:3017812/1‑69 (MQ=255)                cGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcg                             >  1:3154953/1‑69 (MQ=255)                  gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                           >  1:1125018/1‑69 (MQ=255)                  gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                           >  1:3117341/1‑69 (MQ=255)                  gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                           >  1:2010498/1‑69 (MQ=255)                  gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                           >  1:87045/1‑69 (MQ=255)                  gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                           >  1:832853/1‑69 (MQ=255)                  gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                           >  1:1367298/1‑69 (MQ=255)                  gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                           >  1:1711662/1‑69 (MQ=255)                  gCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcggc                           >  1:464509/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:3085536/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:649614/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:3171573/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:331678/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:2954152/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:1219385/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:282264/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:2804174/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:2726916/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:1212644/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:1311457/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:1428314/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:1529389/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:1940176/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:2282454/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:2394072/1‑68 (MQ=255)                    gTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCt                          >  1:2442256/1‑68 (MQ=255)                          cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                    <  1:2697389/68‑1 (MQ=255)                          cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                    <  1:136173/68‑1 (MQ=255)                          cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                    <  1:2760026/68‑1 (MQ=255)                          cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                    <  1:2165796/68‑1 (MQ=255)                          cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                    <  1:3281756/68‑1 (MQ=255)                          cgcgCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCATCGTGGCGGCTGTCAgg                    <  1:1037326/68‑1 (MQ=255)                                 tACCAGCGAATGATTGCGCTGCGTAAGAAAGGTCAGGCGTTGCGTCGTGGCGGCTGTCAgg                    <  1:911256/61‑1 (MQ=255)                                          aTGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGAt   <  1:1239780/69‑1 (MQ=255)                                          aTGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGAt   <  1:1136010/69‑1 (MQ=255)                                          aTGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGAt   <  1:1928150/69‑1 (MQ=255)                                          aTGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGAt   <  1:1350938/69‑1 (MQ=255)                                          aTGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGAt   <  1:2449032/69‑1 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCt                >  1:3360673/1‑55 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTAt            >  1:622308/1‑59 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa  >  1:1437742/1‑69 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa  >  1:1381314/1‑69 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa  >  1:338559/1‑69 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa  >  1:737284/1‑69 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa  >  1:638694/1‑69 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa  >  1:136085/1‑69 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa  >  1:2002061/1‑69 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa  >  1:2125692/1‑69 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa  >  1:3022539/1‑69 (MQ=255)                                           tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATa  >  1:2878361/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)                                                               |                                                 TGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATA  >  minE/252494‑252604

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