Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,086,868	G→A	52.3%	E215K (GAA→AAA)	yqjG →	predicted S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,086,868	0	G	A	52.3%	0.7 / 76.7	86	E215K (GAA→AAA)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (25/16);  new base A (36/9);  total (61/25)
Fisher's exact test for biased strand distribution p-value = 6.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

AAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAG  >  minE/2086812‑2086927                                                         |                                                            aaGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCAc                                                 >  1:1856614/1‑69 (MQ=255) aaGCTGGTTTTGCCA‑CAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACg                                                >  1:500080/1‑69 (MQ=255) aaGCTGGTTTTGCCA‑CAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACg                                                >  1:496257/1‑69 (MQ=255)                  agccagGAAGCTTACGACGAGGCTGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                <  1:297897/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                <  1:2390070/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                <  1:2595925/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                <  1:2174429/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                <  1:799574/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                <  1:593029/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                <  1:3388054/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                <  1:3388683/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:1092562/1‑69 (MQ=255)                         aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                         >  1:254596/1‑69 (MQ=255)                         aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                         >  1:876641/1‑69 (MQ=255)                         aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                         >  1:1836495/1‑69 (MQ=255)                         aaGCCTACGACGAGGCGTTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:3220395/1‑49 (MQ=37)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAAt                                                          >  1:3104729/1‑36 (MQ=37)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:3091027/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:745117/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:2869391/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:3104002/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:855023/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:2601000/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:3180964/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:3349985/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:2517461/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:496504/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:1046140/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:1667206/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:1462278/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:1788599/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:1190831/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:2226783/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                             >  1:120396/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGAGCGACt                                             >  1:1559154/1‑49 (MQ=37)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCt                                                      >  1:488466/1‑36 (MQ=255)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                <  1:1376793/58‑1 (MQ=255)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                      >  1:2482359/1‑68 (MQ=255)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                      >  1:2263092/1‑68 (MQ=255)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                      >  1:2268366/1‑68 (MQ=255)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                      >  1:2570840/1‑68 (MQ=255)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                      >  1:358267/1‑68 (MQ=255)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                      >  1:1800903/1‑68 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGa                                               >  1:549510/1‑39 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                  <  1:1824154/68‑1 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                  <  1:2038609/68‑1 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                  <  1:2107965/68‑1 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                 >  1:1547892/1‑69 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                 >  1:724597/1‑69 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                 >  1:1551505/1‑69 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                 >  1:1504772/1‑69 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                 >  1:123991/1‑69 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                 >  1:3240772/1‑69 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                 >  1:1692510/1‑69 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                 >  1:1800644/1‑69 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc        <  1:2391327/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc        <  1:1393137/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc        <  1:24715/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc        <  1:2239759/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc        <  1:2574738/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc        <  1:366354/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc        <  1:639723/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc        <  1:337538/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc        <  1:783654/69‑1 (MQ=255)                                            ggCGAAAGTGTTTGAATCGCTGGCACGACTCGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCaa      <  1:3276427/69‑1 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCa   >  1:845534/1‑67 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:1118078/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:1010921/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:351842/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:3404085/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:3357949/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:3352757/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:1766405/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:3072137/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:180907/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:2595199/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:2588373/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg  >  1:2541444/1‑68 (MQ=255)                                                         |                                                            AAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAG  >  minE/2086812‑2086927

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