Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 313,863 T→C 17.6% I76V (ATC→GTC) ‡ acrA ← multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE313,8630TC17.6% 194.3 / 33.1 91I76V (ATC→GTC) ‡acrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (42/33);  new base C (10/6);  total (52/39)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTG  >  minE/313793‑313933
                                                                      |                                                                      
gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt                                                                        <  1:772659/71‑1 (MQ=255)
gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt                                                                        <  1:261144/71‑1 (MQ=255)
  aGGATCAATCTGAAAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGTTTCAGGATaa                                                                      <  1:2256437/71‑1 (MQ=255)
      tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATccc                                                                  >  1:3630095/1‑71 (MQ=255)
      tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATcc                                                                   >  1:88970/1‑70 (MQ=255)
            tGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAAc                                                            >  1:5186021/1‑71 (MQ=255)
              aTAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTa                                                              <  1:472724/67‑1 (MQ=255)
              aTAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACtt                                                          <  1:4302614/71‑1 (MQ=255)
              aTAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACt                                                           <  1:4365912/70‑1 (MQ=255)
                 gagagagaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAg                                                       <  1:3355624/71‑1 (MQ=255)
                  agagagaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAgg                                                      <  1:2898845/71‑1 (MQ=255)
                      agaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAAt                                                                     >  1:4562670/1‑52 (MQ=255)
                      agaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGc                                                                >  1:3336319/1‑57 (MQ=255)
                        acacCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAAc                                                <  1:5298/71‑1 (MQ=255)
                           ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt                                                                        <  1:688721/44‑1 (MQ=255)
                           ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTg                                                         <  1:755252/59‑1 (MQ=255)
                           ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTg                                                         <  1:320484/59‑1 (MQ=255)
                           ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACg                                                   <  1:3472635/65‑1 (MQ=255)
                             tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                                      >  1:2894055/1‑44 (MQ=255)
                             tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGa                                                     >  1:1226219/1‑61 (MQ=255)
                             tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTc                                             <  1:390887/69‑1 (MQ=255)
                             tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTc                                             <  1:268844/69‑1 (MQ=255)
                              gCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGaa                                                 >  1:3313676/1‑64 (MQ=255)
                              gCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGaa                                                 >  1:5434545/1‑64 (MQ=255)
                               cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                                      >  1:403855/1‑42 (MQ=255)
                               cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                                      >  1:5225794/1‑42 (MQ=255)
                               cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTg                                           >  1:4422835/1‑69 (MQ=255)
                               cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTg                                           >  1:3209556/1‑69 (MQ=255)
                               cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTg                                           >  1:1906437/1‑69 (MQ=255)
                                  cGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATc                                      >  1:4141687/1‑71 (MQ=255)
                                    acatCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACt                                               <  1:4917131/57‑1 (MQ=255)
                                    acatCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACt                                               <  1:5261606/57‑1 (MQ=255)
                                    aTGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAg                                                       <  1:3123355/52‑1 (MQ=255)
                                     tGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCgg                                   <  1:4271291/71‑1 (MQ=255)
                                      gTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                                      >  1:1187704/1‑35 (MQ=255)
                                      gTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGc                                                                >  1:4953699/1‑41 (MQ=255)
                                      gTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGc                                                                >  1:4344865/1‑41 (MQ=255)
                                      gTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCgg                                   >  1:4521836/1‑70 (MQ=255)
                                      atCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGc                                          <  1:3393494/62‑1 (MQ=255)
                                       tCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACgat                                  >  1:2845283/1‑68 (MQ=255)
                                         gCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACtt                                              <  1:2650266/56‑1 (MQ=255)
                                         gCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAgg                               >  1:4051084/1‑71 (MQ=255)
                                          cTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGc                              <  1:254615/71‑1 (MQ=255)
                                            aCCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGt                                  <  1:1228423/65‑1 (MQ=255)
                                            aCCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGt                                  <  1:5230191/65‑1 (MQ=255)
                                             ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCt                                                               <  1:991978/35‑1 (MQ=255)
                                             ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAAc                                                <  1:2626268/50‑1 (MQ=255)
                                             ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAAc                                                <  1:5564794/50‑1 (MQ=255)
                                             ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACt                           >  1:1876190/1‑71 (MQ=255)
                                             ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACt                           >  1:1869642/1‑71 (MQ=255)
                                             ccTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCatt                           <  1:2516331/71‑3 (MQ=255)
                                               ttcttTGAAATTACGCTTCAGGATAATCCCGCTaa                                                             >  1:3772504/1‑35 (MQ=255)
                                               ttcttTGAAATTACGCTTCAGGATAATCCCGCTAACTTg                                                         <  1:5387738/39‑1 (MQ=255)
                                                tcttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACt                                               >  1:1237771/1‑48 (MQ=255)
                                                 cttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACtt                                              >  1:4605907/1‑48 (MQ=255)
                                                 cttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACtt                                              >  1:3856700/1‑48 (MQ=255)
                                                 cttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTg                       >  1:3984667/1‑71 (MQ=255)
                                                  tttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGc                      >  1:513687/1‑71 (MQ=255)
                                                  tttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGc                      >  1:2760457/1‑71 (MQ=255)
                                                   ttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATcc                                     <  1:2144178/55‑1 (MQ=255)
                                                   gtGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCTGCTATACGATAAGCATTGGTGCg                     <  1:2343094/70‑1 (MQ=255)
                                                    tGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCgg                    >  1:1375392/1‑70 (MQ=255)
                                                     gAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGa                    <  1:5551224/70‑1 (MQ=255)
                                                      aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGt                        >  1:1981234/1‑65 (MQ=255)
                                                      aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg                       >  1:1572776/1‑66 (MQ=255)
                                                      aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg                       >  1:536100/1‑66 (MQ=255)
                                                      aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg                       >  1:3001702/1‑66 (MQ=255)
                                                      aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCggc                   >  1:4205515/1‑68 (MQ=255)
                                                      aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCGGcc                  <  1:4455268/71‑1 (MQ=255)
                                                      aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCGGcc                  >  1:2215143/1‑71 (MQ=255)
                                                      aaaTTACTCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGAc                   <  1:5693475/70‑1 (MQ=255)
                                                      aaaTTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGAcc                  >  1:2856905/1‑71 (MQ=255)
                                                        aTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg                          >  1:5356662/1‑61 (MQ=255)
                                                        aTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg                          >  1:5135467/1‑61 (MQ=255)
                                                        aTTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCg                <  1:2536126/71‑1 (MQ=255)
                                                         ttACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCgg               <  1:4605918/71‑1 (MQ=255)
                                                         ttACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCgg               <  1:2171877/71‑1 (MQ=255)
                                                           aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCg                                    >  1:4701667/1‑48 (MQ=255)
                                                           aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCg                                    >  1:4708280/1‑48 (MQ=255)
                                                           aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCg                                    >  1:1544449/1‑48 (MQ=255)
                                                           aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGa              <  1:2308429/70‑1 (MQ=255)
                                                             gCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGc           >  1:139169/1‑71 (MQ=255)
                                                              cTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCt          <  1:1138572/71‑1 (MQ=255)
                                                                tCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCgg                                   <  1:2520358/44‑1 (MQ=255)
                                                                tCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTc         >  1:4285273/1‑70 (MQ=255)
                                                                 cAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCgg       >  1:3156673/1‑71 (MQ=255)
                                                                 cAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCgg       >  1:723581/1‑71 (MQ=255)
                                                                 cAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCg        <  1:5147570/70‑1 (MQ=255)
                                                                  aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTa                                 >  1:1272766/1‑44 (MQ=255)
                                                                  aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTg                          >  1:2580969/1‑51 (MQ=255)
                                                                  aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTg                          >  1:887874/1‑51 (MQ=255)
                                                                  aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTg                          >  1:1081018/1‑51 (MQ=255)
                                                                   ggATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGtt     >  1:1133632/1‑71 (MQ=255)
                                                                    gATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtg    <  1:663073/71‑1 (MQ=255)
                                                                      tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCa                             >  1:5552838/1‑44 (MQ=255)
                                                                      tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCa                             >  1:480781/1‑44 (MQ=255)
                                                                      tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg  <  1:157456/71‑1 (MQ=255)
                                                                      tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg  <  1:4781529/71‑1 (MQ=255)
                                                                      tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg  <  1:4601654/71‑1 (MQ=255)
                                                                      tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg  <  1:2803890/71‑1 (MQ=255)
                                                                      |                                                                      
GCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTG  >  minE/313793‑313933

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: