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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE313,8010AG13.7% 247.4 / 27.0 102I96I (ATT→ATCacrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (37/51);  new base G (6/8);  total (43/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 106 total aligned reads displayed.

GGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCC  >  minE/313732‑313869
                                                                     |                                                                    
ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                      <  1:2647068/70‑1 (MQ=255)
ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                      <  1:2673492/70‑1 (MQ=255)
         cTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                                              >  1:4880331/1‑69 (MQ=255)
          tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                                              >  1:2990771/1‑68 (MQ=255)
            ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagaga                                                          <  1:3043379/70‑1 (MQ=255)
            ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                                           <  1:1591286/69‑1 (MQ=255)
            ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                                           <  1:3598000/69‑1 (MQ=255)
             gCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagaga                                                          >  1:2471453/1‑69 (MQ=255)
              cTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                                                     <  1:1300814/57‑1 (MQ=255)
              cTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                                                     <  1:2323984/57‑1 (MQ=255)
               tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                             <  1:3241212/64‑1 (MQ=255)
               tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                             <  1:4010275/64‑1 (MQ=255)
               tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                             <  1:5227791/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)
                ttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                             <  1:121238/63‑1 (MQ=255)
                 tCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTgg                                                                <  1:1020961/59‑2 (MQ=255)
                 tCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggactg                                                     <  1:762063/70‑4 (MQ=255)
                 tCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggactg                                                     <  1:1499681/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)
                        aTCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCtt                                             >  1:499912/1‑71 (MQ=255)
                         tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                            <  1:3853621/71‑1 (MQ=255)
                         tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                            <  1:5018371/71‑1 (MQ=255)
                          cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                            >  1:67833/1‑70 (MQ=255)
                          cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                            >  1:3260999/1‑70 (MQ=255)
                          cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                            >  1:3171989/1‑70 (MQ=255)
                          cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                            >  1:2200629/1‑70 (MQ=255)
                            ccTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTTATagagaga                                                        <  1:2828308/56‑1 (MQ=255)
                            ccTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagagaga                                                        <  1:4006488/56‑1 (MQ=255)
                             cTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagagaga                                                        <  1:304061/55‑1 (MQ=255)
                             cTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                                      >  1:3375539/1‑57 (MQ=255)
                              tttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacaggga                                                        <  1:2690401/54‑1 (MQ=255)
                              tttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacaggga                                                        <  1:367153/54‑1 (MQ=255)
                              tttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacaggga                                                        <  1:5173001/54‑1 (MQ=255)
                              tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                            <  1:4040099/50‑1 (MQ=255)
                              tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                            <  1:3441785/50‑1 (MQ=255)
                               ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGc                                               >  1:5028477/1‑62 (MQ=255)
                                tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTgg                                                                >  1:4448477/1‑43 (MQ=255)
                                tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTgg                                                                >  1:3643475/1‑43 (MQ=255)
                                tCGCACTGTCGTATGTCGCCTGATAGGTCGCAGTATCAATCTGATAGAGAGAGacac                                                   >  1:512899/1‑57 (MQ=255)
                                tCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCt                                                                  >  1:4390208/1‑42 (MQ=255)
                                tCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                                           >  1:817059/1‑49 (MQ=255)
                                tCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGacac                                                   >  1:2912821/1‑57 (MQ=255)
                                 cGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCg                                     <  1:5474679/70‑1 (MQ=255)
                                  gCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCt                                                                  <  1:2331018/40‑1 (MQ=255)
                                      tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                            <  1:2973631/42‑1 (MQ=255)
                                      tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                            <  1:40323/42‑1 (MQ=255)
                                      tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTc                                      <  1:2271121/64‑1 (MQ=255)
                                      tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTc                                      <  1:5569640/64‑1 (MQ=255)
                                        tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagaga                                                          >  1:5534865/1‑42 (MQ=255)
                                        tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCt                                   <  1:1916086/65‑1 (MQ=255)
                                        tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCt                                   <  1:5526433/65‑1 (MQ=255)
                                        tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttc                             <  1:3383050/71‑1 (MQ=255)
                                        tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACct                               <  1:3160888/69‑1 (MQ=255)
                                          gTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttctt                           >  1:5369406/1‑71 (MQ=255)
                                          gTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttct                            <  1:4555267/70‑1 (MQ=255)
                                            aTGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt                          >  1:2100569/1‑70 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGGTCAATCGGATaga                                                            <  1:741158/35‑1 (MQ=38)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                            <  1:4241240/35‑1 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                            <  1:3258420/35‑1 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGATAGACACCTGCTTc                                            <  1:680537/51‑1 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                            <  1:4658673/51‑1 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                        <  1:1209712/71‑1 (MQ=255)
                                             tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                        <  1:1669876/71‑1 (MQ=255)
                                              gTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                        >  1:4420989/1‑70 (MQ=255)
                                               tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                            <  1:1296612/49‑1 (MQ=255)
                                               tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                            <  1:3872741/49‑1 (MQ=255)
                                               tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTg                         >  1:4657609/1‑68 (MQ=255)
                                               tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTg                         >  1:4620502/1‑68 (MQ=255)
                                               tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                      >  1:4796847/1‑71 (MQ=255)
                                               tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                      >  1:4803369/1‑71 (MQ=255)
                                               tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                      >  1:884865/1‑71 (MQ=255)
                                               tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                      >  1:4658335/1‑71 (MQ=255)
                                               tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                      >  1:2910723/1‑71 (MQ=255)
                                                 gCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                     <  1:5370912/70‑1 (MQ=255)
                                                 gCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                     <  1:1600183/70‑1 (MQ=255)
                                                 gCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                     >  1:4009006/1‑70 (MQ=255)
                                                   cTGATAGGTCGCGGGATCGATCTGGTACAGGGACTGGCCTGCTTGCACATCGCTGCCTTCAGTGAAAttgc                  <  1:1309565/71‑3 (MQ=255)
                                                   cTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTAc                  <  1:2847420/71‑1 (MQ=255)
                                                   cTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTAc                  <  1:2489528/71‑1 (MQ=255)
                                                     gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                        <  1:1383000/63‑1 (MQ=255)
                                                     gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGc                >  1:5075756/1‑71 (MQ=255)
                                                     gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGc                >  1:1979692/1‑71 (MQ=255)
                                                             gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt        <  1:261144/71‑1 (MQ=255)
                                                             gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt        <  1:772659/71‑1 (MQ=255)
                                                               aGGATCAATCTGAAAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGTTTCAGGATaa      <  1:2256437/71‑1 (MQ=255)
                                                                 gATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                        <  1:4522382/51‑1 (MQ=255)
                                                                  aTCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttctt                           <  1:4295037/47‑1 (MQ=255)
                                                                  aTCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                        <  1:3180622/50‑1 (MQ=255)
                                                                   tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                     >  1:5443208/1‑52 (MQ=255)
                                                                   tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                     >  1:4254792/1‑52 (MQ=255)
                                                                   tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATcc   >  1:88970/1‑70 (MQ=255)
                                                                    cAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCa            >  1:637433/1‑60 (MQ=255)
                                                                     |                                                                    
GGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCC  >  minE/313732‑313869

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: