Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,790 G→A 14.3% N361N (AAC→AAT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7900GA14.3% 245.6 / 23.0 84N361N (AAC→AATacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (30/42);  new base A (8/4);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

GCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGT  >  minE/311722‑311860
                                                                    |                                                                      
gcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGtt                                                                      <  1:2214056/71‑1 (MQ=255)
gcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGtt                                                                      <  1:1888094/71‑1 (MQ=255)
gcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGtt                                                                      <  1:869617/71‑1 (MQ=255)
  aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGTCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt                                                                    <  1:75692/71‑1 (MQ=255)
  aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt                                                                    >  1:173217/1‑71 (MQ=255)
  aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt                                                                    >  1:2433711/1‑71 (MQ=255)
  aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt                                                                    <  1:96267/71‑1 (MQ=255)
  aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt                                                                    <  1:1157785/71‑1 (MQ=255)
    ggACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                                   >  1:1101363/1‑70 (MQ=255)
     gACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                                   >  1:3538603/1‑69 (MQ=255)
     gACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCa                                                                 <  1:980473/71‑1 (MQ=255)
      aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                                                >  1:2153029/1‑71 (MQ=255)
       cGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                                                >  1:1111932/1‑70 (MQ=255)
         gCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                             <  1:2798570/71‑1 (MQ=255)
            aaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                             >  1:92396/1‑68 (MQ=255)
            aaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                             >  1:98439/1‑68 (MQ=255)
             aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGtt                                                                      <  1:1603528/58‑1 (MQ=255)
              ggTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGAt                                                        <  1:2509573/71‑1 (MQ=255)
                  ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCa                                                                 >  1:3491954/1‑58 (MQ=255)
                  ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCa                                                                 >  1:1460361/1‑58 (MQ=255)
                  ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                    <  1:3224834/71‑1 (MQ=255)
                  ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                    <  1:1341076/71‑1 (MQ=255)
                  ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                    <  1:2069847/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:3168855/1‑71 (MQ=255)
                    gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                  >  1:2923309/1‑71 (MQ=255)
                    gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                  >  1:3833828/1‑71 (MQ=255)
                    gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                  >  1:1819928/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:2781428/1‑71 (MQ=255)
                     aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                                                 <  1:1282483/71‑1 (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:1335111/70‑1 (MQ=9)
                          caccacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGa                                            <  1:1648415/70‑1 (MQ=9)
                            ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                           >  1:2530623/1‑70 (MQ=255)
                            ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                           >  1:1643001/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)
                                cGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGt                                                                       <  1:3039447/38‑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:2785062/1‑71 (MQ=255)
                                                                 gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACtt     <  1:1908208/71‑1 (MQ=255)
                                                                 gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACt      >  1:243692/1‑70 (MQ=255)
                                                                  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)
                                                                    |                                                                      
GCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGT  >  minE/311722‑311860

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: