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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4750GA20.8% 275.7 / 52.1 115I466I (ATC→ATTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (38/53);  new base A (9/15);  total (47/68)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 131 total aligned reads displayed.

TCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGA  >  minE/311408‑311545
                                                                   |                                                                      
tCAGGATCAACGCCACCAGTACAGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgatagata                                                                     <  1:4491167/71‑1 (MQ=255)
  agaaTCAATGCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGc                                                                   <  1:4657457/68‑1 (MQ=9)
    gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCa                                                                  >  1:5493319/1‑70 (MQ=255)
    gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                 >  1:2687123/1‑71 (MQ=255)
    gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                 >  1:5538500/1‑71 (MQ=255)
    gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                 >  1:4386025/1‑71 (MQ=255)
    gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                 >  1:5080207/1‑71 (MQ=255)
      tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgataga                                                                       >  1:2619024/1‑63 (MQ=255)
      tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgataga                                                                       >  1:1338692/1‑63 (MQ=255)
      tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                 >  1:454634/1‑69 (MQ=255)
      tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                 >  1:3183534/1‑69 (MQ=255)
        aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGt                                                             <  1:2056951/68‑1 (MQ=9)
          cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                           <  1:5619597/71‑1 (MQ=255)
           gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGa                                                          <  1:228216/71‑1 (MQ=255)
            ccaccaGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                                                         >  1:1947621/1‑71 (MQ=255)
            ccaccaGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCgataaat                                                                      <  1:5341182/58‑1 (MQ=9)
            ccaccaGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCgataaat                                                                      <  1:4087024/58‑1 (MQ=9)
             caccaGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                                                         <  1:1981799/70‑1 (MQ=255)
             caccaGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                                                         <  1:3828326/70‑1 (MQ=255)
                cagaacagaaAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCg                                                      <  1:3048365/70‑1 (MQ=12)
                 aGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGc                                                     >  1:1819492/1‑70 (MQ=255)
                  gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                                                   <  1:93049/71‑1 (MQ=255)
                     ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                           <  1:1057858/60‑1 (MQ=255)
                     ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                                                 <  1:3573257/70‑1 (MQ=255)
                      cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgatag                                                                        <  1:3576800/46‑1 (MQ=255)
                      cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgatag                                                                        <  1:2221125/46‑1 (MQ=255)
                      cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                               <  1:1989027/71‑1 (MQ=255)
                      cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                               <  1:1477252/71‑1 (MQ=255)
                      cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                               <  1:669992/71‑1 (MQ=255)
                       gaaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCTATAAATTGCCCCAGTAGAACCGCCgaagaa                                              <  1:3210604/68‑1 (MQ=12)
                       gaaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCgaagaa                                              <  1:2556783/68‑1 (MQ=12)
                        aaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCgaagaat                                             >  1:4834798/3‑70 (MQ=12)
                        aCAGCGCCATTGCTGAAACAATGGTAATGGAGAACTGACGATAGATAGCAcc                                                                <  1:377635/52‑1 (MQ=255)
                        aCAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaag                                             >  1:3356927/1‑71 (MQ=255)
                         cAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAgg                                            >  1:96034/1‑71 (MQ=255)
                          aGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGc                                           >  1:2207104/1‑71 (MQ=12)
                          aGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGc                                           >  1:454850/1‑71 (MQ=12)
                            cgcCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                                              >  1:4182131/1‑66 (MQ=255)
                            cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGccc                                                                 >  1:5450556/1‑45 (MQ=2)
                            cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGccc                                                                 >  1:1488991/1‑45 (MQ=2)
                            cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGccc                                                                 >  1:3311290/1‑45 (MQ=2)
                               cATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                                               >  1:2871521/1‑46 (MQ=255)
                               cATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGccg                                                   <  1:3185743/58‑2 (MQ=12)
                               cATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGccg                                                   <  1:1130810/58‑2 (MQ=12)
                                aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                     <  1:5433301/71‑1 (MQ=255)
                                aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                     <  1:5331884/71‑1 (MQ=255)
                                aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                     <  1:4161754/71‑1 (MQ=255)
                                 ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                                                 <  1:3313267/58‑1 (MQ=255)
                                 ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                     >  1:1600655/1‑70 (MQ=255)
                                 ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                     >  1:486730/1‑70 (MQ=255)
                                 ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt                                    <  1:3403910/71‑1 (MQ=255)
                                 ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt                                    <  1:64888/71‑1 (MQ=255)
                                  tGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCgga                                    >  1:1026478/1‑69 (MQ=9)
                                   gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                           <  1:5141432/46‑1 (MQ=255)
                                   gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                  >  1:1691408/1‑71 (MQ=255)
                                   gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                  >  1:5204418/1‑71 (MQ=255)
                                   gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                  >  1:5294934/1‑71 (MQ=255)
                                   gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                  >  1:2425457/1‑71 (MQ=255)
                                    cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                  <  1:4707146/70‑1 (MQ=255)
                                    cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                 <  1:4232480/71‑1 (MQ=255)
                                      gAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                               >  1:4973337/1‑71 (MQ=255)
                                       aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                 <  1:4185394/68‑1 (MQ=255)
                                       aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                 <  1:744121/68‑1 (MQ=255)
                                       aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                 <  1:771359/68‑1 (MQ=255)
                                          cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                                                         <  1:2282536/41‑1 (MQ=255)
                                             tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                        <  1:115177/71‑2 (MQ=9)
                                             tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                        <  1:4585766/71‑2 (MQ=9)
                                             tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                        <  1:4220291/71‑2 (MQ=9)
                                             tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                        >  1:4012000/1‑70 (MQ=9)
                                              ggTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                               <  1:2447143/47‑1 (MQ=255)
                                               gTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGa                      <  1:2636719/71‑1 (MQ=255)
                                                  aTAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                               >  1:3368717/1‑59 (MQ=255)
                                                  aTAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                   >  1:5367023/1‑71 (MQ=255)
                                                  aTAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                   >  1:4973468/1‑71 (MQ=255)
                                                    agagAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTa                 >  1:1350901/1‑71 (MQ=255)
                                                     gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                  >  1:5282410/1‑53 (MQ=255)
                                                      agaACTGACGATAGATAGCACCACTAGAACCCCCAaag                                                <  1:2297252/38‑1 (MQ=38)
                                                       gaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCa              <  1:1460261/71‑1 (MQ=255)
                                                       gaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCa              <  1:2545083/71‑1 (MQ=255)
                                                       gaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCa              <  1:4286376/71‑1 (MQ=255)
                                                        aaCTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCAt             <  1:3447088/71‑1 (MQ=9)
                                                        aaCTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCAt             <  1:950705/71‑1 (MQ=9)
                                                         aCTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCAt             <  1:2281518/70‑1 (MQ=9)
                                                          cTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcg           >  1:3512703/1‑71 (MQ=255)
                                                            gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCAt                                        >  1:753224/1‑40 (MQ=255)
                                                            gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCAt                                        >  1:4791495/1‑40 (MQ=255)
                                                            gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                     >  1:2983662/1‑43 (MQ=255)
                                                             gcGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGa        <  1:2261693/70‑1 (MQ=12)
                                                             aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                <  1:13856/63‑1 (MQ=255)
                                                             aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                <  1:4097488/63‑1 (MQ=255)
                                                             aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                <  1:4104221/63‑1 (MQ=255)
                                                             aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                <  1:2462337/63‑1 (MQ=255)
                                                             aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgcg         <  1:372901/70‑1 (MQ=255)
                                                             aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgcg         <  1:61448/70‑1 (MQ=255)
                                                             aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa        >  1:2999704/1‑71 (MQ=255)
                                                             aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa        >  1:3423818/1‑71 (MQ=255)
                                                               gatagataGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                   >  1:2907225/1‑42 (MQ=255)
                                                               gatagataGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                   >  1:3021514/1‑42 (MQ=255)
                                                               gataaattGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGATa      <  1:4363444/71‑1 (MQ=255)
                                                                   gatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGa  >  1:2544643/1‑71 (MQ=255)
                                                                   |                                                                      
TCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGA  >  minE/311408‑311545

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: