Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,492,955	A→G	18.9%	V584V (GTA→GTG)	gyrB →	DNA gyrase, subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,492,955	0	A	G	18.9%	116.2 / 17.8	79	V584V (GTA→GTG)	gyrB	DNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base A (46/18);  new base G (7/8);  total (53/26)
Fisher's exact test for biased strand distribution p-value = 7.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGCAATGCTGA  >  minE/2492887‑2493022                                                                     |                                                                    gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:1379139/1‑69 (MQ=255) gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:1582637/1‑69 (MQ=255) gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:3189671/1‑69 (MQ=255) gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:1404583/1‑69 (MQ=255) gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:1400622/1‑69 (MQ=255) gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:3122200/1‑69 (MQ=255) gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:725153/1‑69 (MQ=255) gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:2223759/1‑69 (MQ=255) gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:739642/1‑69 (MQ=255) gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:97016/1‑69 (MQ=255) gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                                                                     >  1:267595/1‑69 (MQ=255)   cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                                                                    >  1:66919/1‑68 (MQ=255)   cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                                                                    >  1:2914662/1‑68 (MQ=255)   cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                                                                    >  1:2793041/1‑68 (MQ=255)   cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                                                                    >  1:3183530/1‑68 (MQ=255)   cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                                                                    >  1:3211778/1‑68 (MQ=255)   cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                                                                    >  1:1631405/1‑68 (MQ=255)               gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc                                                       <  1:981838/69‑1 (MQ=255)               gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc                                                       <  1:2367122/69‑1 (MQ=255)               gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc                                                       <  1:3286607/69‑1 (MQ=255)               gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc                                                       <  1:3322203/69‑1 (MQ=255)               gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc                                                       <  1:374868/69‑1 (MQ=255)               gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc                                                       <  1:350869/69‑1 (MQ=255)               gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc                                                       <  1:1159763/69‑1 (MQ=255)                cacaCCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc                                                       <  1:310988/68‑1 (MQ=255)                      aaCGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCa                                                 >  1:2905026/1‑68 (MQ=255)                      aaCGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCa                                                 >  1:1511550/1‑68 (MQ=255)                      aaCGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCa                                                 >  1:2515853/1‑68 (MQ=255)                      aaCGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCa                                                 >  1:26627/1‑68 (MQ=255)                        cGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGaa                                              <  1:2219407/69‑1 (MQ=255)                        cGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGaa                                              <  1:635604/69‑1 (MQ=255)                        cGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGaa                                              <  1:2164694/69‑1 (MQ=255)                        cGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGaa                                              <  1:2153411/69‑1 (MQ=255)                           cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGaaaa                                            >  1:3217499/1‑68 (MQ=255)                           cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt                                           >  1:1600848/1‑69 (MQ=255)                           cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt                                           >  1:2562448/1‑69 (MQ=255)                           cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt                                           >  1:2170665/1‑69 (MQ=255)                           cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt                                           >  1:2709993/1‑69 (MQ=255)                           cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt                                           >  1:1866227/1‑69 (MQ=255)                           cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt                                           >  1:1690240/1‑69 (MQ=255)                               gCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGc                                                  <  1:3024683/58‑1 (MQ=255)                               gCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGc                                                  <  1:2015828/58‑1 (MQ=255)                               gCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGc                                                  <  1:737614/58‑1 (MQ=255)                               gCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGc                                                  <  1:1177094/58‑1 (MQ=255)                                         tggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATgg                             <  1:3320736/69‑1 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACa                                                           >  1:2152950/1‑38 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCg                                  >  1:458341/1‑63 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCg                                  >  1:188018/1‑63 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCg                                  >  1:184047/1‑63 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCg                                  >  1:2539940/1‑63 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:3331534/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:2375930/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:335994/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:2721687/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:288109/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:2039557/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:1420806/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:1430905/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:1523148/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:1548322/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:156465/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:1727612/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:1868393/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:2026054/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:1355970/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:2149722/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:133637/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:1230898/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:2428365/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:2482266/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:1059903/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:2486126/1‑69 (MQ=255)                                          ggctggcGAAGCGTTAGAGAAACTGGTA‑‑TGAGTACAACGCGACGCAGAAAATGATCAATAGTATGGAGc                          >  1:80757/1‑69 (MQ=255)                                          ggctggaGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGa                            >  1:440774/1‑69 (MQ=255)                                              ggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGcgt                        >  1:1056597/1‑69 (MQ=255)                                              ggcGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGcgt                        >  1:1671523/1‑69 (MQ=255)                                                         agagAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGaa             <  1:649068/69‑1 (MQ=255)                                                         agagAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGaa             <  1:2732449/69‑1 (MQ=255)                                                         agagAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGaa             <  1:579165/69‑1 (MQ=255)                                                            gaAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          <  1:1631770/69‑1 (MQ=255)                                                            gaAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          >  1:1381196/1‑69 (MQ=255)                                                            gaAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          >  1:137101/1‑69 (MQ=255)                                                            gaAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          >  1:3026540/1‑69 (MQ=255)                                                            gaAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          <  1:3019241/69‑1 (MQ=255)                                                            gaAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          <  1:2217525/69‑1 (MQ=255)                                                            gaAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          >  1:2856244/1‑69 (MQ=255)                                                            gaAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          >  1:1112390/1‑69 (MQ=255)                                                            gaAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          <  1:2789986/69‑1 (MQ=255)                                                             aaaCTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          <  1:169355/68‑1 (MQ=255)                                                             aaaCTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGc          <  1:2045400/68‑1 (MQ=255)                                                                     aTCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGCAATGCTGa  <  1:1223161/68‑1 (MQ=255)                                                                     |                                                                    GACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAATGATCAATCGTATGGAGCGTCGTTATCCGAAAGCAATGCTGA  >  minE/2492887‑2493022

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