Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,697 C→T 14.4% T392T (ACG→ACA acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,6970CT14.4% 318.0 / 30.6 118T392T (ACG→ACAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (52/49);  new base T (8/9);  total (60/58)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

CGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAA  >  minE/311627‑311764
                                                                      |                                                                   
cGTTTTCTACCACAACGATGGCGTCATCCCCCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                                                     >  1:2161141/1‑71 (MQ=255)
    ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                                                                 >  1:5433983/1‑71 (MQ=255)
    ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                                                                 >  1:1397721/1‑71 (MQ=255)
    ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtg                                                                  <  1:4225786/70‑1 (MQ=255)
     tCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                                                                 <  1:3700339/70‑1 (MQ=255)
       caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                                              >  1:4149736/2‑69 (MQ=9)
       caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttg                                                               <  1:4994313/69‑2 (MQ=9)
         ccaccaCTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAtg                                                            <  1:477367/71‑2 (MQ=2)
          caccacTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAtg                                                            <  1:5562986/70‑2 (MQ=2)
           acaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCga                                                          <  1:2423164/71‑1 (MQ=255)
                  tGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                   <  1:1071670/71‑1 (MQ=255)
                   ggCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                   >  1:958283/1‑70 (MQ=255)
                   ggCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                   >  1:2725324/1‑70 (MQ=255)
                    gCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                   >  1:3622721/1‑69 (MQ=255)
                      gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                               <  1:2435204/71‑1 (MQ=255)
                      gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                               >  1:5049028/1‑71 (MQ=255)
                      gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                               >  1:5396240/1‑71 (MQ=255)
                       tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggc                                              <  1:4854868/71‑1 (MQ=255)
                       tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggc                                              <  1:1420939/71‑1 (MQ=255)
                         aTCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcgg                                            >  1:539611/1‑71 (MQ=255)
                          tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                  <  1:5550556/64‑1 (MQ=255)
                          tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                           <  1:5510265/71‑1 (MQ=255)
                          tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcgg                                            >  1:4525407/1‑70 (MQ=255)
                          tCCACCAACAGGCAGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                  <  1:4851666/64‑1 (MQ=255)
                            caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTGTCGAGAAGCCAAAg                                                >  1:4383075/1‑64 (MQ=255)
                            caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                                             >  1:2497605/1‑51 (MQ=255)
                            caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCga                                                          >  1:2975313/1‑54 (MQ=255)
                            caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAg                                                >  1:2174257/1‑64 (MQ=255)
                             accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                        <  1:3755918/71‑1 (MQ=255)
                             accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                        <  1:3421285/71‑1 (MQ=255)
                             accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                        <  1:2964428/71‑1 (MQ=255)
                               gagcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagc                                              <  1:600496/60‑1 (MQ=9)
                               caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                                             <  1:5616912/48‑1 (MQ=255)
                               caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                                             <  1:540908/48‑1 (MQ=255)
                               caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                                             <  1:4071180/48‑1 (MQ=255)
                               caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCga                                                          <  1:510460/51‑1 (MQ=255)
                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGc                                                     >  1:2718452/1‑55 (MQ=255)
                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                                                 >  1:793641/1‑59 (MQ=255)
                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                                                 >  1:1893672/1‑59 (MQ=255)
                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                   >  1:5318445/1‑57 (MQ=255)
                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                     >  1:469348/1‑71 (MQ=255)
                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                     >  1:3244469/1‑71 (MQ=255)
                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                     >  1:4240022/1‑71 (MQ=255)
                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                     >  1:5483622/1‑71 (MQ=255)
                                 gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAgtgt                                                                   <  1:161080/39‑1 (MQ=18)
                                 aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                        <  1:5028315/67‑1 (MQ=255)
                                   aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                                  <  1:841289/71‑1 (MQ=255)
                                    ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                 <  1:572928/71‑1 (MQ=255)
                                    ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                 <  1:4191323/71‑1 (MQ=255)
                                     cccGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCaa                                <  1:1437865/70‑1 (MQ=16)
                                      ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                           <  1:3557547/59‑1 (MQ=255)
                                       cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                              <  1:1094316/71‑1 (MQ=255)
                                        gATGGCTAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                    <  1:3399714/48‑1 (MQ=255)
                                        gATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                    <  1:1626237/48‑1 (MQ=255)
                                        gATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                    <  1:3037034/48‑1 (MQ=255)
                                         aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                           <  1:3707842/56‑1 (MQ=255)
                                         aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                           <  1:4838550/56‑1 (MQ=255)
                                         aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                           <  1:219989/56‑1 (MQ=255)
                                         aTGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGt                            <  1:1344227/71‑1 (MQ=12)
                                          tGGCGAGCACCATCCCTAACATTGTTATCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                 <  1:4981848/65‑1 (MQ=255)
                                           ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga                                                        >  1:4550352/1‑41 (MQ=255)
                                            gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                         <  1:4340837/71‑1 (MQ=255)
                                            gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                         <  1:1865527/71‑1 (MQ=255)
                                            gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                         <  1:1676517/71‑1 (MQ=255)
                                              aaGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTccct                        >  1:85032/2‑69 (MQ=12)
                                               aGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                  <  1:3260133/43‑1 (MQ=255)
                                               aGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                  <  1:646078/43‑1 (MQ=255)
                                               aGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAg                      >  1:2721745/1‑71 (MQ=255)
                                               aGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                     >  1:4923939/1‑55 (MQ=12)
                                               aGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTccctaa                      >  1:1227868/1‑68 (MQ=12)
                                                gCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAg                      >  1:3962273/1‑70 (MQ=255)
                                                 cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                          >  1:4627673/1‑65 (MQ=255)
                                                 cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                     >  1:1881803/1‑70 (MQ=255)
                                                  aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                      >  1:68116/1‑52 (MQ=255)
                                                  aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                      >  1:785974/1‑52 (MQ=255)
                                                  aCCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaa                   <  1:1302844/71‑1 (MQ=12)
                                                   ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCa                    >  1:2161714/1‑69 (MQ=255)
                                                   ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCa                    >  1:2964362/1‑69 (MQ=255)
                                                   ccATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaac                  <  1:1962639/71‑2 (MQ=12)
                                                    cATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTc                           >  1:128890/1‑61 (MQ=12)
                                                    cATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaaca                 >  1:4166017/1‑69 (MQ=12)
                                                      tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcgg                                            >  1:3876650/1‑42 (MQ=255)
                                                      tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCaa                                         >  1:2579569/1‑45 (MQ=255)
                                                      tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                <  1:4952290/54‑1 (MQ=255)
                                                      tCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaaca                 >  1:91188/1‑67 (MQ=12)
                                                        ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                     <  1:4266106/63‑1 (MQ=255)
                                                          gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg           >  1:4877978/1‑71 (MQ=255)
                                                           aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg          <  1:1693640/71‑1 (MQ=255)
                                                            aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg          >  1:4653480/1‑70 (MQ=255)
                                                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                      >  1:2130522/1‑41 (MQ=255)
                                                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                      >  1:3814848/1‑41 (MQ=255)
                                                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                 <  1:1352518/46‑1 (MQ=255)
                                                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg           >  1:1447142/1‑68 (MQ=255)
                                                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGAGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg           >  1:1376820/1‑68 (MQ=255)
                                                                 gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                         >  1:1032302/1‑50 (MQ=255)
                                                                 gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                         >  1:727980/1‑50 (MQ=255)
                                                                 gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                         >  1:1202656/1‑50 (MQ=255)
                                                                 gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACgg     >  1:1724965/1‑70 (MQ=255)
                                                                  ttAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                      <  1:4969280/36‑1 (MQ=255)
                                                                      cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                <  1:800495/38‑1 (MQ=255)
                                                                      |                                                                   
CGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAA  >  minE/311627‑311764

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: