Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,469 A→G 10.0% R468R (CGT→CGC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4690AG10.0% 251.3 / 16.5 90R468R (CGT→CGCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (34/47);  new base G (5/4);  total (39/51)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGC  >  minE/311399‑311537
                                                                      |                                                                    
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:3684563/1‑70 (MQ=255)
           aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag                                                            >  1:362022/1‑70 (MQ=255)
               tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                         >  1:3731864/1‑69 (MQ=255)
               tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                         >  1:1124834/1‑69 (MQ=255)
               tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                                       <  1:3324/71‑1 (MQ=255)
                cAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                                                                   >  1:941713/1‑58 (MQ=255)
                 aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAg                                                      >  1:1952394/4‑70 (MQ=9)
                 aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGAGATAAATTGCCCCAg                                                      >  1:239808/4‑70 (MQ=9)
                 aaCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg                                                      <  1:3192083/70‑1 (MQ=255)
                  aCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACg                                                                    <  1:2226893/55‑1 (MQ=255)
                   cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                   <  1:2575960/71‑1 (MQ=255)
                    gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                                       >  1:398186/1‑66 (MQ=255)
                    gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTa                                                    >  1:3591795/1‑69 (MQ=255)
                         caGTACCGACAGCGCCATTGCTGAAACAATGGTAATATAGAACTGACgat                                                                  <  1:2672565/50‑1 (MQ=255)
                         caGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgat                                                                  <  1:2252326/50‑1 (MQ=255)
                         caGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                                                                   >  1:976084/1‑49 (MQ=255)
                          aGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                                     >  1:3109767/1‑46 (MQ=255)
                          aGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                                     >  1:1847989/1‑46 (MQ=255)
                           gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                                           <  1:798959/71‑1 (MQ=255)
                            tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                                     <  1:742265/44‑1 (MQ=255)
                            tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                   <  1:3384108/62‑1 (MQ=255)
                            tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                   <  1:1599318/62‑1 (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:4008394/71‑1 (MQ=255)
                              ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag                                        <  1:2055908/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:884535/1‑71 (MQ=255)
                                            gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                          >  1:3313949/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:3604915/71‑1 (MQ=255)
                                             cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                         <  1:2696760/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:3733192/1‑37 (MQ=255)
                                                        gTAATAGAGAACTGACGATAGATAGCACCAGTAGAAc                                                >  1:760279/1‑37 (MQ=255)
                                                        gTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGa              <  1:73839/71‑1 (MQ=255)
                                                           aTAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                                         <  1:1344794/41‑1 (MQ=255)
                                                             agagAACTGACGATAGATAGCACCAGTAGGACCGCCAAAGAAGGCCAt                                >  1:2305752/1‑48 (MQ=255)
                                                              gagaACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCAc        >  1:2879698/1‑71 (MQ=9)
                                                              gagaACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                                      >  1:2737850/1‑41 (MQ=255)
                                                              gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                           >  1:2858637/1‑52 (MQ=255)
                                                              gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                      <  1:3472467/57‑1 (MQ=255)
                                                              gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                      <  1:2858723/57‑1 (MQ=255)
                                                              gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                      <  1:54172/57‑1 (MQ=255)
                                                              gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                      <  1:69656/57‑1 (MQ=255)
                                                               agaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                         <  1:497705/53‑1 (MQ=255)
                                                               agaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAcc       <  1:2908535/71‑1 (MQ=255)
                                                                gaACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGa              >  1:3653763/1‑63 (MQ=9)
                                                                 aaCTGACGATAGATAGCACCAGTAGAACCGCCaaa                                         >  1:1190044/1‑35 (MQ=255)
                                                                 aaCTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAccgcc                <  1:1774125/60‑1 (MQ=255)
                                                                   cTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATcg   >  1:3074943/1‑71 (MQ=255)
                                                                   cTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                       >  1:2611265/1‑51 (MQ=255)
                                                                     gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc  >  1:2217854/1‑70 (MQ=255)
                                                                     gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc  >  1:2164353/1‑70 (MQ=255)
                                                                      aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                       >  1:1042407/1‑48 (MQ=255)
                                                                      aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc  <  1:3675328/69‑1 (MQ=255)
                                                                      aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc  <  1:2649993/69‑1 (MQ=255)
                                                                      aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc  <  1:2122659/69‑1 (MQ=255)
                                                                      aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc  <  1:135907/69‑1 (MQ=255)
                                                                      |                                                                    
GAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGC  >  minE/311399‑311537

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: