Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,481 A→C 16.9% G464G (GGT→GGG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4810AC16.9% 260.8 / 50.2 118G464G (GGT→GGGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (42/56);  new base C (5/15);  total (47/71)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 119 total aligned reads displayed.

GATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCA  >  minE/311412‑311548
                                                                     |                                                                   
gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCa                                                                     >  1:5493319/1‑70 (MQ=255)
gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                    >  1:2687123/1‑71 (MQ=255)
gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                    >  1:4386025/1‑71 (MQ=255)
gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                    >  1:5080207/1‑71 (MQ=255)
gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                    >  1:5538500/1‑71 (MQ=255)
  tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                    >  1:3183534/1‑69 (MQ=255)
  tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                                                    >  1:454634/1‑69 (MQ=255)
   cAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTATTAGAGAACTGACGATAGATAGCACCAg                                                                 >  1:3482320/1‑71 (MQ=255)
      cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                              <  1:5619597/71‑1 (MQ=255)
      cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                              <  1:5271663/71‑1 (MQ=255)
       gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGa                                                             <  1:228216/71‑1 (MQ=255)
        ccaccaGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                                                            >  1:1947621/1‑71 (MQ=255)
         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)
              gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                                                      <  1:1905006/71‑1 (MQ=255)
                 ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                              <  1:1057858/60‑1 (MQ=255)
                 ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                                                    <  1:3573257/70‑1 (MQ=255)
                  cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                                  <  1:1989027/71‑1 (MQ=255)
                  cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                                  <  1:1477252/71‑1 (MQ=255)
                  cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                                  <  1:2455357/71‑1 (MQ=255)
                  cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                                  <  1:669992/71‑1 (MQ=255)
                   gaaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCTATAAATTGCCCCAGTAGAACCGCCgaagaa                                                 <  1:3210604/68‑1 (MQ=12)
                   gaaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCgaagaa                                                 <  1:3393074/68‑1 (MQ=12)
                   gaaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCgaagaa                                                 <  1:2556783/68‑1 (MQ=12)
                   gACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag                                                   >  1:1456554/1‑69 (MQ=255)
                    aaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCgaagaat                                                >  1:4834798/3‑70 (MQ=12)
                    aCAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaag                                                >  1:3356927/1‑71 (MQ=255)
                     cAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAgg                                               >  1:96034/1‑71 (MQ=255)
                      aGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGc                                              >  1:2207104/1‑71 (MQ=12)
                        cgcCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                                                 >  1:4182131/1‑66 (MQ=255)
                           cATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                                                  >  1:2871521/1‑46 (MQ=255)
                           cATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                                                  >  1:2991884/1‑46 (MQ=255)
                           cATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGccg                                                      <  1:1130810/58‑2 (MQ=12)
                           cATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGccg                                                      <  1:3185743/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:64888/71‑1 (MQ=255)
                             ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt                                       <  1:5332686/71‑1 (MQ=255)
                             ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt                                       <  1:3403910/71‑1 (MQ=255)
                              tGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCgga                                       >  1:1026478/1‑69 (MQ=9)
                               gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                                              <  1:260148/46‑1 (MQ=255)
                               gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                     >  1:2425457/1‑71 (MQ=255)
                               gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                     >  1:5204418/1‑71 (MQ=255)
                               gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                     >  1:5294934/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:5077240/68‑1 (MQ=255)
                                   aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                    <  1:771359/68‑1 (MQ=255)
                                   aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                    <  1:744121/68‑1 (MQ=255)
                                   aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                    <  1:4191929/68‑1 (MQ=255)
                                   aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                    <  1:4185394/68‑1 (MQ=255)
                                      cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGaa                                                            <  1:2282536/41‑1 (MQ=255)
                                      cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAcc                              <  1:45062/71‑1 (MQ=255)
                                         tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                           <  1:4220291/71‑2 (MQ=9)
                                         tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                           >  1:4012000/1‑70 (MQ=9)
                                         tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                           <  1:115177/71‑2 (MQ=9)
                                         tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                           <  1:4585766/71‑2 (MQ=9)
                                          ggTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                                  <  1:2447143/47‑1 (MQ=255)
                                           gTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGa                         <  1:2636719/71‑1 (MQ=255)
                                            tAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                        <  1:400823/55‑1 (MQ=255)
                                              aTAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                                  >  1:3368717/1‑59 (MQ=255)
                                              aTAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                      >  1:4973468/1‑71 (MQ=255)
                                               tAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACa                       <  1:679377/69‑1 (MQ=255)
                                                agagAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGa                         <  1:788602/66‑1 (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)
                                                    aaCTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCAt                <  1:3447088/71‑1 (MQ=9)
                                                    aaCTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCAt                <  1:950705/71‑1 (MQ=9)
                                                     aCTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCAt                <  1:2281518/70‑1 (MQ=9)
                                                     aCTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                      <  1:1229471/48‑1 (MQ=255)
                                                      cTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcg              >  1:3512703/1‑71 (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:2462337/63‑1 (MQ=255)
                                                         aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                   <  1:4097488/63‑1 (MQ=255)
                                                         aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                   <  1:4104221/63‑1 (MQ=255)
                                                         aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgcg            <  1:372901/70‑1 (MQ=255)
                                                         aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgcg            <  1:61448/70‑1 (MQ=255)
                                                         aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc             <  1:4589323/69‑1 (MQ=255)
                                                         aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa           >  1:3048474/1‑71 (MQ=255)
                                                           gatagataGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                      >  1:2097921/1‑42 (MQ=255)
                                                           gatagataGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                      >  1:2907225/1‑42 (MQ=255)
                                                           gatagataGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                      >  1:3021514/1‑42 (MQ=255)
                                                           gatagataGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGAt          <  1:5466771/70‑1 (MQ=255)
                                                           gataaattGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGATa         <  1:4363444/71‑1 (MQ=255)
                                                                  agCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGAc    >  1:2910826/1‑69 (MQ=255)
                                                                   gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCa  <  1:4451474/70‑1 (MQ=255)
                                                                   gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCa  <  1:1459092/70‑1 (MQ=255)
                                                                     |                                                                   
GATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCA  >  minE/311412‑311548

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: