Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 313,801 A→G 7.2% I96I (ATT→ATC acrA ← multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE313,8010AG 7.2% 238.7 / 8.6 84I96I (ATT→ATCacrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (41/36);  new base G (3/3);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCG  >  minE/313732‑313870
                                                                     |                                                                     
ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                                                      >  1:3447952/1‑71 (MQ=255)
ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                       <  1:2953146/70‑1 (MQ=255)
ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                       <  1:3812138/70‑1 (MQ=255)
     cAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                                                 <  1:2514709/71‑1 (MQ=255)
        ccTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                                                 <  1:3646586/68‑1 (MQ=255)
         cTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                                               >  1:2684716/1‑69 (MQ=255)
          tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                       <  1:999900/60‑1 (MQ=255)
          tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                       <  1:430368/60‑1 (MQ=255)
          tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCATTCTGATa                                                               >  1:1592507/1‑68 (MQ=255)
          tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGAt                                                                >  1:575804/1‑67 (MQ=255)
            ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                     >  1:3178778/1‑60 (MQ=255)
            ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                                     >  1:3545509/1‑60 (MQ=255)
            ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                             >  1:1097008/1‑68 (MQ=255)
             gCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagagaga                                                         <  1:2717014/71‑1 (MQ=255)
               tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                                            <  1:3297119/66‑1 (MQ=255)
               tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                                            <  1:1907769/66‑1 (MQ=255)
                ttCGCCAGATCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggactg                                                      >  1:1811369/1‑68 (MQ=255)
                ttCGCCAGATCACCTTTCGCACTGTCGTATGTCCCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                                       >  1:2040572/1‑70 (MQ=255)
                 tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                                       <  1:316824/69‑1 (MQ=255)
                 tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                                       <  1:2229540/69‑1 (MQ=255)
                 tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                                                      <  1:1184254/70‑1 (MQ=255)
                 tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                                                      <  1:3670336/70‑1 (MQ=255)
                 tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                                                      <  1:1322179/70‑1 (MQ=255)
                   gCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                              <  1:4120862/60‑1 (MQ=255)
                     cAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                                                 <  1:1232969/55‑1 (MQ=255)
                     cAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                                       >  1:3914454/1‑65 (MQ=255)
                     cAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCt                                                  >  1:514151/1‑70 (MQ=255)
                       gATCACCTTTCGCACTGTCGTATGTCGCCTTATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCt                                               <  1:441301/71‑1 (MQ=255)
                       gATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCt                                               <  1:2203957/71‑1 (MQ=255)
                        aTCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCt                                                                   <  1:3090864/50‑1 (MQ=255)
                         tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                             <  1:2338972/71‑1 (MQ=255)
                         tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                             <  1:536999/71‑1 (MQ=255)
                          cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                              >  1:775562/1‑53 (MQ=255)
                          cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                             >  1:713602/1‑70 (MQ=255)
                            ccTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCg                                            <  1:2225833/69‑1 (MQ=255)
                             cTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATg                                         <  1:2848394/71‑1 (MQ=255)
                               ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                                                  <  1:2144366/44‑1 (MQ=255)
                               ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                                                  <  1:1586394/44‑1 (MQ=255)
                               ttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                                                  <  1:3671282/44‑1 (MQ=255)
                               ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTc                                       >  1:1162506/1‑71 (MQ=255)
                               ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTc                                       >  1:1552036/1‑71 (MQ=255)
                                tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTa                                                               >  1:1278226/1‑46 (MQ=255)
                                 cgcgCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTggt                                                                >  1:1986979/1‑42 (MQ=255)
                                   cACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                                                      <  1:2013728/36‑1 (MQ=255)
                                   cACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                                                      <  1:2159283/36‑1 (MQ=255)
                                    aCTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCt                                                  >  1:1316272/1‑55 (MQ=255)
                                    aCTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCt                                                  >  1:3968059/1‑55 (MQ=255)
                                     cTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATc                                                                    >  1:4027660/1‑36 (MQ=255)
                                     cTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTAcc                                 >  1:2066732/1‑71 (MQ=255)
                                     cTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTAc                                  >  1:3722110/1‑70 (MQ=255)
                                      tGTCGTATGTCGCCTGATAGGTCGCAGGATCTATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTa                                   <  1:428906/68‑1 (MQ=255)
                                      tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCg                                            <  1:2301160/59‑1 (MQ=255)
                                      tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTa                                   <  1:1771386/68‑1 (MQ=255)
                                        tCGTATGTCGCCTGATAGGTCTCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                             <  1:1791863/56‑1 (MQ=255)
                                        tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTa                                   >  1:3876393/1‑66 (MQ=255)
                                        tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttc                              >  1:1580177/1‑71 (MQ=255)
                                          gTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                                               >  1:2030076/1‑36 (MQ=255)
                                          gTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttctt                            >  1:1450486/1‑71 (MQ=255)
                                           tATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTTGCTACCTTCttt                           >  1:1451079/1‑71 (MQ=255)
                                           tATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt                           <  1:2620215/71‑1 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATagagagag                                                        >  1:1421552/1‑40 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATagaga                                                           >  1:858950/1‑37 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGc                                     >  1:2204003/1‑59 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTg                          >  1:1325165/1‑70 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                         <  1:1259639/71‑1 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                         <  1:2436995/71‑1 (MQ=255)
                                               tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt                           <  1:434208/67‑1 (MQ=255)
                                                cGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt                                        >  1:4047242/1‑53 (MQ=255)
                                                 gCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGc                                     >  1:1688524/1‑55 (MQ=255)
                                                  ccTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGa                                           >  1:2331952/1‑48 (MQ=255)
                                                    tGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                                                      <  1:2989797/35‑1 (MQ=255)
                                                     gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGc                 >  1:86362/1‑71 (MQ=255)
                                                      aTAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                         <  1:2817295/62‑1 (MQ=255)
                                                      aTAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTa                    >  1:2662533/1‑67 (MQ=255)
                                                       tAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCt                >  1:3936640/1‑70 (MQ=255)
                                                         ggTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTc              <  1:225172/70‑1 (MQ=255)
                                                         ggTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCa             <  1:3618414/71‑1 (MQ=255)
                                                             gCAGGATCAATCTGATAGAGAGAGACACCTTCTTCGAt                                          <  1:1585523/38‑1 (MQ=255)
                                                             gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGAt                                          <  1:2488304/38‑1 (MQ=255)
                                                                ggATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAAt      >  1:3438701/1‑71 (MQ=255)
                                                                ggATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAAt      >  1:272985/1‑71 (MQ=255)
                                                                 gATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                         <  1:3633657/51‑1 (MQ=255)
                                                                 gATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                         <  1:2715499/51‑1 (MQ=255)
                                                                 gATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCtt               <  1:1619768/61‑1 (MQ=255)
                                                                 gATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATc     >  1:2637854/1‑71 (MQ=255)
                                                                 gATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATc     >  1:2941385/1‑71 (MQ=255)
                                                                  aacAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                      <  1:338849/51‑1 (MQ=255)
                                                                  aTCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTAc                                  >  1:1366369/1‑41 (MQ=255)
                                                                  aTCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTAc                                  >  1:1320875/1‑41 (MQ=255)
                                                                   tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt         >  1:218816/1‑65 (MQ=255)
                                                                   tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATccc   >  1:913302/1‑71 (MQ=255)
                                                                     aaTCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCg  <  1:322994/70‑1 (MQ=255)
                                                                     |                                                                     
GGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCG  >  minE/313732‑313870

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: