Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,086,879	A→G	47.3%	A218A (GCA→GCG)	yqjG →	predicted S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,086,879	0	A	G	47.3%	2.7 / 82.5	93	A218A (GCA→GCG)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base G (35/9);  total (59/34)
Fisher's exact test for biased strand distribution p-value = 2.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

AAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGA  >  minE/2086812‑2086941                                                                    |                                                               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:593029/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:2174429/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:2390070/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:3388054/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:799574/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:3388683/69‑1 (MQ=255)                  agccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:2595925/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:254596/1‑69 (MQ=255)                         aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:876641/1‑69 (MQ=255)                         aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:1836495/1‑69 (MQ=255)                         aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:591809/1‑69 (MQ=255)                         aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:1092562/1‑69 (MQ=255)                         aaGCCTACGACGAGGCGTTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:3220395/1‑49 (MQ=37)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:1046140/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:2601000/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:2869391/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:3091027/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:496504/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:2517461/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:3104002/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:3180964/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:2226783/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:3349985/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:1190831/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:120396/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:855023/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:1462278/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:745117/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:1667206/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGCGCGACt                                                           >  1:1788599/1‑49 (MQ=38)                         aaGCCTACGACGAGGCGGTGGCGAAAGTGTTTAAATCGCTGGAGCGACt                                                           >  1:1559154/1‑49 (MQ=37)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTa                                              <  1:1376793/58‑1 (MQ=255)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:2570840/1‑68 (MQ=255)                             ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:358267/1‑68 (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:1800903/1‑68 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGa                                                             >  1:549510/1‑39 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                                <  1:2038609/68‑1 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                                <  1:2107965/68‑1 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                                <  1:1824154/68‑1 (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:1692510/1‑69 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:724597/1‑69 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:123991/1‑69 (MQ=255)                                 gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1800644/1‑69 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:2574738/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:337538/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:366354/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:783654/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:639723/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:1393137/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:2239759/69‑1 (MQ=255)                                          gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:2391327/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:3357949/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1118078/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1327848/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1766405/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:180907/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:519403/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1010921/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3352757/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3072137/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3404085/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:351842/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:2595199/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:2588373/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:366492/1‑68 (MQ=255)                                                 aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:2541444/1‑68 (MQ=255)                                                          aaTCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGaa       <  1:2104850/68‑1 (MQ=255)                                                          aaTCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGaa       <  1:851755/68‑1 (MQ=255)                                                          aaTCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGaa       <  1:930283/68‑1 (MQ=255)                                                              gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:2030927/69‑1 (MQ=255)                                                              gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:1931762/69‑1 (MQ=255)                                                              gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:1394879/69‑1 (MQ=255)                                                              gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:473901/69‑1 (MQ=255)                                                              gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:1309153/69‑1 (MQ=255)                                                              gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCCGCAACCAGCTAACCGAAGCCGa  <  1:3099656/69‑1 (MQ=255)                                                                    |                                                               AAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGA  >  minE/2086812‑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: ‑