Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,579 A→G 6.4% T98T (ACT→ACC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5790AG 6.4% 358.6 / 12.9 124T98T (ACT→ACCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (64/52);  new base G (3/5);  total (67/57)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 128 total aligned reads displayed.

TCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACAT  >  minE/312511‑312647
                                                                    |                                                                    
tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                                     >  1:35390/1‑70 (MQ=255)
 cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                                   >  1:92574/1‑71 (MQ=255)
    gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                  >  1:1210264/1‑69 (MQ=255)
    gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                  >  1:3028572/1‑69 (MQ=255)
    gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                  >  1:2360078/1‑69 (MQ=255)
     cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                                                <  1:1260947/70‑1 (MQ=255)
       gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                             <  1:304415/71‑1 (MQ=255)
          gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                           >  1:3535914/1‑69 (MQ=12)
          gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                           >  1:3531102/1‑69 (MQ=12)
          gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                           >  1:2110872/1‑69 (MQ=12)
          gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                          <  1:3570354/71‑1 (MQ=255)
          gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                          <  1:2788957/71‑1 (MQ=255)
             aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                       >  1:2376829/1‑71 (MQ=255)
             aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                       >  1:2351372/1‑71 (MQ=255)
               gggaTCGACAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                     <  1:3347239/68‑1 (MQ=255)
               ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGa                                                              <  1:2990753/62‑1 (MQ=255)
                 cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                     >  1:2012322/1‑69 (MQ=255)
                 cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                     >  1:3728487/1‑69 (MQ=255)
                  aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                  <  1:2178774/71‑1 (MQ=255)
                  aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                  <  1:1578154/71‑1 (MQ=255)
                  aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                  <  1:841145/71‑1 (MQ=255)
                     gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg                                                >  1:1257401/1‑70 (MQ=255)
                     gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                                                 <  1:3389683/69‑1 (MQ=255)
                     gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                               >  1:3393173/1‑71 (MQ=255)
                     gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                               >  1:3519374/1‑71 (MQ=255)
                     gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                               >  1:1535932/1‑71 (MQ=255)
                     gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                               >  1:3957681/1‑71 (MQ=255)
                     gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCACGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                               >  1:2122566/1‑71 (MQ=255)
                      ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                              <  1:1073045/71‑1 (MQ=255)
                       cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                                             >  1:2646272/1‑71 (MQ=255)
                       aaGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                                             >  1:3665785/2‑71 (MQ=255)
                        aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                            >  1:2523077/1‑71 (MQ=255)
                        aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGAATCAGTACCAGACTCAGAGGTCAGGGTGAt                                            >  1:3466175/1‑71 (MQ=255)
                        aGCTGCAGTTTGTTATGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                            >  1:3491305/1‑71 (MQ=255)
                          cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                                              <  1:734976/51‑1 (MQ=12)
                          cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                                              <  1:3258854/51‑1 (MQ=12)
                          cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                                              <  1:1887770/51‑1 (MQ=12)
                          cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                                              <  1:1655155/51‑1 (MQ=12)
                          cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                            >  1:3189015/1‑69 (MQ=255)
                          cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGAATCAGTACCAGACTCAAAGGTCAGGGTGATCt                                          >  1:2820966/1‑71 (MQ=255)
                          cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCAGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                            >  1:318973/1‑69 (MQ=255)
                           tGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                         >  1:3588475/1‑71 (MQ=255)
                             cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                  <  1:3910271/60‑1 (MQ=255)
                             cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                         <  1:2009352/69‑1 (MQ=255)
                             cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                         <  1:3802011/69‑1 (MQ=255)
                               gTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACg                                     >  1:1432303/1‑71 (MQ=255)
                               gTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACg                                     >  1:3967590/1‑71 (MQ=255)
                                tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                     <  1:2936538/54‑1 (MQ=255)
                                tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                     <  1:1632188/54‑1 (MQ=255)
                                tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                     <  1:1354870/54‑1 (MQ=255)
                                tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                     <  1:633519/54‑1 (MQ=255)
                                tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                              <  1:978370/61‑1 (MQ=255)
                                tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACgg                                    <  1:3506615/71‑1 (MQ=255)
                                tttGTTCTGCACCTGAACCTGCGCGATATCAGCATCAGTACCAGACTCAAAGGTCAGGGTg                                              <  1:3755505/61‑1 (MQ=255)
                                  tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                       >  1:509066/1‑50 (MQ=255)
                                  tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                            >  1:3636047/1‑61 (MQ=255)
                                  tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                            >  1:674896/1‑61 (MQ=255)
                                  tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                                  >  1:2535739/1‑71 (MQ=255)
                                  tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                                  >  1:1927374/1‑71 (MQ=255)
                                      cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                          >  1:507193/1‑43 (MQ=255)
                                      cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                              >  1:945818/1‑71 (MQ=255)
                                       tGCACCTGAAGCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGt                             >  1:1220312/1‑71 (MQ=255)
                                        gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                            >  1:2641954/1‑71 (MQ=255)
                                         cACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc                               >  1:4083189/1‑67 (MQ=255)
                                            cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAg                         >  1:1042034/1‑70 (MQ=255)
                                            cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                        >  1:1760959/1‑71 (MQ=255)
                                            cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                        >  1:1779245/1‑71 (MQ=255)
                                             tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGGACCCGTGGAGt                        >  1:765945/1‑70 (MQ=255)
                                              gAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                                         >  1:138925/1‑36 (MQ=255)
                                              gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa                      >  1:3037525/1‑71 (MQ=255)
                                              gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa                      >  1:2510550/1‑71 (MQ=255)
                                               aaCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                            <  1:922801/64‑1 (MQ=255)
                                               aaCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc                       <  1:4016983/69‑1 (MQ=255)
                                                  cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                  <  1:2621067/71‑1 (MQ=255)
                                                   tGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCAc                     <  1:2003998/67‑1 (MQ=255)
                                                    gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                                       >  1:1874241/1‑48 (MQ=255)
                                                    gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                                       >  1:4011056/1‑48 (MQ=255)
                                                    gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                  >  1:813574/1‑69 (MQ=255)
                                                     cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGATTCACTGTTag               <  1:2926808/71‑1 (MQ=255)
                                                     cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTa                <  1:1194559/70‑1 (MQ=255)
                                                     cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGGGATCTGCACGGTACCCGTGGAGTCACTGTTag               <  1:2374867/71‑1 (MQ=255)
                                                       cgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc                       <  1:1127175/61‑1 (MQ=255)
                                                           atCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                              >  1:3374315/1‑50 (MQ=255)
                                                           atCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTag               <  1:299290/65‑1 (MQ=255)
                                                            tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                            >  1:90525/1‑35 (MQ=255)
                                                            tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACa         <  1:2803921/70‑1 (MQ=255)
                                                              cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc                               >  1:3930323/1‑46 (MQ=255)
                                                              cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg       >  1:598377/1‑70 (MQ=255)
                                                              cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg       >  1:3645555/1‑70 (MQ=255)
                                                              cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg       >  1:2737985/1‑70 (MQ=255)
                                                              cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg       <  1:894640/70‑1 (MQ=255)
                                                               gCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa     <  1:1712299/71‑1 (MQ=255)
                                                                cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc    <  1:4011723/71‑1 (MQ=255)
                                                                 aTCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACa         >  1:3569084/1‑65 (MQ=255)
                                                                 aTCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc    >  1:3707430/1‑70 (MQ=255)
                                                                  tCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACAt  <  1:1892127/71‑1 (MQ=255)
                                                                   cAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                  <  1:3005516/54‑1 (MQ=255)
                                                                   cAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACAt  >  1:2086060/1‑70 (MQ=255)
                                                                    aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa     >  1:951507/1‑66 (MQ=255)
                                                                    aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACAt  >  1:66653/1‑69 (MQ=255)
                                                                    |                                                                    
TCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACAT  >  minE/312511‑312647

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: