Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,759	A→G	10.1%	G110G (GGT→GGC)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,759	0	A	G	10.1%	233.2 / 13.9	90	G110G (GGT→GGC)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (35/45);  new base G (4/5);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

GAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTC  >  minE/313690‑313827                                                                      |                                                                     gAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCAc                                                                     <  1:4211574/71‑1 (MQ=255)       tttCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCg                                                               >  1:2412473/1‑71 (MQ=255)         tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGc                                                              <  1:5119620/70‑1 (MQ=255)         tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGc                                                              <  1:4201905/70‑1 (MQ=255)         tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                             >  1:553005/1‑71 (MQ=255)         tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                             >  1:1300063/1‑71 (MQ=255)         tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                             >  1:3953699/1‑71 (MQ=255)         tCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCa                                                             >  1:1599952/1‑71 (MQ=255)          cTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCAc                                                            >  1:5335662/1‑71 (MQ=255)           tGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACt                                                           >  1:2341343/1‑71 (MQ=255)                aCGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGt                                                      <  1:5612005/71‑1 (MQ=255)                    ttCACAGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGt                                                  >  1:4312247/1‑71 (MQ=255)                      cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                >  1:4604195/1‑71 (MQ=255)                       aCCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCAc                                                            >  1:4115660/1‑58 (MQ=255)                       aCCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGc                                               >  1:984623/1‑71 (MQ=255)                         cGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCt                                             <  1:3861191/71‑1 (MQ=255)                         cGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCt                                             <  1:4053658/71‑1 (MQ=255)                           tCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                           <  1:3068027/71‑1 (MQ=255)                             aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                           >  1:1642031/1‑69 (MQ=255)                              aTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAg                                        <  1:1018523/71‑1 (MQ=255)                                tGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                           >  1:4796092/1‑66 (MQ=255)                                      atatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                <  1:3929537/71‑1 (MQ=255)                                      atatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                <  1:4099394/71‑1 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                       <  1:1788707/47‑1 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCg                                                       <  1:3997587/47‑1 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGcc                                              >  1:2795501/1‑56 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTg                                            <  1:1569277/58‑1 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTg                                            <  1:2777007/58‑1 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGc                                   <  1:3847243/67‑1 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                <  1:4106053/70‑1 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                >  1:3301504/1‑70 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                               >  1:5692641/1‑71 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                               >  1:3319957/1‑71 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                               >  1:488604/1‑71 (MQ=255)                                       tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                               >  1:2168727/1‑71 (MQ=255)                                         ttGGCTGCAGCCTGGGCTTTCGCCAGATCCCCTTTCGCACTGTCGTa                                                     >  1:2059740/1‑47 (MQ=255)                                         ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTa                                                     >  1:5415865/1‑47 (MQ=255)                                         ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                    <  1:245898/48‑1 (MQ=255)                                         ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATg                                                   <  1:2865650/49‑1 (MQ=255)                                         ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGt                                                  <  1:972198/50‑1 (MQ=255)                                         ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                          >  1:2566945/1‑58 (MQ=255)                                         ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                          >  1:1616252/1‑58 (MQ=255)                                          tGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                    <  1:471907/47‑1 (MQ=255)                                          tGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                    <  1:3406594/47‑1 (MQ=255)                                          tGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTAt                                                    <  1:2366677/47‑1 (MQ=255)                                           ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                            <  1:2647068/70‑1 (MQ=255)                                           ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                            <  1:2673492/70‑1 (MQ=255)                                                  gCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAg                                 <  1:726102/58‑1 (MQ=255)                                                    cTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                    >  1:4880331/1‑69 (MQ=255)                                                     tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTa                                                     >  1:843830/1‑35 (MQ=255)                                                     tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTa                                                     >  1:1303687/1‑35 (MQ=255)                                                     tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                    >  1:2990771/1‑68 (MQ=255)                                                      gggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                          <  1:5272139/45‑1 (MQ=255)                                                      gggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                          <  1:5016708/45‑1 (MQ=255)                                                      gggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                          <  1:4795491/45‑1 (MQ=255)                                                      gggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                          <  1:2917242/45‑1 (MQ=255)                                                       ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagaga                <  1:3043379/70‑1 (MQ=255)                                                       ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                 <  1:3598000/69‑1 (MQ=255)                                                       ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                 <  1:1591286/69‑1 (MQ=255)                                                        gCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagaga                >  1:2471453/1‑69 (MQ=255)                                                         ttttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATcg                            <  1:4541670/55‑2 (MQ=255)                                                         cTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                           <  1:2323984/57‑1 (MQ=255)                                                         cTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                           <  1:1300814/57‑1 (MQ=255)                                                          tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                   <  1:4010275/64‑1 (MQ=255)                                                          tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                   <  1:5227791/64‑1 (MQ=255)                                                          tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                   <  1:3241212/64‑1 (MQ=255)                                                          tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                   <  1:2361173/64‑1 (MQ=255)                                                           ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGAt                          >  1:3474709/1‑56 (MQ=255)                                                           ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGAt                          >  1:4953910/1‑56 (MQ=255)                                                           ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggactg           >  1:1668823/1‑68 (MQ=255)                                                           ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggact            <  1:4523712/70‑3 (MQ=255)                                                           ttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                   <  1:121238/63‑1 (MQ=255)                                                            tCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTgg                      <  1:1020961/59‑2 (MQ=255)                                                            tCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggactg           <  1:1499681/70‑4 (MQ=255)                                                            tCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggactg           <  1:762063/70‑4 (MQ=255)                                                            tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGaca          >  1:3964632/1‑71 (MQ=255)                                                             cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                  <  1:3554642/62‑1 (MQ=255)                                                             cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGacac         >  1:1295142/1‑71 (MQ=255)                                                              gCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTa                    >  1:23704/1‑59 (MQ=255)                                                                 aGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGc     >  1:4880651/1‑71 (MQ=255)                                                                   aTCACCTTTCGCACTGTCGTATGTCGGCTGATAGGTCGCAGGATCAATTTGATAGAGAGAGACACCTGCt    <  1:2381469/70‑1 (MQ=255)                                                                   aTCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCg                                    <  1:2793249/38‑1 (MQ=255)                                                                   aTCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCg                                    <  1:4235960/38‑1 (MQ=255)                                                                   aTCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCtt   >  1:499912/1‑71 (MQ=255)                                                                    tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc  <  1:3853621/71‑1 (MQ=255)                                                                    tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc  <  1:5018371/71‑1 (MQ=255)                                                                     cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc  >  1:3260999/1‑70 (MQ=255)                                                                     cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc  >  1:3171989/1‑70 (MQ=255)                                                                     cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc  >  1:67833/1‑70 (MQ=255)                                                                     cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc  >  1:2200629/1‑70 (MQ=255)                                                                      |                                                                     GAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTC  >  minE/313690‑313827

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑