Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,703 T→G 12.2% I390I (ATA→ATC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7030TG12.2% 309.1 / 33.6 123I390I (ATA→ATCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (55/53);  new base G (7/8);  total (62/61)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

CCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGG  >  minE/311636‑311771
                                                                   |                                                                    
ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATc                                                                   >  1:1276196/1‑71 (MQ=255)
 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)
           gcatCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                                         >  1:1607906/4‑70 (MQ=9)
           gCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                          >  1:3622721/1‑69 (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)
                 tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                         <  1:2642586/64‑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:2964428/71‑1 (MQ=255)
                    accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                               <  1:3421285/71‑1 (MQ=255)
                    accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                               <  1:3755918/71‑1 (MQ=255)
                      gagcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagc                                                     <  1:600496/60‑1 (MQ=9)
                      gagcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGa                                             <  1:3877372/68‑1 (MQ=12)
                      caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                                                    <  1:4071180/48‑1 (MQ=255)
                      caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCga                                                                 <  1:510460/51‑1 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGc                                                            >  1:2718452/1‑55 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                                                        >  1:497563/1‑59 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                                                        >  1:1893672/1‑59 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                                                        >  1:793641/1‑59 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                          >  1:5318445/1‑57 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                          >  1:560383/1‑57 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                            >  1:4240022/1‑71 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                            >  1:3244469/1‑71 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                            >  1:5483622/1‑71 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                            >  1:469348/1‑71 (MQ=255)
                           ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                        <  1:4191323/71‑1 (MQ=255)
                           ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                        <  1:572928/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:82209/48‑1 (MQ=255)
                               gATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                           <  1:3037034/48‑1 (MQ=255)
                               gATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                           <  1:1626237/48‑1 (MQ=255)
                               gATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                           <  1:5422774/48‑1 (MQ=255)
                                aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                  <  1:3707842/56‑1 (MQ=255)
                                aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                  <  1:219989/56‑1 (MQ=255)
                                aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                  <  1:4838550/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)
                                  ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                      >  1:5263517/1‑66 (MQ=255)
                                   gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                <  1:1865527/71‑1 (MQ=255)
                                   gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                <  1:4340837/71‑1 (MQ=255)
                                   gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                <  1:3276581/71‑1 (MQ=255)
                                   gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                <  1:1676517/71‑1 (MQ=255)
                                     aaGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTccct                               >  1:85032/2‑69 (MQ=12)
                                      aGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                         <  1:646078/43‑1 (MQ=255)
                                      aGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                            >  1:4923939/1‑55 (MQ=12)
                                       gCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAg                             >  1:3962273/1‑70 (MQ=255)
                                        cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                            >  1:1881803/1‑70 (MQ=255)
                                         aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                             >  1:68116/1‑52 (MQ=255)
                                         aCCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaa                          <  1:1302844/71‑1 (MQ=12)
                                          ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCa                           >  1:2964362/1‑69 (MQ=255)
                                          ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAAt                         <  1:4201267/71‑1 (MQ=255)
                                           cATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                            <  1:1733292/67‑1 (MQ=255)
                                           cATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTc                                  >  1:128890/1‑61 (MQ=12)
                                             tCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaaca                        >  1:91188/1‑67 (MQ=12)
                                               ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                            <  1:4266106/63‑1 (MQ=255)
                                                 gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                  >  1:4877978/1‑71 (MQ=255)
                                                  aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                            <  1:2224202/60‑1 (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)
                                                    cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGAGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                  >  1:1376820/1‑68 (MQ=255)
                                                     aTTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTAc              <  1:4654934/71‑1 (MQ=255)
                                                        gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                >  1:1032302/1‑50 (MQ=255)
                                                        gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                >  1:4456264/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)
                                                          tAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCaa         >  1:196454/1‑71 (MQ=255)
                                                             cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                       <  1:800495/38‑1 (MQ=255)
                                                              gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                     >  1:2100296/1‑39 (MQ=255)
                                                              gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                     >  1:753075/1‑39 (MQ=255)
                                                              gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                     >  1:5462723/1‑39 (MQ=255)
                                                                gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                 <  1:2718617/57‑1 (MQ=255)
                                                                gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCg   >  1:3594483/1‑71 (MQ=255)
                                                                 tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                      >  1:2547344/1‑35 (MQ=255)
                                                                 tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                      >  1:2446163/1‑35 (MQ=255)
                                                                 tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                      >  1:5407435/1‑35 (MQ=255)
                                                                 tttATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                 <  1:141612/56‑1 (MQ=255)
                                                                 tttATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCgg  <  1:494426/71‑1 (MQ=255)
                                                                   tATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                <  1:4050109/39‑1 (MQ=255)
                                                                   |                                                                    
CCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGG  >  minE/311636‑311771

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: