Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,558 A→C 17.0% V105V (GTT→GTG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5580AC17.0% 220.3 / 44.8 106V105V (GTT→GTGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (55/33);  new base C (12/6);  total (67/39)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 112 total aligned reads displayed.

CGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCA  >  minE/312488‑312627
                                                                      |                                                                     
cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                                                                       <  1:4243671/71‑1 (MQ=255)
cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                                                                       <  1:5650579/71‑1 (MQ=255)
cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                                                                       <  1:1460206/71‑1 (MQ=255)
 gCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGaa                                                                      >  1:2302376/1‑71 (MQ=255)
  cTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                                     >  1:5181787/1‑71 (MQ=255)
  cTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                                     >  1:4415100/1‑71 (MQ=255)
  cTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                                     >  1:5188030/1‑71 (MQ=255)
   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:3131442/1‑71 (MQ=255)
        ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                                               <  1:1253011/71‑1 (MQ=255)
        ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcg                                                                >  1:630402/1‑70 (MQ=255)
          ttGCTGCTCAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                                                                       <  1:2550526/61‑1 (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)
                      ttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                                                                       <  1:2553563/49‑1 (MQ=255)
                       tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                                     >  1:1774050/1‑50 (MQ=255)
                       tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                                     >  1:1260472/1‑50 (MQ=255)
                        cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                 <  1:4296053/69‑1 (MQ=255)
                         ttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:5200020/1‑71 (MQ=255)
                          tGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:2379489/1‑70 (MQ=255)
                          tGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              >  1:4264842/1‑70 (MQ=255)
                           gCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtc                                               <  1:2828131/68‑2 (MQ=16)
                           gCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                            >  1:233171/1‑70 (MQ=16)
                           gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                              <  1:5343728/69‑1 (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)
                                     cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                  >  1:5432914/1‑71 (MQ=255)
                                      ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   >  1:3842056/1‑69 (MQ=255)
                                      ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 >  1:1006632/1‑71 (MQ=255)
                                      ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                 <  1:1776993/71‑1 (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: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)
                                               aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                        <  1:3479724/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:5134130/71‑1 (MQ=255)
                                                    cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                   <  1:3328822/71‑1 (MQ=255)
                                                     aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg                              >  1:444222/1‑59 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                   <  1:165212/52‑1 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg                            <  1:159424/59‑1 (MQ=255)
                                                       tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg                            <  1:208285/59‑1 (MQ=255)
                                                         tGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                            >  1:1841353/1‑40 (MQ=12)
                                                         tGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                                            >  1:5196927/1‑40 (MQ=12)
                                                         tGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTccc                                             >  1:2778589/1‑40 (MQ=12)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                      <  1:809661/47‑1 (MQ=255)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:692386/1‑71 (MQ=255)
                                                         tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa              >  1:4119164/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)
                                                            tCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACgg                <  1:1061429/66‑1 (MQ=255)
                                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc           <  1:3835512/70‑1 (MQ=255)
                                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc           >  1:2931164/1‑70 (MQ=255)
                                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg          >  1:515794/1‑71 (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:1963865/1‑49 (MQ=255)
                                                                cACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg                             >  1:1186796/1‑49 (MQ=255)
                                                                cACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg        <  1:4528965/70‑1 (MQ=255)
                                                                  ccTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       <  1:4420204/53‑1 (MQ=255)
                                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       >  1:2754127/1‑51 (MQ=255)
                                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       >  1:5405150/1‑51 (MQ=255)
                                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       >  1:342486/1‑51 (MQ=255)
                                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATc                       >  1:1548633/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)
                                                                      |                                                                     
CGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCA  >  minE/312488‑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: