Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,798 G→A 17.5% L359L (CTG→TTG)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7980GA17.5% 214.9 / 29.3 80L359L (CTG→TTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (25/41);  new base A (9/5);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 8.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

ACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAG  >  minE/311728‑311866
                                                                      |                                                                    
aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                                                      >  1:2153029/1‑71 (MQ=255)
 cGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                                                      >  1:1111932/1‑70 (MQ=255)
   gCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                                   <  1:2798570/71‑1 (MQ=255)
      aaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                                   >  1:98439/1‑68 (MQ=255)
      aaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                                   >  1:92396/1‑68 (MQ=255)
        ggTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGAt                                                              <  1:2509573/71‑1 (MQ=255)
            ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                          <  1:2069847/71‑1 (MQ=255)
            ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                          <  1:3224834/71‑1 (MQ=255)
            ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                          <  1:1341076/71‑1 (MQ=255)
             cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                          >  1:1429650/1‑70 (MQ=255)
             cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                                                         <  1:3872070/71‑1 (MQ=255)
              gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCTCGCGGAAGTTCTGCAGGAACAGATACATaa                                                        >  1:119129/1‑71 (MQ=255)
              gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                        >  1:3833828/1‑71 (MQ=255)
              gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                        >  1:1819928/1‑71 (MQ=255)
              gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                        >  1:3168855/1‑71 (MQ=255)
              gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                        >  1:2923309/1‑71 (MQ=255)
               aacAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                       >  1:3808043/3‑71 (MQ=2)
               aacAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                       >  1:571554/3‑71 (MQ=2)
               aacAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                       >  1:1167036/3‑71 (MQ=2)
               aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                                                       <  1:1282483/71‑1 (MQ=255)
               aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                                                       >  1:2781428/1‑71 (MQ=255)
                gCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                                                       >  1:3407060/1‑70 (MQ=255)
                acAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                       <  1:1004605/69‑1 (MQ=2)
                 cAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                                                       >  1:542339/1‑69 (MQ=255)
                 cAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCa                                                     <  1:1035660/71‑1 (MQ=255)
                 cAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                       >  1:861068/1‑69 (MQ=2)
                    tACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACa                                                           <  1:2039763/62‑1 (MQ=255)
                    caccacgGGTACCGCAATGGTGGGGATCAGCGTTTCTCGCATATTCTGCAAGAACAGATACATCACCAgg                                                   >  1:198794/2‑70 (MQ=9)
                    caccacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGa                                                  <  1:1648415/70‑1 (MQ=9)
                    caccacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGa                                                  <  1:1335111/70‑1 (MQ=9)
                      ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                 >  1:1643001/1‑70 (MQ=255)
                      ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                 >  1:2530623/1‑70 (MQ=255)
                      ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAac                                                <  1:2281146/71‑1 (MQ=255)
                       caccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                 >  1:2588772/1‑69 (MQ=255)
                        accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                                              <  1:1917300/71‑1 (MQ=255)
                        accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACAGGAacac                                              >  1:376607/1‑70 (MQ=255)
                         ccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGa                                                  <  1:600954/66‑1 (MQ=255)
                          cGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                                              >  1:3109271/1‑69 (MQ=255)
                               ccgCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGcat                                        <  1:3386862/68‑1 (MQ=9)
                               cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                                   <  1:3751181/59‑1 (MQ=255)
                               cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                       <  1:2025039/71‑1 (MQ=255)
                                ggCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatga                                      <  1:361876/71‑1 (MQ=255)
                                cgCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGcataa                                      >  1:363357/2‑69 (MQ=9)
                                 gCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACAtca                                                        >  1:726613/1‑50 (MQ=2)
                                  cAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATc                                    >  1:2326278/1‑71 (MQ=255)
                                   aaTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATc                                    <  1:1798241/70‑1 (MQ=255)
                                   aaTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCg                                   <  1:1150541/71‑1 (MQ=255)
                                    aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                                                                 <  1:2911672/40‑1 (MQ=255)
                                     tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCt                                 >  1:207882/1‑71 (MQ=255)
                                      ggTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCtt                                >  1:3127576/1‑71 (MQ=2)
                                        tCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACa                                                           <  1:148009/42‑1 (MQ=255)
                                        tCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCg                              <  1:2622455/71‑1 (MQ=255)
                                            aaTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                                             <  1:432043/52‑1 (MQ=255)
                                               cAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                                   <  1:2055783/43‑1 (MQ=255)
                                               cAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                                   <  1:2428334/43‑1 (MQ=255)
                                                aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                 <  1:3018264/44‑1 (MQ=255)
                                                aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                       <  1:3276463/70‑1 (MQ=255)
                                                 aCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatgat                                     >  1:3532092/1‑55 (MQ=255)
                                                 aCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtt                     <  1:884341/71‑1 (MQ=255)
                                                  cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                    <  1:3800745/71‑1 (MQ=255)
                                                   gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAg                                                    >  1:239906/1‑38 (MQ=255)
                                                    tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                  <  1:2852763/71‑1 (MQ=255)
                                                        gcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                       <  1:1585492/62‑1 (MQ=255)
                                                           gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                  <  1:568395/64‑1 (MQ=255)
                                                           gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACtt           <  1:1908208/71‑1 (MQ=255)
                                                           gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACtt           >  1:2785062/1‑71 (MQ=255)
                                                           gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACt            >  1:243692/1‑70 (MQ=255)
                                                            atattCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTc          >  1:688874/4‑71 (MQ=16)
                                                            aaGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc          <  1:2326430/71‑1 (MQ=255)
                                                             aGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCg         <  1:576963/71‑1 (MQ=255)
                                                              gTTCTGCAGGAACAGATACATAACCAGGAACACGAgg                                          <  1:1324104/37‑1 (MQ=255)
                                                              gTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt        <  1:3177601/71‑1 (MQ=255)
                                                               ttCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTg       <  1:2640254/71‑1 (MQ=21)
                                                                 cTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGc                                  >  1:3766991/1‑42 (MQ=255)
                                                                 cTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                >  1:2822574/1‑60 (MQ=255)
                                                                   gCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                  <  1:2295321/56‑1 (MQ=255)
                                                                    cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa   <  1:2749706/70‑1 (MQ=255)
                                                                    cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCCTGAATa   <  1:912443/70‑1 (MQ=255)
                                                                     aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                             >  1:23756/1‑43 (MQ=255)
                                                                     aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATag  <  1:998935/70‑1 (MQ=255)
                                                                     aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATag  <  1:1467086/70‑1 (MQ=255)
                                                                     aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATag  <  1:41066/70‑1 (MQ=255)
                                                                      ggAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa   >  1:3036528/1‑68 (MQ=255)
                                                                      |                                                                    
ACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAG  >  minE/311728‑311866

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