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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7720AG15.9% 277.8 / 33.9 107I367I (ATT→ATCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (42/48);  new base G (10/7);  total (52/55)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

ATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCA  >  minE/311704‑311842
                                                                    |                                                                      
aTCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAAt                                                                      >  1:2996679/1‑71 (MQ=255)
 tCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATc                                                                     <  1:2303522/71‑1 (MQ=255)
    agaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCa                                                                    <  1:208985/69‑1 (MQ=255)
     gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                                                                 <  1:2329970/71‑1 (MQ=255)
     gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                                                                 <  1:113161/71‑1 (MQ=255)
      aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAAc                                                                  >  1:2228641/1‑69 (MQ=255)
      aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGt                                                                >  1:2290489/1‑71 (MQ=255)
      aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAAAACCACCGGTACGGCAATGGTCGGAATCaa                                                                   >  1:2715058/1‑68 (MQ=255)
          cAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcg                                                            >  1:2071402/1‑71 (MQ=255)
             aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcgcg                                                          >  1:4846186/1‑70 (MQ=255)
             aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCgg                                                         <  1:3663923/71‑1 (MQ=255)
              ggcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGa                                                        <  1:2655864/71‑1 (MQ=255)
               gcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCACCGTCGCGCgg                                                         >  1:3497152/1‑69 (MQ=255)
                   cAAGGACGGCAAAGGGCCCGAGCAATACCACCGGTACGGCAATGGTCGGa                                                                        <  1:4347077/50‑1 (MQ=255)
                    aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt                                                  <  1:1388741/71‑1 (MQ=255)
                     aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                 >  1:2930044/1‑71 (MQ=255)
                     aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                 >  1:2082518/1‑71 (MQ=255)
                     aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGGCGGAATCAACGTCGCGCGGAAGTTCTg                                                 <  1:1104062/71‑1 (MQ=255)
                     aGGACGGCAAAGGTCCCGAGAAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                 >  1:3748501/1‑71 (MQ=255)
                        aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                              >  1:612465/1‑71 (MQ=255)
                          ggCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGa                                            <  1:3631523/71‑1 (MQ=255)
                            cAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                          <  1:5176462/71‑1 (MQ=255)
                              aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                        >  1:5237249/1‑71 (MQ=11)
                              aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                        >  1:2637157/1‑71 (MQ=11)
                              aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                        >  1:5215938/1‑71 (MQ=11)
                               aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                      <  1:2961535/56‑1 (MQ=255)
                               aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                      <  1:4473877/56‑1 (MQ=255)
                               aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                      <  1:4773430/56‑1 (MQ=255)
                               aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                      <  1:2566969/56‑1 (MQ=255)
                               aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                      <  1:2496179/56‑1 (MQ=255)
                                ggTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGa                                       >  1:1236659/1‑70 (MQ=255)
                                    ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                       >  1:4441074/1‑50 (MQ=255)
                                    ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                                         <  1:122985/64‑1 (MQ=255)
                                    ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACa                                   <  1:3948426/70‑1 (MQ=255)
                                     cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAcatacata                                 <  1:128063/71‑1 (MQ=255)
                                     cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                                 <  1:2566942/71‑1 (MQ=255)
                                       aacAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                               >  1:4922309/3‑71 (MQ=2)
                                       aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                           >  1:2609659/1‑59 (MQ=255)
                                       aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                           >  1:2780819/1‑59 (MQ=255)
                                       aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                           >  1:2727410/1‑59 (MQ=255)
                                       aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                               <  1:2970212/71‑1 (MQ=255)
                                       aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                               <  1:966593/71‑1 (MQ=255)
                                        gCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAcc                              <  1:1008852/71‑1 (MQ=255)
                                        acAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                               >  1:1661494/2‑70 (MQ=2)
                                        acAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                               >  1:1975921/2‑70 (MQ=2)
                                          aacaccacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAg                            >  1:4596896/4‑71 (MQ=2)
                                          aaTACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                  <  1:511256/65‑1 (MQ=255)
                                             accacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAgg                           <  1:1106224/69‑1 (MQ=9)
                                             accaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                       >  1:822193/1‑41 (MQ=255)
                                             accaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                       >  1:2786853/1‑41 (MQ=255)
                                             accaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                       >  1:295208/1‑41 (MQ=255)
                                              ccaccGGTACGGCAATGGTCGGATTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                         >  1:4568543/1‑70 (MQ=255)
                                              ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                       >  1:1083101/1‑40 (MQ=255)
                                              ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                       >  1:1165736/1‑40 (MQ=255)
                                              ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                  >  1:4393492/1‑61 (MQ=255)
                                              ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                         >  1:4377344/1‑70 (MQ=255)
                                              ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAac                        <  1:403968/71‑1 (MQ=255)
                                              ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGAAGGAACAg                                        >  1:2907685/1‑55 (MQ=255)
                                               caccGGTACGGCAATGGTCGGAATCAACGTCGCGCgg                                                         >  1:5433506/1‑37 (MQ=255)
                                               caccGGTACGGCAATGGTCGGAATCAACGTCGCGCgg                                                         >  1:4014668/1‑37 (MQ=255)
                                               caccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                          <  1:5398554/52‑1 (MQ=255)
                                               caccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                          <  1:3971975/52‑1 (MQ=255)
                                                acgggTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAg                            >  1:3124164/4‑65 (MQ=2)
                                                acgggTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAacac                      <  1:5493422/68‑1 (MQ=9)
                                                accGGTACGGCAATGGTCGGAATCAACTTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                      <  1:4566776/71‑1 (MQ=255)
                                                accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                           <  1:4882378/66‑1 (MQ=255)
                                                accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                      <  1:3584368/71‑1 (MQ=255)
                                                accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                      <  1:297754/71‑1 (MQ=255)
                                                accGGCACCGCGATCGTCGGGATCAGCGTCGCCCGCCAATTCTGCAGGAACAGGAACATCACCAGGAagac                      <  1:4785712/71‑1 (MQ=2)
                                                 ccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAg                                        >  1:266931/1‑52 (MQ=255)
                                                 ccGGTACGGCAATGGTAGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                                         >  1:4796944/1‑51 (MQ=255)
                                                   ggTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                   >  1:1448485/1‑71 (MQ=255)
                                                      aCGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgat                >  1:4248622/1‑71 (MQ=255)
                                                       ccgCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGcat                <  1:3775086/68‑1 (MQ=9)
                                                       ccgCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGcat                <  1:786931/68‑1 (MQ=9)
                                                       cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg               <  1:109107/71‑1 (MQ=255)
                                                         gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                          >  1:1191143/1‑42 (MQ=255)
                                                         gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatgat             <  1:1848565/71‑1 (MQ=255)
                                                         gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatgat             <  1:5089969/71‑1 (MQ=255)
                                                          cAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAcc                              <  1:1768456/53‑1 (MQ=2)
                                                           aaTGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGACCACCAGcataatg            >  1:5560352/1‑69 (MQ=2)
                                                           aaTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                >  1:2928644/1‑50 (MQ=255)
                                                            aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGa                                       >  1:1807280/1‑42 (MQ=255)
                                                            aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                  <  1:1821335/47‑1 (MQ=255)
                                                            aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                  <  1:1387417/47‑1 (MQ=255)
                                                            aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                                 <  1:2646836/48‑1 (MQ=255)
                                                            aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                     <  1:2180989/60‑1 (MQ=255)
                                                             tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                           >  1:1701577/1‑37 (MQ=255)
                                                             tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                           <  1:1423131/53‑1 (MQ=255)
                                                             tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                   <  1:1524122/61‑1 (MQ=255)
                                                             tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                   <  1:5267957/61‑1 (MQ=255)
                                                             tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                   <  1:3730336/61‑1 (MQ=255)
                                                             tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                   <  1:3283076/61‑1 (MQ=255)
                                                             tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                   <  1:3154050/61‑1 (MQ=255)
                                                             tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                   <  1:2015780/61‑1 (MQ=255)
                                                              ggTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                   >  1:4168476/1‑60 (MQ=255)
                                                              ggTCGGAATCAACGTCGCGAGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCtt        >  1:5645169/1‑71 (MQ=255)
                                                               gTCGGATTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                   <  1:4376172/59‑1 (MQ=255)
                                                                 cGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa     >  1:1014996/1‑71 (MQ=255)
                                                                    aaTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa  <  1:5031656/71‑1 (MQ=255)
                                                                    |                                                                      
ATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCA  >  minE/311704‑311842

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: