Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,724 A→G 9.0% L383L (CTT→CTC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7240AG 9.0% 257.9 / 12.8 93L383L (CTT→CTCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (33/51);  new base G (3/5);  total (36/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCT  >  minE/311658‑311794
                                                                  |                                                                      
caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                                    <  1:2470315/71‑1 (MQ=255)
 agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                                                   <  1:2299073/69‑2 (MQ=12)
 agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                                                   >  1:2319320/3‑70 (MQ=12)
 aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                                   <  1:2742190/71‑1 (MQ=255)
  gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGatg                                                                  >  1:3815727/2‑69 (MQ=12)
  gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGatg                                                                  >  1:576654/2‑69 (MQ=12)
  aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                                                                  >  1:3497243/1‑71 (MQ=255)
    aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACgg                                                                 >  1:3534739/1‑70 (MQ=255)
      gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                                              <  1:2363271/71‑1 (MQ=255)
      gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                                              <  1:1796073/71‑1 (MQ=255)
      gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                                              <  1:3864917/71‑1 (MQ=255)
      gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                                              <  1:2124814/71‑1 (MQ=255)
      gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                               <  1:3006533/70‑1 (MQ=255)
      cccGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCa                                                               <  1:690757/69‑1 (MQ=12)
       ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                               <  1:1565651/69‑1 (MQ=255)
        cGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCaaac                                                            <  1:3652167/71‑2 (MQ=16)
          aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt                                                          >  1:1668109/1‑71 (MQ=255)
          aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt                                                          >  1:820374/1‑71 (MQ=255)
          aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt                                                          <  1:3732859/71‑1 (MQ=255)
           tGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTc                                                         >  1:2077061/1‑71 (MQ=255)
           tGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTc                                                         >  1:1306873/1‑71 (MQ=255)
            ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                                                        >  1:3435692/1‑71 (MQ=255)
             gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                                    >  1:2907973/1‑58 (MQ=255)
             gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                                                                  >  1:762918/1‑60 (MQ=255)
             gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                       <  1:3526958/71‑1 (MQ=255)
             gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                       <  1:2262432/71‑1 (MQ=255)
               aaGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                                                   <  1:663174/56‑2 (MQ=9)
                 gCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAg                                                    >  1:2540071/1‑70 (MQ=255)
                  cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                                                                  >  1:3730817/1‑55 (MQ=255)
                  cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                                                                  >  1:992243/1‑55 (MQ=255)
                   aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCaa                                                 <  1:3568196/71‑1 (MQ=255)
                   aCCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaa                                                 <  1:3435123/71‑1 (MQ=12)
                      aTCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGa                                                     >  1:3509140/1‑64 (MQ=255)
                       tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                                                                <  1:865345/52‑1 (MQ=255)
                       tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                                                                <  1:2946875/52‑1 (MQ=255)
                       tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATac                                              <  1:1866087/70‑1 (MQ=255)
                         ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                                                        <  1:1104027/58‑1 (MQ=255)
                         ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATacca                                            <  1:742993/70‑1 (MQ=255)
                          cGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                            <  1:3615318/53‑1 (MQ=255)
                          cGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                            <  1:832086/53‑1 (MQ=255)
                          cGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                            <  1:2476207/53‑1 (MQ=255)
                           gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGGCGGCAAAGGTCCCGAGCAATACCACCg                                         >  1:3583996/1‑71 (MQ=255)
                           gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                                         <  1:2780201/71‑1 (MQ=255)
                           gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                                         <  1:3451732/71‑1 (MQ=255)
                            aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                        <  1:2634466/71‑1 (MQ=255)
                            aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                        <  1:1878588/71‑1 (MQ=255)
                             aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                        >  1:1436073/1‑70 (MQ=255)
                             aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                        <  1:899232/70‑1 (MQ=255)
                              cATTGTTTGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTa                                      <  1:716698/71‑1 (MQ=255)
                              cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTa                                      >  1:41513/1‑71 (MQ=255)
                              cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTa                                      >  1:313318/1‑71 (MQ=255)
                               aTTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTAc                                     >  1:3030104/1‑71 (MQ=255)
                                 tGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGa                                                     >  1:774512/1‑53 (MQ=255)
                                  gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACgg                                   >  1:537101/1‑70 (MQ=255)
                                   ttAGCGTGTTTATCGAGAAGCCAAAGGCGGCTAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGc                                  <  1:1383837/70‑1 (MQ=255)
                                   ttAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                                                        <  1:1533185/48‑1 (MQ=255)
                                    tAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCaa                                                 >  1:1543152/1‑54 (MQ=255)
                                     aGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATaccac                                           >  1:1478909/1‑59 (MQ=255)
                                       cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATgg                             <  1:2288479/71‑1 (MQ=255)
                                       cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATgg                             <  1:1681432/71‑1 (MQ=255)
                                       cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATgg                             <  1:3819116/71‑1 (MQ=255)
                                        gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGt                            >  1:1594944/1‑71 (MQ=255)
                                         tgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                               <  1:621452/35‑1 (MQ=255)
                                             tATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGa                                                     <  1:2700496/41‑1 (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:573145/68‑1 (MQ=255)
                                                   gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCa                    <  1:811350/68‑1 (MQ=255)
                                                   gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                 <  1:3945945/71‑1 (MQ=255)
                                                   gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                 <  1:2862872/71‑1 (MQ=255)
                                                   gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                 >  1:2081468/1‑71 (MQ=255)
                                                    aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCaa                                                 >  1:1138492/1‑38 (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:2247893/68‑1 (MQ=255)
                                                          aaGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgc             <  1:3907619/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:869617/71‑1 (MQ=255)
                                                                gcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGtt    <  1:2214056/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:1939158/50‑1 (MQ=255)
                                                                  aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGaa                       <  1:2324492/50‑1 (MQ=255)
                                                                  aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt  >  1:173217/1‑71 (MQ=255)
                                                                  aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt  <  1:1157785/71‑1 (MQ=255)
                                                                  aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt  <  1:96267/71‑1 (MQ=255)
                                                                  aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt  >  1:2433711/1‑71 (MQ=255)
                                                                  |                                                                      
CAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCT  >  minE/311658‑311794

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: