Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 313,525 G→A 6.0% N188N (AAC→AAT) ‡ acrA ← multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE313,5250GA 6.0% 302.2 / 5.5 100N188N (AAC→AAT) ‡acrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (54/40);  new base A (2/4);  total (56/44)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGT  >  minE/313456‑313593
                                                                     |                                                                    
cGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCTCCTGACCGTTCTGTACCAATGCGCCATCCGTCACGt                                                                     >  1:1003895/1‑71 (MQ=255)
 ggATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGt                                                                     <  1:4225191/70‑1 (MQ=255)
  gATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTc                                                                   <  1:1810057/71‑1 (MQ=255)
  gATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTc                                                                   <  1:1680452/71‑1 (MQ=255)
   aTCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTc                                                                   <  1:5359403/70‑1 (MQ=255)
    tCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGa                                                                 <  1:1690203/71‑1 (MQ=255)
    tCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGa                                                                 >  1:3070621/1‑71 (MQ=255)
    tCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGa                                                                 >  1:4460138/1‑71 (MQ=255)
     cAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGa                                                                 >  1:5087788/1‑70 (MQ=255)
      aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGAc                                                                >  1:4574650/1‑70 (MQ=255)
      aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGAc                                                                >  1:2868680/1‑70 (MQ=255)
      aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                               >  1:2749826/1‑71 (MQ=255)
      aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                               >  1:4640144/1‑71 (MQ=255)
      aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                               >  1:1856073/1‑71 (MQ=255)
      aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                               >  1:3475514/1‑71 (MQ=255)
      aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                               >  1:1450698/1‑71 (MQ=255)
      aaGTTGCTGCACGGGTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                               <  1:1732072/71‑1 (MQ=255)
           gctgcACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTAcc                                                           <  1:540531/70‑1 (MQ=255)
           gctgcACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTAcc                                                           <  1:3693696/70‑1 (MQ=255)
             tgcACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCg                                                                  <  1:3565365/61‑1 (MQ=255)
             tgcACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCg                                                                  <  1:310165/61‑1 (MQ=255)
              gcACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAAt                                                        >  1:4312313/1‑70 (MQ=255)
               cACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGc                                                      >  1:5399426/1‑71 (MQ=255)
                aCGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGc                                                      >  1:1981536/1‑70 (MQ=255)
                aCGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCg                                                     <  1:1841616/71‑1 (MQ=255)
                 cGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTc                                                                   <  1:5621440/56‑1 (MQ=255)
                  ggTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGAc                                                   >  1:1172146/1‑71 (MQ=255)
                   gTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGc                                                      <  1:3544861/67‑1 (MQ=255)
                   gTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGc                                                      <  1:1699255/67‑1 (MQ=255)
                     tGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGa                                                    >  1:2442808/1‑67 (MQ=255)
                     tGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGa                                                    >  1:2772511/1‑67 (MQ=255)
                     tGCCAGCGCAGTCGCCTGACCGCTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGa                                                    >  1:3024802/1‑67 (MQ=255)
                       ccAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCaa                                                         >  1:3156769/1‑60 (MQ=255)
                       ccAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTa                                              <  1:5544201/71‑1 (MQ=255)
                         aGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAAt                                            >  1:3808787/1‑71 (MQ=255)
                          gcgcAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATc                                           >  1:5354033/1‑71 (MQ=255)
                            gcAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATc                                           >  1:3150097/1‑69 (MQ=255)
                            gcAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCgg                                         <  1:2841279/71‑1 (MQ=255)
                             cAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGga                                        <  1:2631577/71‑1 (MQ=255)
                             cAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGA‑CGCTAATCGga                                        <  1:4390026/70‑1 (MQ=255)
                               gTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGga                                        >  1:1120198/1‑69 (MQ=255)
                                  gCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAgg                                    <  1:4347586/70‑1 (MQ=255)
                                       aCCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCt                                               <  1:5126496/54‑1 (MQ=255)
                                       aCCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGa                                 >  1:1387643/1‑68 (MQ=255)
                                       aCCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACtt                              >  1:1000774/1‑71 (MQ=255)
                                       aCCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACt                               <  1:4078213/70‑1 (MQ=255)
                                        ccGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCgg                                         <  1:2773031/59‑1 (MQ=255)
                                        ccGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCgg                                         <  1:2058600/59‑1 (MQ=255)
                                         cGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGc                                                      >  1:2714395/1‑45 (MQ=255)
                                         cGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTaa                                             >  1:1363686/1‑54 (MQ=255)
                                          gTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGgag                                       >  1:2922155/1‑59 (MQ=255)
                                           ttCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACttt                             <  1:3908390/68‑1 (MQ=255)
                                           ttCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACttt                             <  1:1507394/68‑1 (MQ=255)
                                           ttCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgt                          >  1:1044453/1‑71 (MQ=255)
                                           ttCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgt                          >  1:4168834/1‑71 (MQ=255)
                                           ttCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgt                          >  1:1079468/1‑71 (MQ=255)
                                            tCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtg                         >  1:284897/1‑71 (MQ=255)
                                            tCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtg                         >  1:278007/1‑71 (MQ=255)
                                             cTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtg                         >  1:3060659/1‑70 (MQ=255)
                                                tACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtgt                        <  1:3958640/68‑1 (MQ=255)
                                                tACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTa                       <  1:1234323/69‑1 (MQ=255)
                                                tACCAATGCGCCTTCCGTCACGTGCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtgt                        <  1:2560667/68‑1 (MQ=255)
                                                 aCCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGc                    >  1:4053291/1‑71 (MQ=255)
                                                    aaTGCGGCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAg                 <  1:4957902/71‑1 (MQ=255)
                                                      tGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtg                         >  1:1280405/1‑61 (MQ=255)
                                                      tGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtg                         >  1:2424949/1‑61 (MQ=255)
                                                       gcgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGga                                        >  1:877363/1‑45 (MQ=255)
                                                        cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGga                                        >  1:885729/1‑44 (MQ=255)
                                                        cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGga                                        >  1:609626/1‑44 (MQ=255)
                                                        cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGAtt              >  1:1815974/1‑70 (MQ=255)
                                                        cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGAtt              <  1:555686/70‑1 (MQ=255)
                                                        cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGAtt              >  1:5604312/1‑70 (MQ=255)
                                                        cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTg             <  1:3291585/71‑1 (MQ=255)
                                                         gcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTg             <  1:345774/70‑1 (MQ=255)
                                                           cTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATc          >  1:5429045/1‑71 (MQ=255)
                                                            ttCGGTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGc                    <  1:410310/60‑1 (MQ=255)
                                                            ttCGGTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGc                    <  1:5609237/60‑1 (MQ=255)
                                                            ttCGGTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGc                    <  1:4492323/60‑1 (MQ=255)
                                                            ttCGGTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGCAAGATTGAt           <  1:2698399/69‑1 (MQ=255)
                                                            ttCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAg                 <  1:333139/63‑1 (MQ=255)
                                                            ttCCGTCACGTTCGACTTACCAATACGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATc          >  1:1768711/1‑70 (MQ=255)
                                                             tcggTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGCAAGATTGa            >  1:2488414/4‑67 (MQ=255)
                                                             tcggTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGCAAGATTGa            >  1:1674694/4‑67 (MQ=255)
                                                              ccGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgc       >  1:5309105/1‑71 (MQ=255)
                                                              ccGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgc       >  1:747901/1‑71 (MQ=255)
                                                              ccGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcg        >  1:961809/1‑70 (MQ=255)
                                                               cGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgcg      <  1:2819622/71‑1 (MQ=255)
                                                               cGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgcg      <  1:2231851/71‑1 (MQ=255)
                                                               cGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgcg      >  1:3528495/1‑71 (MQ=255)
                                                                gTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTg                            <  1:4627500/48‑1 (MQ=255)
                                                                gTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTg                            <  1:1119593/48‑1 (MQ=255)
                                                                 tCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtgt                        <  1:5474630/51‑1 (MQ=255)
                                                                 tCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCa    <  1:2370381/71‑1 (MQ=255)
                                                                  cACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGa                >  1:105230/1‑58 (MQ=255)
                                                                  cACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAg   >  1:1840239/1‑71 (MQ=255)
                                                                   aCGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTg             >  1:3201007/1‑60 (MQ=255)
                                                                   aCGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGAt           >  1:5252971/1‑62 (MQ=255)
                                                                   aCGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGt  >  1:3249714/1‑71 (MQ=255)
                                                                   aCGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACATTGGTGTAAGcc                   >  1:378175/1‑54 (MQ=255)
                                                                    cGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACtt                              <  1:5043541/42‑1 (MQ=255)
                                                                     |                                                                    
CGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGT  >  minE/313456‑313593

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: ‑