Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,706 C→G 17.0% S389S (TCG→TCC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7060CG17.0% 244.4 / 34.0 94S389S (TCG→TCCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (29/49);  new base G (10/6);  total (39/55)
Fisher's exact test for biased strand distribution p-value = 9.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAAT  >  minE/311636‑311774
                                                                      |                                                                    
ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATc                                                                      <  1:1727597/71‑1 (MQ=255)
 cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                                                                     <  1:564412/71‑1 (MQ=255)
 cacaACGATCGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                                                                     <  1:1458078/71‑1 (MQ=255)
   caACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag                                                                   <  1:2228932/71‑1 (MQ=255)
    aaCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga                                                                  <  1:2189356/71‑1 (MQ=255)
         tGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                              <  1:2971964/70‑1 (MQ=255)
            cGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAg                                                          <  1:3797414/71‑1 (MQ=255)
             atCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaaaag                                                         >  1:2778838/2‑69 (MQ=9)
             atCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaaaag                                                         >  1:3276913/2‑69 (MQ=9)
              tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                                                                     <  1:1899985/58‑1 (MQ=255)
               cATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAg                                                                >  1:227123/1‑62 (MQ=255)
                aTCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                                                           >  1:563275/1‑66 (MQ=255)
                 tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                     <  1:2575103/71‑1 (MQ=255)
                 tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                     <  1:1824706/71‑1 (MQ=255)
                      caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                <  1:2470315/71‑1 (MQ=255)
                       agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagcggc                                                     <  1:1648596/63‑1 (MQ=12)
                       agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                               >  1:2319320/3‑70 (MQ=12)
                       agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                               <  1:2299073/69‑2 (MQ=12)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGaa                                                                 >  1:2471585/1‑53 (MQ=255)
                       aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                               <  1:2742190/71‑1 (MQ=255)
                        gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGatg                                              >  1:3815727/2‑69 (MQ=12)
                        gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGatg                                              >  1:576654/2‑69 (MQ=12)
                        aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                                              >  1:3497243/1‑71 (MQ=255)
                          aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACgg                                             >  1:3534739/1‑70 (MQ=255)
                           ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGc                                                               <  1:3191325/51‑1 (MQ=255)
                            gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                          <  1:1796073/71‑1 (MQ=255)
                            gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                          <  1:3864917/71‑1 (MQ=255)
                            gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                          <  1:2124814/71‑1 (MQ=255)
                            gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                          <  1:2363271/71‑1 (MQ=255)
                            gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                           <  1:3006533/70‑1 (MQ=255)
                            cccGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCa                                           <  1:690757/69‑1 (MQ=12)
                             ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga                                                                  <  1:899035/46‑1 (MQ=255)
                             ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                           <  1:1565651/69‑1 (MQ=255)
                              cGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCaaac                                        <  1:3652167/71‑2 (MQ=16)
                               gATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                                            >  1:1247340/1‑50 (MQ=9)
                               gATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                                            >  1:2075371/1‑50 (MQ=9)
                               gATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                                            >  1:892493/1‑50 (MQ=9)
                                aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt                                      <  1:3732859/71‑1 (MQ=255)
                                aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt                                      >  1:1668109/1‑71 (MQ=255)
                                aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt                                      >  1:820374/1‑71 (MQ=255)
                                 tGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTc                                     >  1:2077061/1‑71 (MQ=255)
                                 tGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTc                                     >  1:1306873/1‑71 (MQ=255)
                                 tGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagcgg                                                      >  1:2130779/1‑54 (MQ=9)
                                 tGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagcgg                                                      >  1:3950389/1‑54 (MQ=9)
                                  ggCTAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                                         >  1:192406/1‑50 (MQ=255)
                                  ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                                    >  1:3435692/1‑71 (MQ=255)
                                   gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                >  1:2907973/1‑58 (MQ=255)
                                   gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                                              >  1:762918/1‑60 (MQ=255)
                                   gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                   <  1:2262432/71‑1 (MQ=255)
                                   gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                   <  1:3526958/71‑1 (MQ=255)
                                     aaGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                               <  1:663174/56‑2 (MQ=9)
                                      aGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                             <  1:772822/42‑1 (MQ=255)
                                       gCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                                         >  1:1582607/1‑45 (MQ=255)
                                       gCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAg                                >  1:2540071/1‑70 (MQ=255)
                                        cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                                              >  1:992243/1‑55 (MQ=255)
                                        cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                                              >  1:3730817/1‑55 (MQ=255)
                                         aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCaa                             <  1:3568196/71‑1 (MQ=255)
                                         aCCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaa                             <  1:3435123/71‑1 (MQ=12)
                                          ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                     <  1:1463810/46‑1 (MQ=255)
                                            aTCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGa                                 >  1:3509140/1‑64 (MQ=255)
                                             tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggc                                                        <  1:170205/40‑1 (MQ=255)
                                             tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                                            <  1:2946875/52‑1 (MQ=255)
                                             tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                                            <  1:865345/52‑1 (MQ=255)
                                             tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATac                          <  1:1866087/70‑1 (MQ=255)
                                               ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                                    <  1:1104027/58‑1 (MQ=255)
                                               ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATacca                        <  1:742993/70‑1 (MQ=255)
                                                cGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                        <  1:2476207/53‑1 (MQ=255)
                                                cGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                        <  1:3615318/53‑1 (MQ=255)
                                                cGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                        <  1:832086/53‑1 (MQ=255)
                                                 gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcg                                                       <  1:2982814/37‑1 (MQ=255)
                                                 gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGGCGGCAAAGGTCCCGAGCAATACCACCg                     >  1:3583996/1‑71 (MQ=255)
                                                 gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                     <  1:3451732/71‑1 (MQ=255)
                                                 gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                     <  1:2780201/71‑1 (MQ=255)
                                                  aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                    <  1:2634466/71‑1 (MQ=255)
                                                  aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                    <  1:1878588/71‑1 (MQ=255)
                                                   aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                    <  1:899232/70‑1 (MQ=255)
                                                   aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                    >  1:1436073/1‑70 (MQ=255)
                                                    cATTGTTTGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTa                  <  1:716698/71‑1 (MQ=255)
                                                    cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTa                  >  1:41513/1‑71 (MQ=255)
                                                    cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTa                  >  1:313318/1‑71 (MQ=255)
                                                     aTTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTAc                 >  1:3030104/1‑71 (MQ=255)
                                                       tGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGa                                 >  1:774512/1‑53 (MQ=255)
                                                        gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACgg               >  1:537101/1‑70 (MQ=255)
                                                         ttAGCGTGTTTATCGAGAAGCCAAAGGCGGCTAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGc              <  1:1383837/70‑1 (MQ=255)
                                                         ttAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                                    <  1:1533185/48‑1 (MQ=255)
                                                          tAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCaa                             >  1:1543152/1‑54 (MQ=255)
                                                           aGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATaccac                       >  1:1478909/1‑59 (MQ=255)
                                                             cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATgg         <  1:2288479/71‑1 (MQ=255)
                                                             cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATgg         <  1:3819116/71‑1 (MQ=255)
                                                             cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATgg         <  1:1681432/71‑1 (MQ=255)
                                                              gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGt        >  1:1594944/1‑71 (MQ=255)
                                                               tgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                           <  1:621452/35‑1 (MQ=255)
                                                                   tATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGa                                 <  1:2700496/41‑1 (MQ=255)
                                                                    aTCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGaa   >  1:300351/1‑70 (MQ=255)
                                                                     tCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAAt  <  1:2453210/70‑1 (MQ=255)
                                                                      |                                                                    
CCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAAT  >  minE/311636‑311774

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: