Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,561 C→T 16.8% 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,5610CT16.8% 250.0 / 33.8 101Q104Q (CAG→CAAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (50/34);  new base T (11/6);  total (61/40)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 106 total aligned reads displayed.

TCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCA  >  minE/312491‑312627
                                                                      |                                                                  
tCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAcc                                                                    >  1:2304822/1‑71 (MQ=255)
  aCCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCt                                                                   <  1:1199429/70‑1 (MQ=255)
  aCCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                                                                  >  1:124613/1‑71 (MQ=255)
  aCCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                                                                  >  1:444578/1‑71 (MQ=255)
    cccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcg                                                                >  1:5373523/1‑71 (MQ=255)
    cccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcg                                                                >  1:3268437/1‑71 (MQ=255)
     ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                                               >  1:1201616/1‑71 (MQ=255)
     ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                                               <  1:1253011/71‑1 (MQ=255)
     ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                                               >  1:3131442/1‑71 (MQ=255)
     ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                                               >  1:1830243/1‑71 (MQ=255)
     ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcg                                                                >  1:630402/1‑70 (MQ=255)
        tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGa                                                             >  1:2052532/1‑70 (MQ=255)
         gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata                                                           >  1:3822746/1‑71 (MQ=255)
         gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata                                                           >  1:5065724/1‑71 (MQ=255)
          ctgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat                                                          >  1:3947693/1‑71 (MQ=255)
             ctgAACTTCTTGCGGCAGCACCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGCACCTGCGCGATATCCg                                                       >  1:2945349/1‑71 (MQ=255)
                 aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                    >  1:448469/1‑70 (MQ=255)
                 aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                    >  1:4642006/1‑70 (MQ=255)
                 aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                    >  1:515958/1‑70 (MQ=255)
                  cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCt                                                                   >  1:2598920/1‑54 (MQ=255)
                      ttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:5200020/1‑71 (MQ=255)
                       tGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:4264842/1‑70 (MQ=255)
                       tGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:2379489/1‑70 (MQ=255)
                        gCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtc                                               <  1:2828131/68‑2 (MQ=16)
                        gCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                            >  1:233171/1‑70 (MQ=16)
                        gCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                            >  1:3609367/1‑70 (MQ=16)
                        gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                               <  1:4587399/68‑1 (MQ=255)
                        gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              <  1:5343728/69‑1 (MQ=255)
                        gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:4298404/1‑69 (MQ=255)
                        gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:2560355/1‑69 (MQ=255)
                         cGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccgg                                           <  1:4162183/71‑3 (MQ=16)
                         cGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                            <  1:4544287/70‑2 (MQ=16)
                         cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                            <  1:3411071/70‑1 (MQ=255)
                             agcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                                       <  1:4679745/71‑2 (MQ=16)
                             agcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGa                                          <  1:409901/68‑1 (MQ=255)
                              gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                      <  1:4628002/71‑1 (MQ=255)
                              gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                      <  1:3845462/71‑1 (MQ=255)
                               caACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtgga                                     <  1:3355794/71‑4 (MQ=12)
                               caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATc                                                         >  1:2660822/1‑51 (MQ=255)
                               caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGa                                          <  1:3765029/66‑1 (MQ=255)
                               caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                     <  1:1095740/71‑1 (MQ=255)
                                aaCGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtgg                                      >  1:4200181/1‑67 (MQ=12)
                                aaCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                     <  1:3928987/70‑1 (MQ=255)
                                 aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   >  1:1332078/1‑71 (MQ=255)
                                  cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                  >  1:5432914/1‑71 (MQ=255)
                                   ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   >  1:3842056/1‑69 (MQ=255)
                                   ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 <  1:1776993/71‑1 (MQ=255)
                                   ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 >  1:1006632/1‑71 (MQ=255)
                                    gCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTc                                <  1:480879/71‑1 (MQ=255)
                                      aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                              >  1:4784177/1‑71 (MQ=255)
                                       tCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                             >  1:4466570/1‑71 (MQ=255)
                                         gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                           >  1:2117113/1‑71 (MQ=255)
                                         gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt                           >  1:5243350/1‑71 (MQ=255)
                                          ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                          <  1:1625862/71‑1 (MQ=255)
                                            aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa                         <  1:3901193/70‑1 (MQ=255)
                                            aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                        <  1:2939575/71‑1 (MQ=255)
                                             gCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtc                                               >  1:5688232/1‑46 (MQ=16)
                                                ggaGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGAc                                         >  1:1986088/3‑50 (MQ=12)
                                                gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCa                                                  >  1:3125337/1‑41 (MQ=255)
                                                 cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGc                    <  1:4687738/70‑1 (MQ=255)
                                                 cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                   <  1:3328822/71‑1 (MQ=255)
                                                 cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                   <  1:5134130/71‑1 (MQ=255)
                                                  aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                              >  1:444222/1‑59 (MQ=255)
                                                    tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   <  1:165212/52‑1 (MQ=255)
                                                    tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg                            <  1:208285/59‑1 (MQ=255)
                                                    tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg                            <  1:159424/59‑1 (MQ=255)
                                                      tGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                            >  1:5196927/1‑40 (MQ=12)
                                                      tGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTccc                                             >  1:2778589/1‑40 (MQ=12)
                                                      tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                      <  1:809661/47‑1 (MQ=255)
                                                      tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:4119164/1‑71 (MQ=255)
                                                      tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:692386/1‑71 (MQ=255)
                                                      tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:1367844/1‑71 (MQ=255)
                                                      tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:2858834/1‑71 (MQ=255)
                                                      tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:2388103/1‑71 (MQ=255)
                                                        ttCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGTAATTGTCACGCTAcc            >  1:183813/1‑71 (MQ=255)
                                                         tCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGTAATTGTCACGCTAccg           >  1:3060230/1‑70 (MQ=255)
                                                         tCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACgg                <  1:1061429/66‑1 (MQ=255)
                                                          cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc           <  1:3835512/70‑1 (MQ=255)
                                                          cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc           <  1:5346171/70‑1 (MQ=255)
                                                          cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc           >  1:2931164/1‑70 (MQ=255)
                                                           tGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCa                               <  1:2028869/49‑1 (MQ=255)
                                                           tGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCa                               <  1:2099292/49‑1 (MQ=255)
                                                            gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg        >  1:2156103/1‑71 (MQ=255)
                                                             cACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                             >  1:1186796/1‑49 (MQ=255)
                                                             cACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                             >  1:1963865/1‑49 (MQ=255)
                                                             cACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg        <  1:4528965/70‑1 (MQ=255)
                                                               ccTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       <  1:4420204/53‑1 (MQ=255)
                                                                 tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       >  1:342486/1‑51 (MQ=255)
                                                                 tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       >  1:5405150/1‑51 (MQ=255)
                                                                 tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       >  1:1548633/1‑51 (MQ=255)
                                                                 tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       >  1:2754127/1‑51 (MQ=255)
                                                                 tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       >  1:4096186/1‑51 (MQ=255)
                                                                 tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt    >  1:4291963/1‑70 (MQ=255)
                                                                  gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa  >  1:5314578/1‑71 (MQ=255)
                                                                  gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa  >  1:5319846/1‑71 (MQ=255)
                                                                    aCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                          <  1:5263897/45‑1 (MQ=255)
                                                                      cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTgg       <  1:3187097/62‑1 (MQ=255)
                                                                      cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTgg       <  1:2092730/62‑1 (MQ=255)
                                                                      cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTgg       <  1:3717452/62‑1 (MQ=255)
                                                                      cTGCGCGATATCAGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTgg       <  1:2477612/62‑1 (MQ=255)
                                                                      |                                                                  
TCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCA  >  minE/312491‑312627

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: