Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,561 C→T 8.2% Q104Q (CAG→CAA acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5610CT 8.2% 360.6 / 14.9 110Q104Q (CAG→CAAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (59/41);  new base T (3/6);  total (63/47)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 112 total aligned reads displayed.

CCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGT  >  minE/312494‑312631
                                                                   |                                                                      
ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                                                     >  1:2161056/1‑71 (MQ=255)
ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                                                     >  1:3584377/1‑71 (MQ=255)
ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                                                     <  1:3909276/71‑1 (MQ=255)
   cTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                                                   >  1:144054/1‑70 (MQ=255)
      gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                                <  1:1439362/70‑1 (MQ=255)
      gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACATGCGCGata                                                               >  1:4121072/1‑71 (MQ=255)
       ctgctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGat                                                                <  1:3652408/69‑1 (MQ=255)
         gctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATcc                                                            >  1:3016428/1‑71 (MQ=255)
          ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                                <  1:882196/66‑1 (MQ=255)
          ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg                                                           >  1:2171676/1‑71 (MQ=255)
            gAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                                                                      <  1:510482/58‑1 (MQ=255)
             aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                        <  1:3180559/71‑1 (MQ=255)
              aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                        >  1:4026884/1‑70 (MQ=255)
               cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                                >  1:2182012/1‑61 (MQ=255)
               cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                                >  1:1338591/1‑61 (MQ=255)
               cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                                >  1:1272280/1‑61 (MQ=255)
                 tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                     >  1:35390/1‑70 (MQ=255)
                  cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                   >  1:92574/1‑71 (MQ=255)
                     gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                  >  1:3028572/1‑69 (MQ=255)
                     gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                  >  1:1210264/1‑69 (MQ=255)
                      cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                                <  1:1260947/70‑1 (MQ=255)
                        gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                             <  1:304415/71‑1 (MQ=255)
                           gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                           >  1:3531102/1‑69 (MQ=12)
                           gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                           >  1:3535914/1‑69 (MQ=12)
                           gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                           >  1:2110872/1‑69 (MQ=12)
                           gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata                                                               <  1:1416572/50‑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)
                               cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTTAACCTGCGCGATATCCGCAt                                                        <  1:791388/53‑1 (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: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:3519374/1‑71 (MQ=255)
                                      gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                               >  1:1535932/1‑71 (MQ=255)
                                      gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                               >  1:3393173/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)
                                        cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat                                                              >  1:1442397/1‑38 (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:3258854/51‑1 (MQ=12)
                                           cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                              <  1:734976/51‑1 (MQ=12)
                                           cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                              <  1:1655155/51‑1 (MQ=12)
                                           cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                              <  1:1887770/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)
                                           cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCAGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                            >  1:318973/1‑69 (MQ=255)
                                            tGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAtg                                                       <  1:3866301/41‑2 (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)
                                              cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                         <  1:3802011/69‑1 (MQ=255)
                                               aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                      <  1:2507194/71‑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)
                                                 tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                     <  1:1354870/54‑1 (MQ=255)
                                                 tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                              <  1:978370/61‑1 (MQ=255)
                                                 tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACgg                    <  1:3506615/71‑1 (MQ=255)
                                                   tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                       >  1:509066/1‑50 (MQ=255)
                                                   tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                            >  1:3636047/1‑61 (MQ=255)
                                                   tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                            >  1:674896/1‑61 (MQ=255)
                                                   tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                       <  1:625446/66‑1 (MQ=255)
                                                   tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                  >  1:1927374/1‑71 (MQ=255)
                                                   tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                  >  1:2535739/1‑71 (MQ=255)
                                                     ttCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAcc                >  1:3602160/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)
                                                             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)
                                                               gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa      >  1:3037525/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)
                                                                   |                                                                      
CCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGT  >  minE/312494‑312631

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: