Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,829 G→A 19.1% I348I (ATC→ATT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8290GA19.1% 273.4 / 48.7 120I348I (ATC→ATTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (51/46);  new base A (12/11);  total (63/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.06e-01

Only 100 of 131 total aligned reads displayed.

CGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGG  >  minE/311759‑311898
                                                                      |                                                                     
cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                                       <  1:4147432/71‑1 (MQ=255)
  gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatgat                                                                     <  1:5089969/71‑1 (MQ=255)
  gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatgat                                                                     <  1:1848565/71‑1 (MQ=255)
    aaTGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGACCACCAGcataatg                                                                    >  1:5560352/1‑69 (MQ=2)
       ggTCGGAATCAACGTCGCGAGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCtt                                                                >  1:5645169/1‑71 (MQ=255)
          cGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                                                             >  1:1014996/1‑71 (MQ=255)
             aaTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                                                          <  1:5031656/71‑1 (MQ=255)
              aTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAg                                                         >  1:2390574/1‑71 (MQ=255)
              aTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAg                                                         >  1:1071591/1‑71 (MQ=255)
                 aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                                                       <  1:3908644/70‑1 (MQ=255)
                   cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                                    >  1:2748303/1‑71 (MQ=255)
                   cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                                    >  1:1417996/1‑71 (MQ=255)
                    gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                                   <  1:1484380/71‑1 (MQ=255)
                    gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                                   <  1:3121986/71‑1 (MQ=255)
                    gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                                   <  1:3234591/71‑1 (MQ=255)
                    gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                                   <  1:4426564/71‑1 (MQ=255)
                    gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                                   <  1:950230/71‑1 (MQ=255)
                     tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                                                       <  1:1937636/66‑1 (MQ=255)
                     tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                  <  1:3918965/71‑1 (MQ=255)
                     tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                  <  1:1121907/71‑1 (MQ=255)
                       gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                 >  1:3370189/1‑70 (MQ=255)
                       gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                 >  1:1197818/1‑70 (MQ=255)
                          cgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCAc                                             <  1:4399059/71‑1 (MQ=255)
                           ggAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACt                                            >  1:3698759/1‑71 (MQ=255)
                              aGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                                                             >  1:3504734/1‑51 (MQ=255)
                                ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                                                          <  1:4535045/52‑1 (MQ=255)
                                ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                                                          <  1:5266851/52‑1 (MQ=255)
                                ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                                                          <  1:4895700/52‑1 (MQ=255)
                                ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                       <  1:4650009/71‑1 (MQ=255)
                                ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                       <  1:5126597/71‑1 (MQ=255)
                                ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                       <  1:1266517/71‑1 (MQ=255)
                                ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                       <  1:1448128/71‑1 (MQ=255)
                                ttCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCg                                                              >  1:4140710/1‑48 (MQ=12)
                                ttCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGt                                        >  1:2819597/1‑70 (MQ=18)
                                 tCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTg                                       >  1:3099907/1‑70 (MQ=18)
                                  cTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGc                                                        >  1:2516665/1‑52 (MQ=255)
                                  cTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                     >  1:2686001/1‑71 (MQ=255)
                                     caagAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGaa                                     >  1:4058890/4‑68 (MQ=21)
                                     caagAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGaa                                     >  1:716852/4‑68 (MQ=21)
                                     cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAAt                                    <  1:3558846/69‑1 (MQ=255)
                                     cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAAt                                    <  1:207952/69‑1 (MQ=255)
                                      aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                 <  1:836204/71‑1 (MQ=255)
                                      aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                 <  1:1660673/71‑1 (MQ=255)
                                       agAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATaga                                 >  1:224138/2‑70 (MQ=21)
                                       agAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATaga                                 >  1:2064462/2‑70 (MQ=21)
                                        gAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATagaaa                               <  1:4424764/71‑3 (MQ=18)
                                          aCAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATagaaa                               >  1:3916900/1‑67 (MQ=18)
                                          aCAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCa                                              <  1:2979806/54‑1 (MQ=255)
                                          aCAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtt                             >  1:4570445/1‑71 (MQ=255)
                                          aCAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtt                             >  1:4664282/1‑71 (MQ=255)
                                              aTACATCACCAGGAACACCAGCATAATGGCTTCGaa                                                            <  1:3142984/36‑2 (MQ=12)
                                              aTACATCACCAGGAACACCAGCATAATGGCTTCGaa                                                            <  1:2692604/36‑2 (MQ=12)
                                              aTACATAACGAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                         >  1:1882947/1‑71 (MQ=255)
                                              aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                     <  1:1074778/59‑1 (MQ=255)
                                                aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc                        >  1:4986943/1‑70 (MQ=255)
                                                 cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                      <  1:1791893/71‑1 (MQ=255)
                                                 cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                      <  1:4763334/71‑1 (MQ=255)
                                                 cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                      <  1:2122378/71‑1 (MQ=255)
                                                 cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                      <  1:3941275/71‑1 (MQ=255)
                                                  aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                         >  1:4590069/1‑67 (MQ=255)
                                                  aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTACGTGAATAGAGATTTTCa                         >  1:4825532/1‑67 (MQ=255)
                                                   tcaCCAGGAACACCAGCATAATGGCTTCGAACAGCGtttt                                                   <  1:3807416/38‑1 (MQ=16)
                                                   tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc                    <  1:4887055/71‑1 (MQ=255)
                                                    aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                 >  1:5133547/1‑57 (MQ=255)
                                                    aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                 >  1:1121544/1‑57 (MQ=255)
                                                    aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                 >  1:573596/1‑57 (MQ=255)
                                                    aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                 >  1:3942716/1‑57 (MQ=255)
                                                     aacAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                 >  1:4258127/3‑40 (MQ=255)
                                                     aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                 >  1:4075661/1‑40 (MQ=255)
                                                     aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                 >  1:3716009/1‑40 (MQ=255)
                                                     aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                 >  1:2845680/1‑40 (MQ=255)
                                                       cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                  <  1:163417/37‑1 (MQ=255)
                                                       cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                  <  1:4088301/37‑1 (MQ=255)
                                                        aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt               <  1:5323838/71‑1 (MQ=255)
                                                         ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg              >  1:2177368/1‑71 (MQ=255)
                                                         ggAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtg              <  1:3748263/71‑1 (MQ=255)
                                                         ggAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtg              >  1:4970918/1‑71 (MQ=255)
                                                          gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                           <  1:1019905/57‑1 (MQ=255)
                                                          gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                           <  1:326925/57‑1 (MQ=255)
                                                          gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAt                              >  1:2514016/1‑54 (MQ=255)
                                                          gAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATag                                  <  1:4688178/50‑1 (MQ=255)
                                                          gAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATag                                  <  1:4258475/50‑1 (MQ=255)
                                                          gAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtgg             >  1:4428081/1‑71 (MQ=255)
                                                          gAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtgg             >  1:5289143/1‑71 (MQ=255)
                                                           aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt            <  1:2953439/71‑1 (MQ=255)
                                                           aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt            >  1:3048319/1‑71 (MQ=255)
                                                           aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg              >  1:2557573/1‑69 (MQ=255)
                                                           aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg              >  1:3123258/1‑69 (MQ=255)
                                                            acacGAGGATGATCGCTTCGACCAGCGTTTTAACCa                                              <  1:4256686/36‑1 (MQ=255)
                                                            acacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggtg           <  1:4196140/71‑1 (MQ=255)
                                                              acGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGgtgt          <  1:1719512/70‑1 (MQ=255)
                                                               cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                     <  1:1091756/42‑1 (MQ=255)
                                                               cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                         >  1:1715554/1‑54 (MQ=255)
                                                                gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg        >  1:2362560/1‑70 (MQ=255)
                                                                 aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                         >  1:1313830/1‑52 (MQ=255)
                                                                 aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa      <  1:3480806/71‑1 (MQ=255)
                                                                 aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGAGTGGTGTCGTa      >  1:2278972/1‑71 (MQ=255)
                                                                     tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                   <  1:2902771/54‑1 (MQ=255)
                                                                     tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATggg  >  1:3992616/1‑71 (MQ=255)
                                                                      gatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                       <  1:2354159/49‑1 (MQ=255)
                                                                      |                                                                     
CGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGG  >  minE/311759‑311898

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: