Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,691 T→C 22.5% T394T (ACA→ACG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,6910TC22.5% 258.9 / 61.6 111T394T (ACA→ACGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (36/50);  new base C (15/10);  total (51/60)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 113 total aligned reads displayed.

CTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTAC  >  minE/311622‑311759
                                                                     |                                                                    
cTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAtt                                                                      <  1:1341490/70‑1 (MQ=255)
   aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGtt                                                                   <  1:2023695/70‑1 (MQ=255)
   aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTa                                                                  <  1:692119/71‑1 (MQ=255)
     cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                                                <  1:848997/71‑1 (MQ=255)
     cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                                                >  1:1053744/1‑71 (MQ=255)
      gTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCg                                                               >  1:992196/1‑71 (MQ=255)
      gTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCg                                                               <  1:2724515/71‑1 (MQ=255)
         ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                                                            >  1:2238548/1‑71 (MQ=255)
         ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                                                            >  1:282127/1‑71 (MQ=255)
          tCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                                           >  1:3860039/1‑71 (MQ=255)
            tACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                                         >  1:3897047/1‑71 (MQ=255)
            tACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                                         >  1:1699142/1‑71 (MQ=255)
            caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                                         >  1:3350701/2‑69 (MQ=9)
            caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                                         >  1:482721/2‑69 (MQ=9)
            caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                                         >  1:93297/2‑69 (MQ=9)
            caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttg                                                          <  1:193640/69‑2 (MQ=9)
            caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttg                                                          <  1:1540903/69‑2 (MQ=9)
             accaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAt                                                        <  1:1366210/71‑1 (MQ=9)
             accaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAt                                                        >  1:2544236/1‑71 (MQ=9)
             aCCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                                        >  1:3452488/1‑71 (MQ=255)
              ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                                        <  1:1304492/70‑1 (MQ=255)
              ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATc                                                       <  1:1727597/71‑1 (MQ=255)
               cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                                         <  1:1515133/68‑1 (MQ=255)
               cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                                                      <  1:564412/71‑1 (MQ=255)
               cacaACGATCGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                                                      <  1:1458078/71‑1 (MQ=255)
                 caACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag                                                    <  1:2228932/71‑1 (MQ=255)
                  aaCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga                                                   <  1:2189356/71‑1 (MQ=255)
                       tGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                               <  1:2971964/70‑1 (MQ=255)
                          cGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAg                                           <  1:3797414/71‑1 (MQ=255)
                           atCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaaaag                                          >  1:3276913/2‑69 (MQ=9)
                           atCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaaaag                                          >  1:2778838/2‑69 (MQ=9)
                            tCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGtt                                                           <  1:1540298/53‑1 (MQ=16)
                            tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                                                      <  1:1899985/58‑1 (MQ=255)
                             cATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAg                                                 >  1:227123/1‑62 (MQ=255)
                              aTCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                                            >  1:563275/1‑66 (MQ=255)
                               tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                      <  1:1824706/71‑1 (MQ=255)
                               tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                      <  1:2575103/71‑1 (MQ=255)
                                   ccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                                           >  1:3401248/1‑46 (MQ=255)
                                     agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagcggc                                      <  1:1648596/63‑1 (MQ=12)
                                     agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                <  1:2299073/69‑2 (MQ=12)
                                     agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                >  1:2319320/3‑70 (MQ=12)
                                     aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                                                >  1:1262751/1‑39 (MQ=255)
                                     aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                                                >  1:426267/1‑39 (MQ=255)
                                     aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                                        <  1:3787862/47‑1 (MQ=255)
                                     aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGaa                                                  >  1:2471585/1‑53 (MQ=255)
                                      gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGatg                               >  1:576654/2‑69 (MQ=12)
                                      gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGatg                               >  1:3815727/2‑69 (MQ=12)
                                      aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                                           >  1:2683269/1‑43 (MQ=255)
                                      aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                                           >  1:1569027/1‑43 (MQ=255)
                                      aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                               >  1:3497243/1‑71 (MQ=255)
                                        aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACgg                              >  1:3534739/1‑70 (MQ=255)
                                         ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGc                                                <  1:3191325/51‑1 (MQ=255)
                                          gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                           <  1:1796073/71‑1 (MQ=255)
                                          gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                           <  1:2124814/71‑1 (MQ=255)
                                          gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                           <  1:2363271/71‑1 (MQ=255)
                                          cccGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCa                            <  1:690757/69‑1 (MQ=12)
                                           ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga                                                   <  1:899035/46‑1 (MQ=255)
                                           ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                            <  1:1565651/69‑1 (MQ=255)
                                            cGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCaaac                         <  1:3652167/71‑2 (MQ=16)
                                             gATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                             >  1:892493/1‑50 (MQ=9)
                                             gATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                             >  1:1247340/1‑50 (MQ=9)
                                             gATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                             >  1:2075371/1‑50 (MQ=9)
                                              aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt                       >  1:1668109/1‑71 (MQ=255)
                                              aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt                       <  1:3732859/71‑1 (MQ=255)
                                              aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt                       >  1:820374/1‑71 (MQ=255)
                                              aTGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAt                                                        >  1:3094098/1‑38 (MQ=17)
                                               tGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTc                      >  1:1306873/1‑71 (MQ=255)
                                               tGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTc                      >  1:2077061/1‑71 (MQ=255)
                                               tGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagcgg                                       >  1:3950389/1‑54 (MQ=9)
                                                ggCTAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                          >  1:192406/1‑50 (MQ=255)
                                                 gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                 >  1:2907973/1‑58 (MQ=255)
                                                 gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                               >  1:762918/1‑60 (MQ=255)
                                                 gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                    <  1:3526958/71‑1 (MQ=255)
                                                 gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                    <  1:2262432/71‑1 (MQ=255)
                                                    aGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                              <  1:772822/42‑1 (MQ=255)
                                                     gCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAg                 >  1:2540071/1‑70 (MQ=255)
                                                      cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                               >  1:3730817/1‑55 (MQ=255)
                                                       aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCaa              <  1:3568196/71‑1 (MQ=255)
                                                       aCCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaa              <  1:3435123/71‑1 (MQ=12)
                                                        ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                      <  1:1463810/46‑1 (MQ=255)
                                                          aTCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGa                  >  1:3509140/1‑64 (MQ=255)
                                                           tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggc                                         <  1:170205/40‑1 (MQ=255)
                                                           tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                             <  1:865345/52‑1 (MQ=255)
                                                           tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                             <  1:2946875/52‑1 (MQ=255)
                                                           tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATac           <  1:1866087/70‑1 (MQ=255)
                                                             ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                     <  1:1104027/58‑1 (MQ=255)
                                                             ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATacca         <  1:742993/70‑1 (MQ=255)
                                                              cGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                         <  1:2476207/53‑1 (MQ=255)
                                                               gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcg                                        <  1:2982814/37‑1 (MQ=255)
                                                               gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGGCGGCAAAGGTCCCGAGCAATACCACCg      >  1:3583996/1‑71 (MQ=255)
                                                               gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg      <  1:2780201/71‑1 (MQ=255)
                                                               gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg      <  1:3451732/71‑1 (MQ=255)
                                                                aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg     <  1:2634466/71‑1 (MQ=255)
                                                                 aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg     <  1:899232/70‑1 (MQ=255)
                                                                 aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg     >  1:1436073/1‑70 (MQ=255)
                                                                  cATTGTTTGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTa   <  1:716698/71‑1 (MQ=255)
                                                                  cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTa   >  1:41513/1‑71 (MQ=255)
                                                                  cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTa   >  1:313318/1‑71 (MQ=255)
                                                                   aTTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTAc  >  1:3030104/1‑71 (MQ=255)
                                                                     tGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGa                  >  1:774512/1‑53 (MQ=255)
                                                                     |                                                                    
CTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTAC  >  minE/311622‑311759

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: