Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,760 G→C 8.4% A371A (GCC→GCG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7600GC 8.4% 296.1 / 12.2 95A371A (GCC→GCGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (40/47);  new base C (5/3);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

GTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGA  >  minE/311692‑311830
                                                                    |                                                                      
gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACgg                                                                       >  1:537101/1‑70 (MQ=255)
 ttAGCGTGTTTATCGAGAAGCCAAAGGCGGCTAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGc                                                                      <  1:1383837/70‑1 (MQ=255)
     cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATgg                                                                 <  1:1681432/71‑1 (MQ=255)
     cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATgg                                                                 <  1:2288479/71‑1 (MQ=255)
     cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATgg                                                                 <  1:3819116/71‑1 (MQ=255)
      gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGt                                                                >  1:1594944/1‑71 (MQ=255)
            aTCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGaa                                                           >  1:300351/1‑70 (MQ=255)
             tCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAAt                                                          <  1:2453210/70‑1 (MQ=255)
                agaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCGATGGTCGGAATCAAc                                                      >  1:1445490/1‑71 (MQ=255)
                agaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAAc                                                      <  1:3123425/71‑1 (MQ=255)
                 gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCa                                                        <  1:811350/68‑1 (MQ=255)
                 gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCa                                                        <  1:573145/68‑1 (MQ=255)
                 gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                                                     <  1:2862872/71‑1 (MQ=255)
                 gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                                                     >  1:2081468/1‑71 (MQ=255)
                 gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                                                     <  1:3945945/71‑1 (MQ=255)
                  aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGt                                                    >  1:464883/1‑71 (MQ=255)
                  aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGt                                                    >  1:460268/1‑71 (MQ=255)
                    gCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCgg                                                             <  1:3917453/60‑1 (MQ=255)
                    gCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcg                                                  <  1:2149287/71‑1 (MQ=255)
                      cAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATATCACCGGTACGGCAATGGTCGGAATCAACGTcgc                                                 <  1:1414779/70‑1 (MQ=255)
                      cAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgc                                                 <  1:413555/70‑1 (MQ=255)
                       aaaGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTc                                                               <  1:1522741/55‑1 (MQ=255)
                        aaGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcgcg                                              >  1:2744453/1‑71 (MQ=255)
                        aaGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcgcg                                              <  1:1799095/71‑1 (MQ=255)
                        aaGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgc                                                 <  1:3907619/68‑1 (MQ=255)
                        aaGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgc                                                 <  1:2247893/68‑1 (MQ=255)
                         aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcgcg                                              >  1:1885244/1‑70 (MQ=255)
                          ggcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAAt                                                          >  1:2360096/1‑57 (MQ=255)
                          ggcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGGCGGAATCAACGTCGCGCGGa                                            <  1:2395896/71‑1 (MQ=255)
                             ggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcg                                                  <  1:1670308/62‑1 (MQ=255)
                              gcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGtt                                        <  1:2214056/71‑1 (MQ=255)
                              gcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGtt                                        <  1:869617/71‑1 (MQ=255)
                              gcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGtt                                        <  1:1888094/71‑1 (MQ=255)
                               cAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcg                                                >  1:1410075/1‑62 (MQ=255)
                                aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGTCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt                                      <  1:75692/71‑1 (MQ=255)
                                aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGaa                                                           <  1:2324492/50‑1 (MQ=255)
                                aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGaa                                                           <  1:1939158/50‑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)
                                    aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGt                                                    >  1:1321248/1‑53 (MQ=255)
                                    aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                  >  1:2153029/1‑71 (MQ=255)
                                     cGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                  >  1:1111932/1‑70 (MQ=255)
                                       gCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                               <  1:2798570/71‑1 (MQ=255)
                                        cAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                                                     >  1:1753044/1‑48 (MQ=255)
                                        cAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                                                     <  1:1122996/48‑1 (MQ=255)
                                          aaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTc                                                   >  1:3666597/1‑48 (MQ=255)
                                          aaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                               >  1:98439/1‑68 (MQ=255)
                                          aaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                               >  1:92396/1‑68 (MQ=255)
                                           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGtt                                        <  1:1603528/58‑1 (MQ=255)
                                            ggTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGAt                          <  1:2509573/71‑1 (MQ=255)
                                                ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCa                                   >  1:1460361/1‑58 (MQ=255)
                                                ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCa                                   >  1:3491954/1‑58 (MQ=255)
                                                ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                      <  1:2069847/71‑1 (MQ=255)
                                                ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                      <  1:1341076/71‑1 (MQ=255)
                                                ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                      <  1:3224834/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:3833828/1‑71 (MQ=255)
                                                  gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                    >  1:1819928/1‑71 (MQ=255)
                                                  gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                    >  1:2923309/1‑71 (MQ=255)
                                                   aacAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                   >  1:571554/3‑71 (MQ=2)
                                                   aacAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                   >  1:1167036/3‑71 (MQ=2)
                                                   aacAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                   >  1:3808043/3‑71 (MQ=2)
                                                   aGCAATACCACCGGTACGGCAATGGTCGGAATCaa                                                       >  1:2534103/1‑35 (MQ=255)
                                                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                           >  1:1292292/1‑47 (MQ=255)
                                                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                           >  1:2579188/1‑47 (MQ=255)
                                                   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: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)
                                                             ccGGTACGCCAATGGTCGGAATCAACGTCGCGCgg                                             <  1:1916262/35‑1 (MQ=255)
                                                              cGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGt                                         <  1:3039447/38‑1 (MQ=255)
                                                              cGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac          >  1:3109271/1‑69 (MQ=255)
                                                                   cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg               <  1:3751181/59‑1 (MQ=255)
                                                                   cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg   <  1:2025039/71‑1 (MQ=255)
                                                                    ggCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatga  <  1:361876/71‑1 (MQ=255)
                                                                    |                                                                      
GTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGA  >  minE/311692‑311830

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: