Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,585 A→G 5.3% S96S (TCT→TCC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5850AG 5.3% 390.7 / 10.0 132S96S (TCT→TCCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (67/58);  new base G (3/4);  total (70/62)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 136 total aligned reads displayed.

CGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAT  >  minE/312516‑312655
                                                                     |                                                                      
cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                                                        <  1:1260947/70‑1 (MQ=255)
  gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                                     <  1:304415/71‑1 (MQ=255)
     gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                                   >  1:3535914/1‑69 (MQ=12)
        aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                               >  1:2351372/1‑71 (MQ=255)
        aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                               >  1:2376829/1‑71 (MQ=255)
          gggaTCGACAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                             <  1:3347239/68‑1 (MQ=255)
          ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGa                                                                      <  1:2990753/62‑1 (MQ=255)
            cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                             >  1:3728487/1‑69 (MQ=255)
             aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                          <  1:841145/71‑1 (MQ=255)
             aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                          <  1:2178774/71‑1 (MQ=255)
             aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                          <  1:1578154/71‑1 (MQ=255)
                gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg                                                        >  1:1257401/1‑70 (MQ=255)
                gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                                                         <  1:3389683/69‑1 (MQ=255)
                gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                                       >  1:3957681/1‑71 (MQ=255)
                gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                                       >  1:1535932/1‑71 (MQ=255)
                gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                                       >  1:3445840/1‑71 (MQ=255)
                gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                                       >  1:3519374/1‑71 (MQ=255)
                gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCACGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                                       >  1:2122566/1‑71 (MQ=255)
                 ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                                      <  1:3038001/71‑1 (MQ=255)
                 ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                                      <  1:1073045/71‑1 (MQ=255)
                  cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                                                     >  1:2646272/1‑71 (MQ=255)
                  cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                                                     >  1:3511589/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:1887770/51‑1 (MQ=12)
                     cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                                                      <  1:1655155/51‑1 (MQ=12)
                     cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                                                      <  1:3258854/51‑1 (MQ=12)
                     cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                                      <  1:1203974/67‑1 (MQ=255)
                     cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                                    >  1:3189015/1‑69 (MQ=255)
                     cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGAATCAGTACCAGACTCAAAGGTCAGGGTGATCt                                                  >  1:2820966/1‑71 (MQ=255)
                      tGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                                 >  1:3588475/1‑71 (MQ=255)
                        cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                          <  1:3910271/60‑1 (MQ=255)
                        cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                                 <  1:3822765/69‑1 (MQ=255)
                        cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                                 <  1:2009352/69‑1 (MQ=255)
                        cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                                 <  1:2848358/69‑1 (MQ=255)
                          gTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACg                                             >  1:1432303/1‑71 (MQ=255)
                          gTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACg                                             >  1:3967590/1‑71 (MQ=255)
                           tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                                              <  1:2998188/53‑1 (MQ=255)
                           tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                             <  1:633519/54‑1 (MQ=255)
                           tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                             <  1:2936538/54‑1 (MQ=255)
                           tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                                      <  1:978370/61‑1 (MQ=255)
                           tttGTTCTGCACCTGAACCTGCGCGATATCAGCATCAGTACCAGACTCAAAGGTCAGGGTg                                                      <  1:3755505/61‑1 (MQ=255)
                             tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                               >  1:509066/1‑50 (MQ=255)
                             tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCa                                                           >  1:3556662/1‑54 (MQ=255)
                             tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                                    >  1:674896/1‑61 (MQ=255)
                             tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                                               <  1:625446/66‑1 (MQ=255)
                             tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                                          >  1:1927374/1‑71 (MQ=255)
                               ttCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAcc                                        >  1:3602160/1‑71 (MQ=255)
                                 cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                                  >  1:507193/1‑43 (MQ=255)
                                       cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAg                                 >  1:1042034/1‑70 (MQ=255)
                                       cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                                >  1:1760959/1‑71 (MQ=255)
                                       cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                                >  1:1779245/1‑71 (MQ=255)
                                         gAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                                                 >  1:138925/1‑36 (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)
                                             cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                          >  1:1798092/1‑71 (MQ=255)
                                              tGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCAc                             <  1:2003998/67‑1 (MQ=255)
                                               gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                          >  1:813574/1‑69 (MQ=255)
                                                cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGATTCACTGTTag                       <  1:2926808/71‑1 (MQ=255)
                                                cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTa                        <  1:1194559/70‑1 (MQ=255)
                                                cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTa                        <  1:932816/70‑1 (MQ=255)
                                                cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGGGATCTGCACGGTACCCGTGGAGTCACTGTTag                       <  1:2374867/71‑1 (MQ=255)
                                                  cgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc                               <  1:1127175/61‑1 (MQ=255)
                                                    atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTag                       <  1:2601741/67‑1 (MQ=255)
                                                     tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACgg                                            <  1:2038381/45‑1 (MQ=255)
                                                      atCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                                      >  1:3374315/1‑50 (MQ=255)
                                                       tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                                    >  1:2764503/1‑35 (MQ=255)
                                                       tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                                    >  1:90525/1‑35 (MQ=255)
                                                       tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGa                                  <  1:842902/53‑1 (MQ=255)
                                                       tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACa                 <  1:2803921/70‑1 (MQ=255)
                                                         cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg               >  1:3645555/1‑70 (MQ=255)
                                                         cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg               >  1:2737985/1‑70 (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:3984926/71‑1 (MQ=255)
                                                             tCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACAt          <  1:1892127/71‑1 (MQ=255)
                                                              cAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                          <  1:3005516/54‑1 (MQ=255)
                                                              cAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACAt          >  1:2086060/1‑70 (MQ=255)
                                                               aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc                               >  1:4122432/1‑48 (MQ=255)
                                                               aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTa                        >  1:1674557/1‑55 (MQ=255)
                                                               aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa             >  1:951507/1‑66 (MQ=255)
                                                                  aCCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACt                            <  1:2519534/48‑1 (MQ=255)
                                                                  aCCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGt     >  1:1452009/1‑71 (MQ=255)
                                                                   ccAGACTCAAAGGTCAGGGTGATCTTCACGGTACCCGTGGAGTCACTGTTag                       <  1:703865/52‑1 (MQ=255)
                                                                   ccAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                          >  1:3242456/1‑49 (MQ=255)
                                                                   ccAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                          >  1:3237600/1‑49 (MQ=255)
                                                                   ccAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTag                       <  1:3831042/52‑1 (MQ=255)
                                                                   ccAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTag                       <  1:2066498/52‑1 (MQ=255)
                                                                   ccAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGtt    <  1:1006709/71‑1 (MQ=255)
                                                                    cAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                          <  1:1239490/48‑1 (MQ=255)
                                                                    cAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg               <  1:3607766/59‑1 (MQ=255)
                                                                     aGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                          <  1:407650/47‑1 (MQ=255)
                                                                     aGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt  >  1:217562/1‑71 (MQ=255)
                                                                     |                                                                      
CGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAT  >  minE/312516‑312655

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: