Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,823 G→C 18.8% L350L (CTC→CTG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8230GC18.8% 329.9 / 59.9 138L350L (CTC→CTGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (52/60);  new base C (13/13);  total (65/73)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 141 total aligned reads displayed.

GGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCG  >  minE/311755‑311892
                                                                    |                                                                     
ggTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                                                     >  1:1448485/1‑71 (MQ=255)
  tACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGga                                                                   >  1:2878751/1‑71 (MQ=255)
    ccgCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGcat                                                                  <  1:3775086/68‑1 (MQ=9)
    cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                                 <  1:109107/71‑1 (MQ=255)
      gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatgat                                                               <  1:1848565/71‑1 (MQ=255)
        aaTGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGACCACCAGcataatg                                                              >  1:5560352/1‑69 (MQ=2)
         aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                                                                       <  1:2180989/60‑1 (MQ=255)
          tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                                                     <  1:3730336/61‑1 (MQ=255)
          tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                                                     <  1:3283076/61‑1 (MQ=255)
          tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                                                     <  1:3154050/61‑1 (MQ=255)
          tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                                                     <  1:5267957/61‑1 (MQ=255)
           ggTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                                                     >  1:4168476/1‑60 (MQ=255)
           ggTCGGAATCAACGTCGCGAGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCtt                                                          >  1:5645169/1‑71 (MQ=255)
            gTCGGATTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                                                     <  1:4376172/59‑1 (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)
                   tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                                                                       <  1:1713518/50‑1 (MQ=255)
                   tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                                 <  1:1918509/56‑1 (MQ=255)
                   tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                                                 <  1:2094437/56‑1 (MQ=255)
                     aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                                                 <  1:3908644/70‑1 (MQ=255)
                       cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                              >  1:2754883/1‑71 (MQ=255)
                       cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                              >  1:5371946/1‑71 (MQ=255)
                        gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                             <  1:3121986/71‑1 (MQ=255)
                        gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                             <  1:1484380/71‑1 (MQ=255)
                        gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                             <  1:4208913/71‑1 (MQ=255)
                        gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                             <  1:3234591/71‑1 (MQ=255)
                        gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                             <  1:950230/71‑1 (MQ=255)
                        gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                             <  1:4426564/71‑1 (MQ=255)
                         tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                                                 <  1:1937636/66‑1 (MQ=255)
                         tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                            <  1:1121907/71‑1 (MQ=255)
                         tCGCGCGGAAGTTCTGAAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                                                 <  1:2248828/66‑1 (MQ=255)
                           gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                           >  1:3370189/1‑70 (MQ=255)
                           gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                           >  1:1197818/1‑70 (MQ=255)
                           gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                          <  1:2932654/71‑1 (MQ=255)
                              cgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCAc                                       <  1:4399059/71‑1 (MQ=255)
                                gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAgg                                                                    <  1:1332773/40‑1 (MQ=255)
                                  aGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                                                       >  1:3504734/1‑51 (MQ=255)
                                   gTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt                                  >  1:3065788/1‑71 (MQ=255)
                                    ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                                                    <  1:4895700/52‑1 (MQ=255)
                                    ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                                                    <  1:4535045/52‑1 (MQ=255)
                                    ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                 <  1:4650009/71‑1 (MQ=255)
                                    ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                 <  1:5126597/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)
                                     tCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                 >  1:4024720/1‑70 (MQ=255)
                                      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)
                                           agAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATaga                           >  1:2064462/2‑70 (MQ=21)
                                              aCAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATagaaa                         >  1:3916900/1‑67 (MQ=18)
                                              aCAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtt                       >  1:4664282/1‑71 (MQ=255)
                                              aCAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtt                       >  1:4570445/1‑71 (MQ=255)
                                                aGATACATAACCAGGAACACGAGGATGATCGCTTCGAcc                                                     >  1:3670653/1‑39 (MQ=255)
                                                  aTACATCACCAGGAACACCAGCATAATGGCTTCGaa                                                      <  1:376063/36‑2 (MQ=12)
                                                  aTACATCACCAGGAACACCAGCATAATGGCTTCGaa                                                      <  1:2625063/36‑2 (MQ=12)
                                                  aTACATCACCAGGAACACCAGCATAATGGCTTCGaa                                                      <  1:2692604/36‑2 (MQ=12)
                                                  aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                               <  1:1598976/59‑1 (MQ=255)
                                                    aCATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGttttt                                            <  1:5185825/44‑2 (MQ=16)
                                                     cATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGttttt                                            <  1:412807/43‑2 (MQ=16)
                                                     cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                <  1:1791893/71‑1 (MQ=255)
                                                     cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                <  1:4763334/71‑1 (MQ=255)
                                                     cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                <  1:3941275/71‑1 (MQ=255)
                                                      aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                               >  1:3922039/1‑55 (MQ=255)
                                                      aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                   >  1:4590069/1‑67 (MQ=255)
                                                       tcaCCAGGAACACCAGCATAATGGCTTCGAACAGCGtttt                                             <  1:3807416/38‑1 (MQ=16)
                                                       tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc              <  1:4887055/71‑1 (MQ=255)
                                                        aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                             >  1:620188/1‑39 (MQ=255)
                                                        aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                           >  1:5133547/1‑57 (MQ=255)
                                                        aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                           >  1:3212469/1‑57 (MQ=255)
                                                         aacAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                           >  1:4258127/3‑40 (MQ=255)
                                                         aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                           >  1:3716009/1‑40 (MQ=255)
                                                         aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                           >  1:4075661/1‑40 (MQ=255)
                                                           cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                            <  1:163417/37‑1 (MQ=255)
                                                            aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt         <  1:5056854/71‑1 (MQ=255)
                                                             ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg        >  1:1388699/1‑71 (MQ=255)
                                                             ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg        >  1:2177368/1‑71 (MQ=255)
                                                             ggAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtg        >  1:4970918/1‑71 (MQ=255)
                                                              gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                     <  1:1019905/57‑1 (MQ=255)
                                                              gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAt                        >  1:2514016/1‑54 (MQ=255)
                                                              gAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtgg       >  1:4428081/1‑71 (MQ=255)
                                                              gAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtgg       >  1:5289143/1‑71 (MQ=255)
                                                               aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt      >  1:3048319/1‑71 (MQ=255)
                                                               aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt      <  1:2953439/71‑1 (MQ=255)
                                                               aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg        >  1:4267877/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)
                                                                   cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                     >  1:1144417/1‑52 (MQ=255)
                                                                   cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTc                    >  1:4707099/1‑53 (MQ=255)
                                                                   cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                   >  1:1715554/1‑54 (MQ=255)
                                                                    gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg  >  1:2362560/1‑70 (MQ=255)
                                                                    |                                                                     
GGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCG  >  minE/311755‑311892

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: