Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,457 A→G 10.6% I472I (ATT→ATC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4570AG10.6% 230.8 / 17.3 85I472I (ATT→ATCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (33/43);  new base G (3/6);  total (36/49)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

GTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGA  >  minE/311389‑311525
                                                                    |                                                                    
gtggCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTaa                                                                     <  1:3944650/70‑1 (MQ=255)
   gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                                   >  1:2363516/1‑69 (MQ=255)
   gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                                   >  1:1001928/1‑69 (MQ=255)
    cacaAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTaa                                                                     <  1:2640482/66‑1 (MQ=255)
          gcgCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGa                                                          >  1:181615/3‑71 (MQ=255)
           agCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                         <  1:1056014/71‑1 (MQ=255)
                gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGa                                                          >  1:342043/1‑65 (MQ=255)
                gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgatag                                                    <  1:115343/71‑1 (MQ=255)
                gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgat                                                      <  1:3740230/69‑1 (MQ=255)
                     aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag                                                >  1:362022/1‑70 (MQ=255)
                     aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag                                                >  1:3684563/1‑70 (MQ=255)
                       gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                                   <  1:2757028/49‑1 (MQ=255)
                         tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATaga                                                                 >  1:425825/1‑49 (MQ=255)
                         tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                             >  1:1124834/1‑69 (MQ=255)
                         tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                             >  1:3731864/1‑69 (MQ=255)
                         tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                           <  1:3324/71‑1 (MQ=255)
                          cAATGCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAgaga                                                               <  1:324753/50‑1 (MQ=255)
                          cAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                                                       >  1:941713/1‑58 (MQ=255)
                           aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAg                                          >  1:1952394/4‑70 (MQ=9)
                           aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGAGATAAATTGCCCCAg                                          >  1:239808/4‑70 (MQ=9)
                           aaCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATagag                                                                <  1:1619759/48‑1 (MQ=255)
                           aaCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg                                          <  1:3192083/70‑1 (MQ=255)
                            aCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACg                                                        <  1:2226893/55‑1 (MQ=255)
                             tgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAgaga                                                               <  1:3478613/46‑1 (MQ=19)
                             cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                       <  1:2575960/71‑1 (MQ=255)
                              gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                           >  1:398186/1‑66 (MQ=255)
                              gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTa                                        >  1:3591795/1‑69 (MQ=255)
                                  ccaGTACCGACAGCGCCATTGCTGAAACAATGGTAATagag                                                                <  1:1045963/41‑1 (MQ=255)
                                   cagaacagaaAGCGCCATTGCCGAAACGATGGTGATAGAGAAc                                                             >  1:76246/1‑43 (MQ=16)
                                   caGTACCGACAGCGCCATTGCTGAAACAATGGTAATATAGAACTGACgat                                                      <  1:2672565/50‑1 (MQ=255)
                                   caGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgat                                                      <  1:2252326/50‑1 (MQ=255)
                                   caGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                                                       >  1:976084/1‑49 (MQ=255)
                                    aGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                         >  1:1847989/1‑46 (MQ=255)
                                    aGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                         >  1:3109767/1‑46 (MQ=255)
                                     gTACCGACAGCGCCATTGCTGAAACAATGGTAATag                                                                  <  1:3117052/36‑1 (MQ=255)
                                     gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                               <  1:798959/71‑1 (MQ=255)
                                      tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                         <  1:742265/44‑1 (MQ=255)
                                      tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                       <  1:1599318/62‑1 (MQ=255)
                                      tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                       <  1:3384108/62‑1 (MQ=255)
                                        ccGACAGCGCCATTGCTGAAACAATGGTAATAGCGAACt                                                            >  1:1596407/1‑39 (MQ=255)
                                        ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                                                       >  1:3956652/1‑44 (MQ=255)
                                        ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCg                                  >  1:2230114/1‑65 (MQ=255)
                                        ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                               >  1:4080458/1‑68 (MQ=255)
                                        ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag                            <  1:2055908/71‑1 (MQ=255)
                                        ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag                            <  1:4008394/71‑1 (MQ=255)
                                         cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgat                                                      <  1:3272867/44‑1 (MQ=255)
                                           aCAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaag                         >  1:1475278/1‑71 (MQ=255)
                                           aCAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                          >  1:942825/1‑70 (MQ=255)
                                               cgcCATTGCTGAAACAATTGTAATAGAGAACTGACGATAGATAGCACCa                                           <  1:434122/49‑1 (MQ=255)
                                                gcCATTGCTGAAACAATGGTAATAGAGAACTGACgata                                                     <  1:2975359/38‑1 (MQ=255)
                                                gcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGaag                            <  1:3311858/63‑1 (MQ=12)
                                                 ccATTGCTGAAACAATGGTAATAGAGAACTGACgata                                                     <  1:2301310/37‑1 (MQ=255)
                                                   aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                             <  1:217482/59‑1 (MQ=255)
                                                    ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                       <  1:2775942/48‑1 (MQ=255)
                                                    ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                 >  1:1240749/1‑70 (MQ=255)
                                                    ttGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAAcc                                   <  1:2349282/52‑1 (MQ=12)
                                                    ttGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGaag                            <  1:2683586/59‑1 (MQ=9)
                                                     tGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                             <  1:4143381/57‑1 (MQ=255)
                                                      gCTGAAACAATGGTAATAGAGAACTGACGATAGATATCACCAg                                          <  1:5109/43‑1 (MQ=255)
                                                      gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg                                          <  1:3609478/43‑1 (MQ=255)
                                                      gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg                                          <  1:3441354/43‑1 (MQ=255)
                                                      gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                                     >  1:2858793/1‑48 (MQ=255)
                                                      gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                                     >  1:3383940/1‑48 (MQ=255)
                                                      gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc              >  1:3593586/1‑71 (MQ=255)
                                                      gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc              >  1:3313949/1‑71 (MQ=255)
                                                      gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc              >  1:884535/1‑71 (MQ=255)
                                                       cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc              >  1:1052156/1‑70 (MQ=255)
                                                       cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc              >  1:1075745/1‑70 (MQ=255)
                                                       cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg             <  1:2634234/71‑1 (MQ=255)
                                                       cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg             <  1:2696760/71‑1 (MQ=255)
                                                       cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg             <  1:3604915/71‑1 (MQ=255)
                                                       cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg             <  1:445949/71‑1 (MQ=255)
                                                         gAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGaaga                           <  1:612709/55‑1 (MQ=12)
                                                          aaaCAATGGTAATAGAGAACTGAGGATAGATAGCAc                                             <  1:1026948/36‑1 (MQ=255)
                                                          aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa           <  1:3294276/70‑1 (MQ=255)
                                                          aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa           >  1:1214635/1‑70 (MQ=255)
                                                          aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa           >  1:2967883/1‑70 (MQ=255)
                                                          aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt          <  1:456851/71‑1 (MQ=255)
                                                          aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt          <  1:726494/71‑1 (MQ=255)
                                                           aaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                               <  1:3808073/49‑1 (MQ=255)
                                                            aCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg             <  1:928449/66‑1 (MQ=255)
                                                             cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAcc       <  1:603014/71‑1 (MQ=255)
                                                             cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAcc       <  1:2620980/71‑1 (MQ=255)
                                                                tGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAccgcc    <  1:1802132/71‑1 (MQ=255)
                                                                  gTAATAGAGAACTGACGATAGATAGCACCAGTAGAAc                                    >  1:760279/1‑37 (MQ=255)
                                                                  gTAATAGAGAACTGACGATAGATAGCACCAGTAGAAc                                    >  1:3733192/1‑37 (MQ=255)
                                                                  gTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGa  <  1:73839/71‑1 (MQ=255)
                                                                    |                                                                    
GTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGA  >  minE/311389‑311525

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