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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4570AG19.2% 299.3 / 53.0 125I472I (ATT→ATCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (43/58);  new base G (9/15);  total (52/73)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 126 total aligned reads displayed.

GGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGA  >  minE/311388‑311525
                                                                     |                                                                    
ggtggCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTaa                                                                     >  1:1744286/1‑71 (MQ=255)
 gtggCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAAt                                                                    <  1:1732544/71‑1 (MQ=255)
  tggCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                                   <  1:1762660/71‑1 (MQ=255)
    gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                                   >  1:4996799/1‑69 (MQ=255)
         aaGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAACGCCATTGCTGAAACAATGGTAATAGAGAACt                                                            >  1:3903711/1‑71 (MQ=255)
          agagCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTg                                                           <  1:4718153/71‑1 (MQ=255)
           gagCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGa                                                          <  1:4798511/71‑1 (MQ=255)
           gagCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGa                                                          <  1:5594981/71‑1 (MQ=255)
            agCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                         <  1:156881/71‑1 (MQ=255)
             gCTGGAGTCAGGATCAACGCCACCAGTACCTACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                         <  1:303446/70‑1 (MQ=255)
             gCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                         <  1:4281210/70‑1 (MQ=255)
              cTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                         >  1:2590486/1‑69 (MQ=255)
                ggAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                                                       >  1:1241653/1‑69 (MQ=255)
                ggAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                                                       >  1:4453816/1‑69 (MQ=255)
                 gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACt                                                            <  1:688034/63‑1 (MQ=255)
                 gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgata                                                     <  1:4413214/70‑1 (MQ=255)
                 gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgata                                                     <  1:1889252/70‑1 (MQ=255)
                  aGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                                                       >  1:1595649/1‑67 (MQ=255)
                    tCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                                   <  1:302890/53‑1 (MQ=255)
                    tCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                                   <  1:5185310/53‑1 (MQ=255)
                    tCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                                   <  1:4649763/53‑1 (MQ=255)
                    tCAGGATCAACGCCACCAGTACAGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgatagata                                                 <  1:4491167/71‑1 (MQ=255)
                      agaaTCAATGCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGc                                               <  1:4657457/68‑1 (MQ=9)
                      aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                         <  1:3808007/61‑1 (MQ=255)
                      aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag                                                >  1:2271140/1‑70 (MQ=255)
                        gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCa                                              >  1:5493319/1‑70 (MQ=255)
                        gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                             >  1:5080207/1‑71 (MQ=255)
                        gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                             >  1:5538500/1‑71 (MQ=255)
                        gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                             >  1:2687123/1‑71 (MQ=255)
                         aTCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                                         <  1:1375858/58‑1 (MQ=255)
                          tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgataga                                                   >  1:1338692/1‑63 (MQ=255)
                          tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgataga                                                   >  1:2619024/1‑63 (MQ=255)
                          tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                             >  1:3183534/1‑69 (MQ=255)
                           cAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTATTAGAGAACTGACGATAGATAGCACCAg                                          >  1:3482320/1‑71 (MQ=255)
                            aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGt                                         <  1:2056951/68‑1 (MQ=9)
                            aaCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg                                          <  1:483869/70‑1 (MQ=255)
                             aCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGaa                                                              >  1:2156626/1‑49 (MQ=255)
                              cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGaa                                                              >  1:2794765/1‑48 (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:1538083/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)
                                      gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                               <  1:1905006/71‑1 (MQ=255)
                                       tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgataga                                                   >  1:3704843/1‑50 (MQ=255)
                                         ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                       <  1:1057858/60‑1 (MQ=255)
                                         ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                             <  1:3573257/70‑1 (MQ=255)
                                          cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgatag                                                    <  1:2221125/46‑1 (MQ=255)
                                          cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                           <  1:669992/71‑1 (MQ=255)
                                          cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                           <  1:1989027/71‑1 (MQ=255)
                                          cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                           <  1:2455357/71‑1 (MQ=255)
                                           gaaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCTATAAATTGCCCCAGTAGAACCGCCgaagaa                          <  1:3210604/68‑1 (MQ=12)
                                           gaaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCgaagaa                          <  1:2556783/68‑1 (MQ=12)
                                           gACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag                            >  1:1456554/1‑69 (MQ=255)
                                            aaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCgata                                                     <  1:4078824/41‑1 (MQ=14)
                                            aCAGCGCCATTGCTGAAACAATGGTAATGGAGAACTGACGATAGATAGCAcc                                            <  1:377635/52‑1 (MQ=255)
                                            aCAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaag                         >  1:3356927/1‑71 (MQ=255)
                                             cAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAgg                        >  1:96034/1‑71 (MQ=255)
                                              aGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGc                       >  1:454850/1‑71 (MQ=12)
                                                cgcCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                          >  1:4182131/1‑66 (MQ=255)
                                                cgcCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                          >  1:9675/1‑66 (MQ=255)
                                                cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGccc                                             >  1:1488991/1‑45 (MQ=2)
                                                cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGccc                                             >  1:5450556/1‑45 (MQ=2)
                                                cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGccc                                             >  1:3311290/1‑45 (MQ=2)
                                                cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGcca                                             >  1:4314905/1‑45 (MQ=2)
                                                   cATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                           >  1:2871521/1‑46 (MQ=255)
                                                   cATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGccg                               <  1:1130810/58‑2 (MQ=12)
                                                    aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                 <  1:4161754/71‑1 (MQ=255)
                                                    aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                 <  1:5331884/71‑1 (MQ=255)
                                                     ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                             <  1:3313267/58‑1 (MQ=255)
                                                     ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                 >  1:486730/1‑70 (MQ=255)
                                                     ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt                <  1:3403910/71‑1 (MQ=255)
                                                     ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt                <  1:64888/71‑1 (MQ=255)
                                                     ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt                <  1:5332686/71‑1 (MQ=255)
                                                      tGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCgga                >  1:1026478/1‑69 (MQ=9)
                                                       gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg                                       <  1:5141432/46‑1 (MQ=255)
                                                       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)
                                                       gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc              >  1:1691408/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)
                                                           aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                          <  1:2727620/55‑1 (MQ=255)
                                                           aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg             <  1:744121/68‑1 (MQ=255)
                                                           aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg             <  1:4185394/68‑1 (MQ=255)
                                                             aCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                 >  1:2289126/1‑62 (MQ=255)
                                                              cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAcc       <  1:45062/71‑1 (MQ=255)
                                                                aTGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAccgc     >  1:4106045/1‑71 (MQ=255)
                                                                 tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct    <  1:115177/71‑2 (MQ=9)
                                                                 tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct    >  1:4012000/1‑70 (MQ=9)
                                                                 tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct    <  1:4585766/71‑2 (MQ=9)
                                                                 tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct    <  1:4220291/71‑2 (MQ=9)
                                                                  ggTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                           <  1:2447143/47‑1 (MQ=255)
                                                                    tAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                 <  1:400823/55‑1 (MQ=255)
                                                                     |                                                                    
GGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGA  >  minE/311388‑311525

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: