Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,811 A→C 7.8% 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,8110AC 7.8% 242.7 / 14.3 90V354V (GTT→GTGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (42/41);  new base C (2/5);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACG  >  minE/311741‑311881
                                                                      |                                                                      
cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                                                                        <  1:918311/71‑1 (MQ=255)
  aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                                                                      >  1:3695114/1‑71 (MQ=255)
   gCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAcc                                                                     <  1:171374/71‑1 (MQ=255)
   acAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                                      <  1:613047/69‑1 (MQ=2)
      aTACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                                                  <  1:2620909/71‑1 (MQ=255)
       tACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                                       >  1:1574379/1‑65 (MQ=255)
       tACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                                       >  1:3207598/1‑65 (MQ=255)
        accacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGaa                                                                >  1:3657808/1‑71 (MQ=2)
         ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTGTGCAGGAACAGATACATAACCa                                                                    >  1:1858654/1‑66 (MQ=255)
         ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                                >  1:1292228/1‑70 (MQ=255)
         ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                                >  1:3087128/1‑70 (MQ=255)
         ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                                >  1:22376/1‑70 (MQ=255)
         ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                                >  1:3295309/1‑70 (MQ=255)
         ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                                >  1:75649/1‑70 (MQ=255)
           accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                                                             <  1:2716727/71‑1 (MQ=255)
           accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                                                             <  1:2975405/71‑1 (MQ=255)
            ccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAaca                                                              >  1:2776294/1‑69 (MQ=255)
            ccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                                                            <  1:3403257/71‑1 (MQ=255)
               gTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAgg                                                         <  1:2146543/71‑1 (MQ=255)
                 aCGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAg                                                                   >  1:1099032/1‑59 (MQ=255)
                  cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCa                                                                    >  1:130830/1‑57 (MQ=255)
                  cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                      <  1:3088395/71‑1 (MQ=255)
                     cAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                                       >  1:935959/1‑51 (MQ=255)
                     cAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                                       >  1:4087146/1‑51 (MQ=255)
                     cAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                                       >  1:1229588/1‑51 (MQ=255)
                        tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                                                  <  1:998136/53‑1 (MQ=255)
                        tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAgg                                                         <  1:2365597/62‑1 (MQ=255)
                        tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                      <  1:2181292/65‑1 (MQ=255)
                        tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                      <  1:1804984/65‑1 (MQ=255)
                        tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCt                                                <  1:1241613/71‑1 (MQ=255)
                         ggTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGc                                                 >  1:546188/1‑69 (MQ=255)
                          gTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAaca                                                              >  1:3097006/1‑55 (MQ=255)
                          gTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAaca                                                              >  1:229335/1‑55 (MQ=255)
                           tCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTc                                              >  1:3433327/1‑70 (MQ=255)
                             ggAATCAACGTCTCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGAc                                           >  1:3416269/1‑71 (MQ=255)
                             ggAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                                            <  1:352620/70‑1 (MQ=255)
                             ggAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                                            <  1:3150179/70‑1 (MQ=255)
                             ggAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                                            <  1:2647656/70‑1 (MQ=255)
                               aaTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                                <  1:3400413/48‑1 (MQ=255)
                               aaTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                                          >  1:2835448/1‑54 (MQ=255)
                                aTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAgg                                                         <  1:2680220/54‑1 (MQ=255)
                                aTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                      <  1:4032425/57‑1 (MQ=255)
                                 tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                                                                        <  1:289060/38‑1 (MQ=255)
                                 tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAg                                                                   >  1:4103127/1‑43 (MQ=255)
                                 tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                      <  1:711709/56‑1 (MQ=255)
                                 tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                      <  1:784997/56‑1 (MQ=255)
                                   aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                                          >  1:524253/1‑50 (MQ=255)
                                   aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatgat                                                    >  1:3307076/1‑56 (MQ=255)
                                   aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATc                                                   >  1:2900565/1‑57 (MQ=255)
                                       tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCa                                                                    >  1:3309955/1‑36 (MQ=255)
                                       tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                      <  1:3841210/50‑1 (MQ=255)
                                        cgcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAgg                                                         >  1:2766538/1‑46 (MQ=255)
                                        cgcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                >  1:706391/1‑71 (MQ=255)
                                        cgcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                >  1:3432227/1‑71 (MQ=255)
                                         gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGga                                                        >  1:622806/1‑46 (MQ=255)
                                         gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                >  1:3644790/1‑70 (MQ=255)
                                         gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                               >  1:3917369/1‑71 (MQ=255)
                                         gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                               >  1:2901924/1‑71 (MQ=255)
                                         gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                               <  1:2472535/71‑1 (MQ=255)
                                            cgGAAGTTCTGCAGGAACAGATACATAACCCGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCAc                            <  1:3786173/71‑1 (MQ=255)
                                             ggAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACt                           >  1:2842000/1‑71 (MQ=255)
                                              gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCt                                                <  1:2356451/49‑1 (MQ=255)
                                               atattCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTc                         <  1:1257570/68‑1 (MQ=16)
                                               atattCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTc                         <  1:4132534/68‑1 (MQ=16)
                                                 attCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGt                       >  1:3806065/2‑71 (MQ=18)
                                                  ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                      <  1:3151805/71‑1 (MQ=255)
                                                  ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                      <  1:3051962/71‑1 (MQ=255)
                                                    cTGCAGGAACAGAAACATAACCAGGAACACGAGGAAGATCGCTTCGACCAGCGTTTTAACCACt                           >  1:2082536/1‑64 (MQ=255)
                                                       caagAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGaa                    <  1:2376885/65‑1 (MQ=21)
                                                       cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                    >  1:3926964/1‑68 (MQ=255)
                                                        aagAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGaa                    <  1:3641581/65‑1 (MQ=21)
                                                        aGGAACAGATACATAACCAGGAACACGAGgatgat                                                    >  1:2525127/1‑35 (MQ=255)
                                                         ggAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATagag               <  1:3727808/71‑1 (MQ=255)
                                                            aCAGATACATAACCAGGAACACGAGGATGATCGCt                                                <  1:709071/35‑1 (MQ=255)
                                                            aCAGATACATAACCAGGAACACGAGGATGATCGCt                                                <  1:2457115/35‑1 (MQ=255)
                                                             cAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTAACCACTTCGTGAATAGAGAtttt          <  1:282269/71‑1 (MQ=255)
                                                              aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt          <  1:1906737/71‑1 (MQ=255)
                                                              aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAttt           >  1:3224490/1‑70 (MQ=255)
                                                               gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                   >  1:3809363/1‑45 (MQ=255)
                                                               gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTc         <  1:1770351/71‑1 (MQ=255)
                                                                aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa                  >  1:3316468/1‑61 (MQ=255)
                                                                aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt          <  1:360248/69‑1 (MQ=255)
                                                                 tACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc       <  1:3823302/71‑1 (MQ=255)
                                                                  aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                         <  1:1019695/52‑1 (MQ=255)
                                                                  aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                         <  1:2960730/52‑1 (MQ=255)
                                                                  aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                         <  1:2870579/52‑1 (MQ=255)
                                                                  aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt                       <  1:3108121/54‑1 (MQ=255)
                                                                  aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt                       <  1:683765/54‑1 (MQ=255)
                                                                  aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                      <  1:2335233/55‑1 (MQ=255)
                                                                  aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt          >  1:2972066/1‑67 (MQ=255)
                                                                  aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc       >  1:4121586/1‑70 (MQ=255)
                                                                    aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                      <  1:2657994/53‑1 (MQ=255)
                                                                    aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa    >  1:66945/1‑71 (MQ=255)
                                                                    aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa    >  1:442460/1‑71 (MQ=255)
                                                                    aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa    >  1:906935/1‑71 (MQ=255)
                                                                     tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc   >  1:2195688/1‑71 (MQ=255)
                                                                     tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc   >  1:2112705/1‑71 (MQ=255)
                                                                      aaCCAGGAACACGTGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg  >  1:1061898/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: