Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,673 G→A 23.1% L400L (CTC→CTT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,6730GA23.1% 231.9 / 55.2 100L400L (CTC→CTTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (35/41);  new base A (15/8);  total (51/49)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 101 total aligned reads displayed.

CTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCC  >  minE/311603‑311741
                                                                      |                                                                    
cttcttCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCg                                                                      <  1:123282/71‑1 (MQ=255)
 ttcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGa                                                                     >  1:3515405/1‑71 (MQ=255)
  tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                                                    >  1:488728/1‑71 (MQ=255)
  tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                                                    >  1:3675685/1‑71 (MQ=255)
      ccGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCa                                                                  <  1:1651762/69‑1 (MQ=255)
       cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCAcc                                                                <  1:274408/70‑1 (MQ=255)
       cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCAcc                                                                <  1:903897/70‑1 (MQ=255)
       cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCAcc                                                                <  1:3713951/70‑1 (MQ=255)
          cATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATcc                                                            >  1:333152/1‑71 (MQ=255)
              acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                                        >  1:2133125/1‑71 (MQ=255)
              acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGa                                                         >  1:3355033/1‑70 (MQ=255)
               cacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                                        >  1:3239778/1‑70 (MQ=255)
                acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATc                                                             <  1:3175583/64‑1 (MQ=255)
                acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATc                                                             <  1:3904518/64‑1 (MQ=255)
                acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                                      <  1:1348260/71‑1 (MQ=255)
                 cGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAt                                                     <  1:1507946/71‑1 (MQ=255)
                 cGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAt                                                     >  1:3003479/1‑71 (MQ=255)
                   cTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAtt                                                    <  1:1341490/70‑1 (MQ=255)
                      aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGa                                                                     >  1:1735466/1‑50 (MQ=255)
                      aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGtt                                                 <  1:2023695/70‑1 (MQ=255)
                      aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTa                                                <  1:692119/71‑1 (MQ=255)
                        cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                              <  1:848997/71‑1 (MQ=255)
                        cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                              >  1:1053744/1‑71 (MQ=255)
                         gTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCg                                             <  1:2724515/71‑1 (MQ=255)
                          ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAg                                               <  1:1369796/68‑1 (MQ=255)
                            ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                                          >  1:2238548/1‑71 (MQ=255)
                            ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                                          >  1:282127/1‑71 (MQ=255)
                             tCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                         >  1:3860039/1‑71 (MQ=255)
                               tACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                       >  1:3897047/1‑71 (MQ=255)
                               tACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                       >  1:1699142/1‑71 (MQ=255)
                               caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                       >  1:93297/2‑69 (MQ=9)
                               caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                       >  1:3350701/2‑69 (MQ=9)
                               caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                       >  1:482721/2‑69 (MQ=9)
                               caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttg                                        <  1:193640/69‑2 (MQ=9)
                               caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttg                                        <  1:1540903/69‑2 (MQ=9)
                                accaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAt                                      <  1:1366210/71‑1 (MQ=9)
                                accaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAt                                      >  1:2544236/1‑71 (MQ=9)
                                aCCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                      >  1:3452488/1‑71 (MQ=255)
                                 ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                      <  1:1304492/70‑1 (MQ=255)
                                 ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATc                                     <  1:1727597/71‑1 (MQ=255)
                                  cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                                      <  1:3109005/53‑1 (MQ=255)
                                  cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                       <  1:1515133/68‑1 (MQ=255)
                                  cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                                    <  1:564412/71‑1 (MQ=255)
                                  cacaACGATCGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                                    <  1:1458078/71‑1 (MQ=255)
                                    caACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag                                  <  1:2228932/71‑1 (MQ=255)
                                     aaCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATccc                                                           >  1:1344730/1‑45 (MQ=255)
                                     aaCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga                                 <  1:2189356/71‑1 (MQ=255)
                                          tGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                             <  1:2971964/70‑1 (MQ=255)
                                             cGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAg                         <  1:3797414/71‑1 (MQ=255)
                                              atCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaaaag                        >  1:2778838/2‑69 (MQ=9)
                                              atCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaaaag                        >  1:3276913/2‑69 (MQ=9)
                                               tCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGtt                                         <  1:1540298/53‑1 (MQ=16)
                                               tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                                    <  1:1899985/58‑1 (MQ=255)
                                                cATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAg                               >  1:227123/1‑62 (MQ=255)
                                                 aTCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                          >  1:563275/1‑66 (MQ=255)
                                                  tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                    <  1:2575103/71‑1 (MQ=255)
                                                  tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                    <  1:1824706/71‑1 (MQ=255)
                                                      ccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                         >  1:3401248/1‑46 (MQ=255)
                                                       caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa               <  1:2470315/71‑1 (MQ=255)
                                                        agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagcggc                    <  1:1648596/63‑1 (MQ=12)
                                                        agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat              >  1:2319320/3‑70 (MQ=12)
                                                        agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat              <  1:2299073/69‑2 (MQ=12)
                                                        aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                              >  1:1262751/1‑39 (MQ=255)
                                                        aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                              >  1:426267/1‑39 (MQ=255)
                                                        aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                                      <  1:3787862/47‑1 (MQ=255)
                                                        aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGaa                                >  1:2471585/1‑53 (MQ=255)
                                                        aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc              <  1:2742190/71‑1 (MQ=255)
                                                         gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGatg             >  1:576654/2‑69 (MQ=12)
                                                         gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGatg             >  1:3815727/2‑69 (MQ=12)
                                                         aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                         >  1:2683269/1‑43 (MQ=255)
                                                         aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                         >  1:1569027/1‑43 (MQ=255)
                                                         aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg             >  1:3497243/1‑71 (MQ=255)
                                                           aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACgg            >  1:3534739/1‑70 (MQ=255)
                                                            ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGc                              <  1:3191325/51‑1 (MQ=255)
                                                             gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa         <  1:1796073/71‑1 (MQ=255)
                                                             gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa         <  1:3864917/71‑1 (MQ=255)
                                                             gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa         <  1:2363271/71‑1 (MQ=255)
                                                             gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa         <  1:2124814/71‑1 (MQ=255)
                                                             gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa          <  1:3006533/70‑1 (MQ=255)
                                                             cccGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCa          <  1:690757/69‑1 (MQ=12)
                                                              ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga                                 <  1:899035/46‑1 (MQ=255)
                                                              ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa          <  1:1565651/69‑1 (MQ=255)
                                                               cGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCaaac       <  1:3652167/71‑2 (MQ=16)
                                                                gATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                           >  1:1247340/1‑50 (MQ=9)
                                                                gATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                           >  1:892493/1‑50 (MQ=9)
                                                                gATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                           >  1:2075371/1‑50 (MQ=9)
                                                                 aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt     <  1:3732859/71‑1 (MQ=255)
                                                                 aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt     >  1:1668109/1‑71 (MQ=255)
                                                                 aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt     >  1:820374/1‑71 (MQ=255)
                                                                 aTGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAt                                      >  1:3094098/1‑38 (MQ=17)
                                                                  tGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTc    >  1:1306873/1‑71 (MQ=255)
                                                                  tGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTc    >  1:2077061/1‑71 (MQ=255)
                                                                  tGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagcgg                     >  1:2130779/1‑54 (MQ=9)
                                                                  tGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagcgg                     >  1:3950389/1‑54 (MQ=9)
                                                                   ggCTAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                        >  1:192406/1‑50 (MQ=255)
                                                                   ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc   >  1:3435692/1‑71 (MQ=255)
                                                                    gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa               >  1:2907973/1‑58 (MQ=255)
                                                                    gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg             >  1:762918/1‑60 (MQ=255)
                                                                    gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc  <  1:2262432/71‑1 (MQ=255)
                                                                    gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc  <  1:3526958/71‑1 (MQ=255)
                                                                      |                                                                    
CTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCC  >  minE/311603‑311741

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: