Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,481 A→C 11.8% G464G (GGT→GGG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4810AC11.8% 228.3 / 24.7 93G464G (GGT→GGGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (31/51);  new base C (5/6);  total (36/57)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAG  >  minE/311414‑311549
                                                                   |                                                                    
tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                     >  1:3731864/1‑69 (MQ=255)
tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                     >  1:1124834/1‑69 (MQ=255)
tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                                                   <  1:3324/71‑1 (MQ=255)
  aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAg                                                                  >  1:1952394/4‑70 (MQ=9)
  aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGAGATAAATTGCCCCAg                                                                  >  1:239808/4‑70 (MQ=9)
  aaCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg                                                                  <  1:3192083/70‑1 (MQ=255)
    cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                               <  1:2575960/71‑1 (MQ=255)
     gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                                                   >  1:398186/1‑66 (MQ=255)
     gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTa                                                                >  1:3591795/1‑69 (MQ=255)
            gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                                                       <  1:798959/71‑1 (MQ=255)
             tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                               <  1:3384108/62‑1 (MQ=255)
             tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                               <  1:1599318/62‑1 (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)
                  aCAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaag                                                 >  1:1475278/1‑71 (MQ=255)
                  aCAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                                                  >  1:942825/1‑70 (MQ=255)
                      cgcCATTGCTGAAACAATTGTAATAGAGAACTGACGATAGATAGCACCa                                                                   <  1:434122/49‑1 (MQ=255)
                       gcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGaag                                                    <  1:3311858/63‑1 (MQ=12)
                          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:3441354/43‑1 (MQ=255)
                             gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg                                                                  <  1:3609478/43‑1 (MQ=255)
                             gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                                                             >  1:3383940/1‑48 (MQ=255)
                             gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                                                             >  1:2858793/1‑48 (MQ=255)
                             gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                      >  1:884535/1‑71 (MQ=255)
                             gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                      >  1:3313949/1‑71 (MQ=255)
                             gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                      >  1:3593586/1‑71 (MQ=255)
                              cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                      >  1:1052156/1‑70 (MQ=255)
                              cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                      >  1:1075745/1‑70 (MQ=255)
                              cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                     <  1:3604915/71‑1 (MQ=255)
                              cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                     <  1:2634234/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:760279/1‑37 (MQ=255)
                                         gTAATAGAGAACTGACGATAGATAGCACCAGTAGAAc                                                            >  1:3733192/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:69656/57‑1 (MQ=255)
                                               gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                                  <  1:54172/57‑1 (MQ=255)
                                               gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                                  <  1:3472467/57‑1 (MQ=255)
                                               gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                                  <  1:2858723/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)
                                                       gcGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGa            <  1:57355/70‑1 (MQ=255)
                                                       aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                                   >  1:1042407/1‑48 (MQ=255)
                                                       aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc              <  1:135907/69‑1 (MQ=255)
                                                       aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc              <  1:2649993/69‑1 (MQ=255)
                                                       aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc              <  1:3675328/69‑1 (MQ=255)
                                                       aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc              <  1:2122659/69‑1 (MQ=255)
                                                        cGATAAATTGCCCCAGTAGAACCGCCGGAGAATGCCATCGGAATAAATAcc                               <  1:2682735/51‑1 (MQ=255)
                                                          atagatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                       <  1:2630602/57‑1 (MQ=255)
                                                          atagatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                       <  1:566775/57‑1 (MQ=255)
                                                           tagatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                     <  1:510594/42‑1 (MQ=255)
                                                            agatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATAcc        >  1:203348/1‑70 (MQ=255)
                                                             gatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                     <  1:2847626/40‑1 (MQ=255)
                                                              atagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc              <  1:136807/62‑1 (MQ=255)
                                                              atagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc              <  1:440967/62‑1 (MQ=255)
                                                                agCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGAc     >  1:1004187/1‑69 (MQ=255)
                                                                 gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa            >  1:2320443/1‑61 (MQ=255)
                                                                 gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCag  <  1:304715/71‑1 (MQ=255)
                                                                 gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCag  <  1:3208573/71‑1 (MQ=255)
                                                                 gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCag  <  1:1857996/71‑1 (MQ=255)
                                                                  cACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAccgcc                            <  1:1104317/44‑1 (MQ=255)
                                                                   aCCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAcc                   <  1:1521534/52‑1 (MQ=255)
                                                                   |                                                                    
TCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAG  >  minE/311414‑311549

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: