Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,811 A→C 19.7% V354V (GTT→GTG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8110AC19.7% 175.9 / 42.9 81V354V (GTT→GTGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (25/40);  new base C (10/6);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 9.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACG  >  minE/311741‑311881
                                                                      |                                                                      
cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                                                                        <  1:3872070/71‑1 (MQ=255)
 gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCTCGCGGAAGTTCTGCAGGAACAGATACATaa                                                                       >  1:119129/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)
 gAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                                       >  1:3168855/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)
  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)
       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)
                     cAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATc                                                   >  1:2326278/1‑71 (MQ=255)
                      aaTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATc                                                   <  1:1798241/70‑1 (MQ=255)
                      aaTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCg                                                  <  1:1150541/71‑1 (MQ=255)
                        tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCt                                                >  1:207882/1‑71 (MQ=255)
                         ggTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCtt                                               >  1:3127576/1‑71 (MQ=2)
                           tCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCg                                             <  1:2622455/71‑1 (MQ=255)
                               aaTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                                                            <  1:432043/52‑1 (MQ=255)
                                  cAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                                                  <  1:2428334/43‑1 (MQ=255)
                                  cAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                                                  <  1:2055783/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)
                                               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:41066/70‑1 (MQ=255)
                                                        aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATag                 <  1:998935/70‑1 (MQ=255)
                                                        aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATag                 <  1:1467086/70‑1 (MQ=255)
                                                         ggAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa                  >  1:3036528/1‑68 (MQ=255)
                                                          gAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGtt                                    <  1:1705182/49‑1 (MQ=16)
                                                          gAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa                  >  1:3152570/1‑67 (MQ=255)
                                                           aaCATATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATagaaag             >  1:936893/1‑68 (MQ=255)
                                                           aaCAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATagaaag             >  1:2180403/1‑68 (MQ=18)
                                                            aCAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                 <  1:665536/50‑1 (MQ=255)
                                                             cAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAttt           <  1:2916420/71‑1 (MQ=255)
                                                              aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt          >  1:2720296/1‑71 (MQ=255)
                                                              aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt          <  1:2519929/71‑1 (MQ=255)
                                                               gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt          <  1:2376201/70‑1 (MQ=255)
                                                               gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTc         <  1:16216/71‑1 (MQ=255)
                                                                aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                         >  1:1030851/1‑54 (MQ=255)
                                                                aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa        >  1:2419244/1‑71 (MQ=255)
                                                                 tACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtt                                    >  1:1949651/1‑42 (MQ=255)
                                                                   cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                               <  1:1544294/45‑1 (MQ=255)
                                                                   cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg      >  1:3199439/1‑70 (MQ=255)
                                                                   cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg      >  1:3619781/1‑70 (MQ=255)
                                                                    aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa    <  1:458249/71‑1 (MQ=255)
                                                                     tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg      <  1:1422026/68‑1 (MQ=255)
                                                                      aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg  >  1:3135286/1‑71 (MQ=255)
                                                                      aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg  >  1:3267093/1‑71 (MQ=255)
                                                                      |                                                                      
CGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACG  >  minE/311741‑311881

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: