Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,174 G→C 5.4% P368R (CCG→CGG)  mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,1740GC 5.4% 240.4 / 2.9 74P368R (CCG→CGG) mdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base C (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

TTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCTCAAC  >  minE/2586105‑2586244
                                                                     |                                                                      
ttGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCgg                                                                       <  1:2343186/71‑1 (MQ=255)
ttGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCgg                                                                       >  1:1471265/1‑71 (MQ=255)
    tGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGgatg                                                                   <  1:2256472/71‑1 (MQ=255)
      gTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGgatgtt                                                                 <  1:1284896/71‑3 (MQ=255)
        gAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATtg                                                               >  1:3318135/1‑71 (MQ=255)
         aaaCCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATtg                                                               <  1:3101951/70‑1 (MQ=255)
          aaCCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCgg                                                                       >  1:3704587/1‑61 (MQ=255)
          aaCCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCgg                                                                       >  1:1687016/1‑61 (MQ=255)
           accgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTgg                                                            <  1:581848/71‑1 (MQ=255)
           accgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTgg                                                            <  1:2912076/71‑1 (MQ=255)
           accgacagCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTgg                                                            >  1:1641868/1‑71 (MQ=255)
            ccgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATtg                                                               <  1:105942/67‑1 (MQ=255)
            ccgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATtg                                                               <  1:2715280/67‑1 (MQ=255)
            ccgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATtg                                                               <  1:1515990/67‑1 (MQ=255)
            ccgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATtg                                                               <  1:1243369/67‑1 (MQ=255)
            ccgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGc                                                           <  1:292813/71‑1 (MQ=255)
             cgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCa                                                          >  1:327009/1‑71 (MQ=255)
             cgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCa                                                          >  1:19074/1‑71 (MQ=255)
             cgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCa                                                          >  1:331483/1‑71 (MQ=255)
               accgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACg                                                        <  1:2218634/71‑1 (MQ=255)
                 cgccgACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGa                                                      <  1:1128298/71‑1 (MQ=255)
                   ccgACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGaaa                                                    <  1:2914440/71‑1 (MQ=255)
                   ccgACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGaaa                                                    <  1:2910231/71‑1 (MQ=255)
                   ccgACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGaaa                                                    >  1:2427262/1‑71 (MQ=255)
                     gACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAAtt                                                  <  1:3677517/71‑1 (MQ=255)
                     gACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAAtt                                                  <  1:3225248/71‑1 (MQ=255)
                      aCAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                                 >  1:2734120/1‑71 (MQ=255)
                      aCAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                                 >  1:2592374/1‑71 (MQ=255)
                      aCAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                                 <  1:949799/71‑1 (MQ=255)
                       cAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                                 <  1:2786680/70‑1 (MQ=255)
                        aaGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTg                                               >  1:975482/1‑71 (MQ=255)
                        aaGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTg                                               >  1:3164775/1‑71 (MQ=255)
                          gATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGc                                              >  1:3246673/1‑70 (MQ=255)
                            tCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                           <  1:2627765/71‑1 (MQ=255)
                            tCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaa                                            <  1:470624/70‑1 (MQ=255)
                             cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                           >  1:2340945/1‑70 (MQ=255)
                             cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                           >  1:3196857/1‑70 (MQ=255)
                             cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                           >  1:1933909/1‑70 (MQ=255)
                             cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                           >  1:1490418/1‑70 (MQ=255)
                             cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                           >  1:645919/1‑70 (MQ=255)
                             cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                           >  1:3433945/1‑70 (MQ=255)
                             cGCAAAAGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                                 >  1:476352/1‑64 (MQ=255)
                              gCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGTCACGGAAATTCTGCaaaaa                                         <  1:887505/71‑1 (MQ=255)
                              gCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaa                                          <  1:1956740/70‑1 (MQ=255)
                               cAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                         >  1:2202582/1‑70 (MQ=255)
                                aaaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                         >  1:2755802/1‑69 (MQ=255)
                                aaaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                         >  1:1479832/1‑69 (MQ=255)
                                aaaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                         >  1:2612917/1‑69 (MQ=255)
                                aaaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                         >  1:1830737/1‑69 (MQ=255)
                                aaaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                         >  1:2989303/1‑69 (MQ=255)
                                 aaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                         >  1:855571/1‑68 (MQ=255)
                                 aaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAg                                      <  1:3239687/71‑1 (MQ=255)
                                   cGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGAt                                    <  1:1361341/71‑1 (MQ=255)
                                        cGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATg                               <  1:3148438/71‑1 (MQ=255)
                                        cGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATg                               <  1:3208608/71‑1 (MQ=255)
                                           gaATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGAc                             >  1:2973894/1‑70 (MQ=255)
                                           gaATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGAc                             >  1:2884051/1‑70 (MQ=255)
                                            aaTAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                         >  1:307445/1‑57 (MQ=255)
                                               aaCCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAgc                         >  1:1516235/1‑69 (MQ=255)
                                                accaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGaa                       <  1:1832919/71‑1 (MQ=255)
                                                accaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGaa                       <  1:3106758/71‑1 (MQ=255)
                                                accaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGaa                       <  1:908787/71‑1 (MQ=255)
                                                accaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGaa                       <  1:3143398/71‑1 (MQ=255)
                                                 ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                                 <  1:887869/44‑1 (MQ=255)
                                                 ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                                 <  1:1514033/44‑1 (MQ=255)
                                                 ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                                 <  1:3178140/44‑1 (MQ=255)
                                                 ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                                 <  1:287674/44‑1 (MQ=255)
                                                 ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAgg                         >  1:3509186/1‑68 (MQ=255)
                                                 ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAgg                         >  1:913836/1‑68 (MQ=255)
                                                    ccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGAc                    >  1:1235509/1‑70 (MQ=255)
                                                     cGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGAc                    >  1:2067913/1‑69 (MQ=255)
                                                      ggTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAg                 <  1:2921220/71‑1 (MQ=255)
                                                       gTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAg                 >  1:972558/1‑70 (MQ=255)
                                                       gTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAg                 >  1:413418/1‑70 (MQ=255)
                                                           cacAATCATCCCGACGATTGTGGCACGGAAATTCTGCAAAAACAGATACAt                                <  1:813345/49‑1 (MQ=21)
                                                           cacAATCATCCCGACGATTGTGGCACGGAAATTCTGCAAAAACAGATACAt                                <  1:3489275/49‑1 (MQ=21)
                                                             cAATCATCCCGACGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGa                              >  1:2843095/1‑51 (MQ=25)
                                                             cAATCATCCCGACGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGa                              >  1:3098548/1‑51 (MQ=25)
                                                                tcgtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCa                           <  1:1750214/51‑1 (MQ=255)
                                                                tcgtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGcc       <  1:939492/71‑1 (MQ=255)
                                                                 cgtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATTACCAGGAAGACTAGGATGATAGCCt      >  1:68112/1‑71 (MQ=255)
                                                                 cgtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGcc       <  1:968366/70‑1 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTTCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCt      >  1:1127860/1‑70 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAAc                                        >  1:450100/1‑36 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAg                          >  1:1039440/1‑50 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAg                          >  1:783366/1‑50 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAg                          >  1:1096080/1‑50 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAg                          >  1:2784802/1‑50 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAg                          >  1:1362154/1‑50 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGgatg             <  1:170488/63‑1 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGgatg             <  1:2165725/63‑1 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGgatg             <  1:2489156/63‑1 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGgatg             <  1:647807/63‑1 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGgatg             <  1:2828473/63‑1 (MQ=255)
                                                                  gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCt      >  1:2387138/1‑70 (MQ=255)
                                                                     gggATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGc        <  1:2736986/65‑1 (MQ=255)
                                                                     gggATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCTCAAc  >  1:2055405/1‑71 (MQ=255)
                                                                     gggATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCTCAAc  >  1:2041426/1‑71 (MQ=255)
                                                                     gggATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCTCAAc  <  1:311468/71‑1 (MQ=255)
                                                                     |                                                                      
TTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCTCAAC  >  minE/2586105‑2586244

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: