Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,445 T→C 18.7% S476S (TCA→TCG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4450TC18.7% 295.2 / 50.3 123S476S (TCA→TCGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (44/56);  new base C (8/15);  total (52/71)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 124 total aligned reads displayed.

TCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACG  >  minE/311375‑311514
                                                                      |                                                                     
  ggTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGa                                                                     <  1:5086730/71‑1 (MQ=255)
  ggTTTAAGCAGCGTTGCACATAACGCAGGGGTAAGAATCAATGCCACCAGAACAGAAAGCGCCATTGCCGa                                                                     <  1:1040737/71‑1 (MQ=255)
  ggTTTAAGCAGCGTTGCACATAACGCAGGGGTAAGAATCAATGCCACCAGAACAGAAAGCGCCATTGCCGa                                                                     <  1:567275/71‑1 (MQ=255)
   gTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGaa                                                                    <  1:4773737/71‑1 (MQ=255)
     ttCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGaaa                                                                   >  1:1444453/1‑70 (MQ=255)
      tCAGCATGGTGGCACAAAGCGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                                  >  1:2331511/1‑70 (MQ=255)
       cAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACaa                                                                >  1:3589307/1‑71 (MQ=255)
       cAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACa                                                                 <  1:4875638/70‑1 (MQ=255)
       cAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACa                                                                 <  1:1842904/70‑1 (MQ=255)
        aGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACaa                                                                <  1:3463187/70‑1 (MQ=255)
        aGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCAACAGTACCGACAGCGCCATTGCTGAAACaa                                                                <  1:4230990/70‑1 (MQ=255)
           aTGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGt                                                            <  1:4324341/71‑1 (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:5594981/71‑1 (MQ=255)
                        gagCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGa                                               >  1:1974332/1‑71 (MQ=255)
                         agCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                              <  1:156881/71‑1 (MQ=255)
                          gCTGGAGTCAGGATCAACGCCACCAGTACCTACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                              <  1:303446/70‑1 (MQ=255)
                           cTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                              >  1:2590486/1‑69 (MQ=255)
                              gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACt                                                 <  1:688034/63‑1 (MQ=255)
                              gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgata                                          <  1:4413214/70‑1 (MQ=255)
                              gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgata                                          <  1:2513684/70‑1 (MQ=255)
                               aGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                                            >  1:1595649/1‑67 (MQ=255)
                                 tCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                        <  1:4649763/53‑1 (MQ=255)
                                 tCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                        <  1:302890/53‑1 (MQ=255)
                                 tCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                                        <  1:5185310/53‑1 (MQ=255)
                                 tCAGGATCAACGCCACCAGTACAGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgatagata                                      <  1:4491167/71‑1 (MQ=255)
                                   aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag                                     >  1:2271140/1‑70 (MQ=255)
                                     gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCa                                   >  1:5493319/1‑70 (MQ=255)
                                     gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                  >  1:4386025/1‑71 (MQ=255)
                                     gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                  >  1:2687123/1‑71 (MQ=255)
                                     gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                                  >  1:5538500/1‑71 (MQ=255)
                                      aTCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                                              <  1:1375858/58‑1 (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)
                                       tCAACGCCACCAGTACCGACAGCGCCATTGCTCAAACAATGGTAATAGAGAACTGAc                                              <  1:1154881/57‑1 (MQ=255)
                                         aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGt                              <  1:2056951/68‑1 (MQ=9)
                                         aaCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg                               <  1:483869/70‑1 (MQ=255)
                                          aCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGaa                                                   >  1:4705125/1‑49 (MQ=255)
                                          aCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGaa                                                   >  1:2156626/1‑49 (MQ=255)
                                          aCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGaa                                                   >  1:3889680/1‑49 (MQ=255)
                                           cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGt                                                            <  1:721919/39‑1 (MQ=255)
                                           cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGaa                                                   >  1:1789062/1‑48 (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)
                                                  aGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGc                      >  1:1819492/1‑70 (MQ=255)
                                                   gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                    <  1:1905006/71‑1 (MQ=255)
                                                   gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                    <  1:93049/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)
                                                       cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgatag                                         <  1:3576800/46‑1 (MQ=255)
                                                       cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                <  1:2455357/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:3393074/68‑1 (MQ=12)
                                                        gACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag                 >  1:1456554/1‑69 (MQ=255)
                                                         aaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCgata                                          <  1:4078824/41‑1 (MQ=14)
                                                         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:454850/1‑71 (MQ=12)
                                                           aGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGc            >  1:2207104/1‑71 (MQ=12)
                                                             cgcCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa               >  1:9675/1‑66 (MQ=255)
                                                             cgcCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa               >  1:4182131/1‑66 (MQ=255)
                                                             cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGccc                                  >  1:1488991/1‑45 (MQ=2)
                                                             cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGccc                                  >  1:3311290/1‑45 (MQ=2)
                                                             cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGccc                                  >  1:5450556/1‑45 (MQ=2)
                                                             cgcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGcca                                  >  1:4314905/1‑45 (MQ=2)
                                                                cATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                >  1:2991884/1‑46 (MQ=255)
                                                                cATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa                                >  1:2871521/1‑46 (MQ=255)
                                                                cATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGccg                    <  1:1130810/58‑2 (MQ=12)
                                                                cATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGccg                    <  1:3185743/58‑2 (MQ=12)
                                                                 aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg      <  1:4161754/71‑1 (MQ=255)
                                                                 aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg      <  1:5331884/71‑1 (MQ=255)
                                                                  ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg      >  1:1600655/1‑70 (MQ=255)
                                                                  ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg      >  1:486730/1‑70 (MQ=255)
                                                                  ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt     <  1:5332686/71‑1 (MQ=255)
                                                                  ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt     <  1:64888/71‑1 (MQ=255)
                                                                  ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt     <  1:3403910/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:5294934/1‑71 (MQ=255)
                                                                    gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc   >  1:5204418/1‑71 (MQ=255)
                                                                    gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc   >  1:1691408/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)
                                                                      |                                                                     
TCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACG  >  minE/311375‑311514

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: