Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,582 A→C 5.6% G97G (GGT→GGG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5820AC 5.6% 347.5 / 15.0 124G97G (GGT→GGGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (65/52);  new base C (3/4);  total (68/56)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 128 total aligned reads displayed.

CTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGT  >  minE/312512‑312652
                                                                      |                                                                      
cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                                        >  1:92574/1‑71 (MQ=255)
   gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                       >  1:2360078/1‑69 (MQ=255)
   gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                       >  1:1210264/1‑69 (MQ=255)
   gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                       >  1:3028572/1‑69 (MQ=255)
    cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                                                     <  1:1260947/70‑1 (MQ=255)
      gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                                  <  1:304415/71‑1 (MQ=255)
         gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                                >  1:2110872/1‑69 (MQ=12)
         gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                                >  1:3535914/1‑69 (MQ=12)
         gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                                                >  1:3531102/1‑69 (MQ=12)
         gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                               <  1:2788957/71‑1 (MQ=255)
         gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                               <  1:3570354/71‑1 (MQ=255)
          caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                                              <  1:3336376/71‑1 (MQ=255)
            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)
                cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                          >  1:2012322/1‑69 (MQ=255)
                 aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                       <  1:1578154/71‑1 (MQ=255)
                 aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                       <  1:841145/71‑1 (MQ=255)
                 aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                                                       <  1:2178774/71‑1 (MQ=255)
                    gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg                                                     >  1:1257401/1‑70 (MQ=255)
                    gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                                                      <  1:3389683/69‑1 (MQ=255)
                    gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                                                    >  1:3445840/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:3038001/71‑1 (MQ=255)
                     ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                                   <  1:1073045/71‑1 (MQ=255)
                      cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                                                  >  1:3511589/1‑71 (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:1887770/51‑1 (MQ=12)
                         cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                                                   <  1:734976/51‑1 (MQ=12)
                         cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                                                   <  1:1655155/51‑1 (MQ=12)
                         cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                                 >  1:3189015/1‑69 (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:2848358/69‑1 (MQ=255)
                            cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                              <  1:3802011/69‑1 (MQ=255)
                            cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                              <  1:2009352/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:1632188/54‑1 (MQ=255)
                               tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                          <  1:2936538/54‑1 (MQ=255)
                               tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                          <  1:633519/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)
                                 tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                                 >  1:674896/1‑61 (MQ=255)
                                 tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                                 >  1:3636047/1‑61 (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:1779245/1‑71 (MQ=255)
                                           cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                             >  1:1760959/1‑71 (MQ=255)
                                            tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGGACCCGTGGAGt                             >  1:765945/1‑70 (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)
                                                  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)
                                                        atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTag                    <  1:2601741/67‑1 (MQ=255)
                                                          atCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                                   >  1:3374315/1‑50 (MQ=255)
                                                           tCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                                 >  1:2764503/1‑35 (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:3645555/1‑70 (MQ=255)
                                                             cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg            >  1:598377/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)
                                                                      aCCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACt                         <  1:2519534/48‑1 (MQ=255)
                                                                      aCCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGt  >  1:1452009/1‑71 (MQ=255)
                                                                      |                                                                      
CTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGT  >  minE/312512‑312652

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: