Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,086,888	G→A	35.6%	E221E (GAG→GAA)	yqjG →	predicted S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,086,888	0	G	A	35.6%	48.3 / 38.1	70	E221E (GAG→GAA)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (20/25);  new base A (16/9);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGA  >  minE/2086829‑2086941                                                            |                                                      agccagGAAGCTTACGACGAGGCTGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:297897/69‑1 (MQ=255) agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:3388054/69‑1 (MQ=255) agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:2174429/69‑1 (MQ=255) agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:2390070/69‑1 (MQ=255) agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:2595925/69‑1 (MQ=255) agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:3388683/69‑1 (MQ=255) agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:593029/69‑1 (MQ=255) agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:799574/69‑1 (MQ=255) agccagGAAGCTTACGACGAGGCGGTGGCGAAACTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:92310/69‑1 (MQ=255)  gccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:2274246/68‑1 (MQ=255)  gccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAg                                             <  1:872363/69‑1 (MQ=255)        aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAg                                        >  1:2862762/1‑68 (MQ=255)        aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:591809/1‑69 (MQ=255)        aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:876641/1‑69 (MQ=255)        aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:254596/1‑69 (MQ=255)        aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:1836495/1‑69 (MQ=255)        aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:1092562/1‑69 (MQ=255)            ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:1376793/58‑1 (MQ=255)            ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:358267/1‑68 (MQ=255)            ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:2570840/1‑68 (MQ=255)            ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:2482359/1‑68 (MQ=255)            ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:2268366/1‑68 (MQ=255)            ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:2263092/1‑68 (MQ=255)            ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:1800903/1‑68 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                                <  1:2107965/68‑1 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                                <  1:2038609/68‑1 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                                <  1:1824154/68‑1 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1800644/1‑69 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1692510/1‑69 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:3240772/1‑69 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1551505/1‑69 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1547892/1‑69 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1504772/1‑69 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:123991/1‑69 (MQ=255)                gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:724597/1‑69 (MQ=255)                         gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:366354/69‑1 (MQ=255)                         gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:1393137/69‑1 (MQ=255)                         gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:337538/69‑1 (MQ=255)                         gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:639723/69‑1 (MQ=255)                         gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:783654/69‑1 (MQ=255)                         gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:2239759/69‑1 (MQ=255)                         gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:2391327/69‑1 (MQ=255)                         gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:2574738/69‑1 (MQ=255)                         gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:24715/69‑1 (MQ=255)                           ggCGAAAGTGTTTGAATCGCTGGCACGACTCGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCaa                    <  1:3276427/69‑1 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCa                 >  1:845534/1‑67 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:2595199/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1118078/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1327848/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:519403/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:366492/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:180907/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:351842/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3404085/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:2541444/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1010921/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3357949/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3352757/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:2588373/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3072137/1‑68 (MQ=255)                                aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1766405/1‑68 (MQ=255)                                         aaTCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGaa       <  1:851755/68‑1 (MQ=255)                                         aaTCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGaa       <  1:930283/68‑1 (MQ=255)                                         aaTCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGaa       <  1:2104850/68‑1 (MQ=255)                                             gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:1394879/69‑1 (MQ=255)                                             gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:473901/69‑1 (MQ=255)                                             gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:1309153/69‑1 (MQ=255)                                             gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:1931762/69‑1 (MQ=255)                                             gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:2030927/69‑1 (MQ=255)                                             gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCCGCAACCAGCTAACCGAAGCCGa  <  1:3099656/69‑1 (MQ=255)                                                            |                                                      AGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGA  >  minE/2086829‑2086941

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