Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,798 G→A 5.8% L359L (CTG→TTG)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7980GA 5.8% 289.1 / 5.0 86L359L (CTG→TTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (47/34);  new base A (2/3);  total (49/37)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

AAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAG  >  minE/311733‑311868
                                                                 |                                                                      
aaaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                                     >  1:3787433/1‑69 (MQ=255)
aaaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                                     >  1:2888802/1‑69 (MQ=255)
aaaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                                                                   <  1:245396/71‑1 (MQ=255)
aaaGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                                                                   <  1:736998/71‑1 (MQ=255)
    gTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATa                                                               >  1:2987938/1‑71 (MQ=255)
      cacGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACa                                                             >  1:1551980/3‑71 (MQ=255)
        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)
                 ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAg                                                                  >  1:3541749/1‑55 (MQ=255)
                 ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                   >  1:3295309/1‑70 (MQ=255)
                 ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                   >  1:75649/1‑70 (MQ=255)
                 ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                   >  1:22376/1‑70 (MQ=255)
                 ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                   >  1:1292228/1‑70 (MQ=255)
                 ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                   >  1:3087128/1‑70 (MQ=255)
                   accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAgg                                                                       <  1:2537481/48‑1 (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)
                     cGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                                                    >  1:3081165/1‑49 (MQ=255)
                       gTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAgg                                            <  1:2146543/71‑1 (MQ=255)
                        tACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAgg                                                                       <  1:476794/43‑1 (MQ=255)
                        tACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGa                                                                      <  1:3674377/44‑1 (MQ=255)
                        tACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGa                                                                      <  1:1730221/44‑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:4087146/1‑51 (MQ=255)
                             cAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                          >  1:935959/1‑51 (MQ=255)
                             cAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                          >  1:1229588/1‑51 (MQ=255)
                               aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                                                        <  1:2470786/35‑1 (MQ=255)
                               aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGAt                                                                <  1:52716/43‑1 (MQ=255)
                               aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGAt                                                                <  1:9457/43‑1 (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)
                                 ggTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATAc                                                              >  1:3711296/1‑43 (MQ=255)
                                 ggTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGc                                    >  1:546188/1‑69 (MQ=255)
                                  gTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                                     >  1:3065777/1‑35 (MQ=255)
                                  gTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATAc                                                              >  1:2182661/1‑42 (MQ=255)
                                  gTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATAc                                                              >  1:2456779/1‑42 (MQ=255)
                                  gTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAaca                                                 >  1:229335/1‑55 (MQ=255)
                                  gTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAaca                                                 >  1:3097006/1‑55 (MQ=255)
                                   tCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTc                                 >  1:3433327/1‑70 (MQ=255)
                                     ggAATCAACGTCTCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGAc                              >  1:3416269/1‑71 (MQ=255)
                                     ggAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                               <  1:3150179/70‑1 (MQ=255)
                                     ggAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                               <  1:352620/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:784997/56‑1 (MQ=255)
                                         tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                         <  1:711709/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:2472535/71‑1 (MQ=255)
                                                 gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                  >  1:2901924/1‑71 (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:4132534/68‑1 (MQ=16)
                                                       atattCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTc            <  1:1257570/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)
                                                               cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa       >  1:3926964/1‑68 (MQ=255)
                                                                aGGAACAGATACATAACCAGGAACACGAGgatgat                                       >  1:2525127/1‑35 (MQ=255)
                                                                 ggAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATagag  <  1:3727808/71‑1 (MQ=255)
                                                                 |                                                                      
AAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAG  >  minE/311733‑311868

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: