Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,586,171 G→A 5.1% T369M (ACG→ATG)  mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,586,1710GA 5.1% 263.5 / 2.9 78T369M (ACG→ATG) mdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base A (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 101 total aligned reads displayed.

CGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCT  >  Exported/2586101‑2586240
                                                                      |                                                                     
cGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAAtcg                                                                       >  1:1219027/1‑71 (MQ=255)
cGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAAtcg                                                                       <  1:1343396/71‑1 (MQ=255)
 gtgtTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAAtcgt                                                                      <  1:452519/71‑1 (MQ=255)
 gtgtTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAAtcgt                                                                      <  1:1131262/71‑1 (MQ=255)
  tgtTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAAtcgtc                                                                     <  1:2491027/71‑1 (MQ=255)
    ttGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCgg                                                                   >  1:1471265/1‑71 (MQ=255)
    ttGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCgg                                                                   <  1:2343186/71‑1 (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)
              aaCCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAAtcgtc                                                                     >  1:3550194/1‑59 (MQ=255)
              aaCCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCgg                                                                   >  1:3704587/1‑61 (MQ=255)
              aaCCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCgg                                                                   >  1:1687016/1‑61 (MQ=255)
               accgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTgg                                                        <  1:2912076/71‑1 (MQ=255)
               accgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTgg                                                        <  1:581848/71‑1 (MQ=255)
               accgacagCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTgg                                                        >  1:1641868/1‑71 (MQ=255)
                ccgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATtg                                                           <  1:2715280/67‑1 (MQ=255)
                ccgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATtg                                                           <  1:105942/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:331483/1‑71 (MQ=255)
                 cgaccgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCa                                                      >  1:19074/1‑71 (MQ=255)
                   accgCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACg                                                    <  1:2218634/71‑1 (MQ=255)
                     cgccgACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGa                                                  <  1:1128298/71‑1 (MQ=255)
                       ccgACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGaaa                                                <  1:2910231/71‑1 (MQ=255)
                       ccgACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGaaa                                                <  1:2914440/71‑1 (MQ=255)
                       ccgACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGaaa                                                >  1:2427262/1‑71 (MQ=255)
                         gACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAAtt                                              <  1:3225248/71‑1 (MQ=255)
                         gACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAAtt                                              <  1:3677517/71‑1 (MQ=255)
                          aCAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                             <  1:949799/71‑1 (MQ=255)
                          aCAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                             >  1:2592374/1‑71 (MQ=255)
                          aCAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                             >  1:2734120/1‑71 (MQ=255)
                           cAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                             <  1:2786680/70‑1 (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:645919/1‑70 (MQ=255)
                                 cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                       >  1:1490418/1‑70 (MQ=255)
                                 cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                       >  1:1933909/1‑70 (MQ=255)
                                 cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                       >  1:3433945/1‑70 (MQ=255)
                                 cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                       >  1:2340945/1‑70 (MQ=255)
                                 cGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaa                                       >  1:3196857/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:2989303/1‑69 (MQ=255)
                                    aaaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                     >  1:2612917/1‑69 (MQ=255)
                                    aaaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                     >  1:2755802/1‑69 (MQ=255)
                                    aaaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                     >  1:1479832/1‑69 (MQ=255)
                                    aaaCGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                     >  1:1830737/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:3208608/71‑1 (MQ=255)
                                            cGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATg                           <  1:3148438/71‑1 (MQ=255)
                                               gaATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGAc                         >  1:2884051/1‑70 (MQ=255)
                                               gaATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGAc                         >  1:2973894/1‑70 (MQ=255)
                                                aaTAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCaaaaa                                     >  1:307445/1‑57 (MQ=255)
                                                   aaCCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAgc                     >  1:1516235/1‑69 (MQ=255)
                                                    accaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGaa                   <  1:3143398/71‑1 (MQ=255)
                                                    accaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGaa                   <  1:1832919/71‑1 (MQ=255)
                                                    accaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGaa                   <  1:908787/71‑1 (MQ=255)
                                                    accaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGaa                   <  1:3106758/71‑1 (MQ=255)
                                                     ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                             <  1:287674/44‑1 (MQ=255)
                                                     ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                             <  1:3178140/44‑1 (MQ=255)
                                                     ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                             <  1:887869/44‑1 (MQ=255)
                                                     ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTc                                             <  1:1514033/44‑1 (MQ=255)
                                                     ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAgg                     >  1:913836/1‑68 (MQ=255)
                                                     ccaccGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAgg                     >  1:3509186/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:413418/1‑70 (MQ=255)
                                                           gTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAg             >  1:972558/1‑70 (MQ=255)
                                                               cacAATCATCCCGACGATTGTGGCACGGAAATTCTGCAAAAACAGATACAt                            <  1:3489275/49‑1 (MQ=21)
                                                               cacAATCATCCCGACGATTGTGGCACGGAAATTCTGCAAAAACAGATACAt                            <  1:813345/49‑1 (MQ=21)
                                                                 cAATCATCCCGACGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGa                          >  1:3098548/1‑51 (MQ=25)
                                                                 cAATCATCCCGACGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGa                          >  1:2843095/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:1362154/1‑50 (MQ=255)
                                                                      gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAg                      >  1:1096080/1‑50 (MQ=255)
                                                                      gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAg                      >  1:783366/1‑50 (MQ=255)
                                                                      gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAg                      >  1:2784802/1‑50 (MQ=255)
                                                                      gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAg                      >  1:1039440/1‑50 (MQ=255)
                                                                      gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGgatg         <  1:2165725/63‑1 (MQ=255)
                                                                      gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGgatg         <  1:170488/63‑1 (MQ=255)
                                                                      gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGgatg         <  1:2828473/63‑1 (MQ=255)
                                                                      gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGgatg         <  1:647807/63‑1 (MQ=255)
                                                                      gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGgatg         <  1:2489156/63‑1 (MQ=255)
                                                                      gtcgGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCt  >  1:2387138/1‑70 (MQ=255)
                                                                      |                                                                     
CGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCGAGAATAACCACCGGTACGGCAATCGTCGGGATGATTGTGGCACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCT  >  Exported/2586101‑2586240

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: