Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,053 G→A 13.7% D408D (GAC→GAT mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0530GA13.7% 291.2 / 26.9 102D408D (GAC→GATmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base G (34/54);  new base A (8/6);  total (42/60)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

CTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGA  >  minE/2585989‑2586120
                                                                |                                                                   
cttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcc                                                               <  1:2360006/71‑1 (MQ=255)
cttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcc                                                               <  1:576957/71‑1 (MQ=255)
  tccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                             >  1:625026/1‑71 (MQ=255)
  tccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                             >  1:1285240/1‑71 (MQ=255)
   ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTTATccac                                                             <  1:511122/70‑1 (MQ=255)
   ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                             <  1:2412808/70‑1 (MQ=255)
      ttggcggGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAg                                                          <  1:2026811/68‑1 (MQ=12)
       cggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCGCCAGt                                                         >  1:1686169/1‑70 (MQ=255)
        ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa                                                       >  1:3302608/1‑71 (MQ=255)
        ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                                                        <  1:1024032/70‑1 (MQ=255)
        ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                                                        <  1:1783259/70‑1 (MQ=255)
        ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                                                        <  1:3740850/70‑1 (MQ=255)
         gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                      >  1:3736467/1‑71 (MQ=255)
           ggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                           <  1:1971060/64‑1 (MQ=255)
            ggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                                                    >  1:1364286/3‑70 (MQ=255)
             tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                  <  1:3595756/71‑1 (MQ=255)
             tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                  <  1:290303/71‑1 (MQ=255)
              aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                 >  1:3489502/1‑71 (MQ=255)
               gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                <  1:1975963/71‑1 (MQ=255)
               gCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGAtg                                                >  1:3295347/1‑70 (MQ=255)
                 ttATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                              >  1:2878324/1‑71 (MQ=255)
                  tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                              >  1:3234193/1‑70 (MQ=255)
                  tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                              >  1:1274685/1‑70 (MQ=255)
                   aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                                    <  1:1540280/63‑1 (MQ=255)
                   aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCg                                               <  1:2119080/68‑1 (MQ=255)
                   aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGcc                                             <  1:3554745/70‑1 (MQ=255)
                    tCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCa                                            <  1:2884450/70‑1 (MQ=255)
                      ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATACACCAGTAACCCTATCGCCAGCa                                         >  1:3488866/1‑71 (MQ=255)
                       tCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                        >  1:521560/1‑71 (MQ=255)
                        ccGCAATGACACTCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                         <  1:1871226/69‑1 (MQ=255)
                        ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                         <  1:3080823/69‑1 (MQ=255)
                         cGCAATTACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                       <  1:2019901/70‑1 (MQ=255)
                         cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                       <  1:59734/70‑1 (MQ=255)
                         cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                      <  1:3776617/71‑1 (MQ=255)
                          gCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                 >  1:179953/1‑59 (MQ=255)
                            aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa                                                       >  1:2727636/1‑51 (MQ=255)
                            aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                  <  1:1279990/56‑1 (MQ=255)
                            aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                   >  1:253318/1‑71 (MQ=255)
                            aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                   >  1:3364156/1‑71 (MQ=255)
                              tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                >  1:1467828/1‑56 (MQ=255)
                              tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                >  1:3748009/1‑56 (MQ=255)
                              tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc                                  >  1:416224/1‑70 (MQ=255)
                               gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                        <  1:786780/63‑1 (MQ=255)
                                acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATGCCGaa                               >  1:2101319/4‑71 (MQ=2)
                                acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                               >  1:1284789/4‑71 (MQ=2)
                                acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                               >  1:1833409/4‑71 (MQ=2)
                                acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                               >  1:3790593/4‑71 (MQ=2)
                                   cGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACa                             >  1:1706856/1‑70 (MQ=2)
                                   cGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                            <  1:3795494/71‑1 (MQ=255)
                                    gCTCGACGTTCTCCACCACCGCGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg                                 >  1:3353477/1‑65 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                         <  1:1200652/57‑1 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                         <  1:3711481/57‑1 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                         <  1:148099/57‑1 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                                                         <  1:1123023/41‑1 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                             >  1:1872219/1‑69 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                             >  1:2023796/1‑69 (MQ=255)
                                     cTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                                    <  1:2353192/45‑1 (MQ=255)
                                      tCGACGTTCTCCACCACCACGATGGCGTCATccacc                                                            <  1:2801586/36‑1 (MQ=40)
                                      tCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                         >  1:187113/1‑71 (MQ=255)
                                       cGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGt                         >  1:3416428/1‑70 (MQ=2)
                                       cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                              <  1:928251/65‑1 (MQ=255)
                                       cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                              <  1:413724/65‑1 (MQ=255)
                                        gACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGtta                       <  1:85629/71‑3 (MQ=2)
                                        gACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGtta                       <  1:2084177/71‑3 (MQ=2)
                                        gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                         <  1:2651753/69‑1 (MQ=255)
                                        gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                        <  1:2405175/70‑1 (MQ=255)
                                        gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCa                       <  1:2225889/71‑1 (MQ=255)
                                         aCGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAg                          >  1:2130414/1‑67 (MQ=255)
                                          cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                              >  1:2417964/1‑62 (MQ=255)
                                          cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                     <  1:3045170/71‑1 (MQ=255)
                                          cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                     >  1:417058/1‑71 (MQ=255)
                                           gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                    <  1:2210775/71‑1 (MQ=255)
                                           gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                    <  1:2441673/71‑1 (MQ=255)
                                           gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                    >  1:1140588/1‑71 (MQ=255)
                                            ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                     <  1:1396246/69‑1 (MQ=255)
                                            ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                   >  1:1738287/1‑71 (MQ=255)
                                             tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                           <  1:3877786/46‑1 (MQ=255)
                                             tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtg                  >  1:142341/1‑71 (MQ=255)
                                             tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                   <  1:526951/70‑1 (MQ=255)
                                                  accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                            <  1:1956328/56‑1 (MQ=255)
                                                  accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                            <  1:836523/56‑1 (MQ=255)
                                                  accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt             >  1:3757272/1‑71 (MQ=255)
                                                  accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt             <  1:916688/71‑1 (MQ=255)
                                                  accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt             >  1:2685255/1‑71 (MQ=255)
                                                    caccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                   <  1:1421996/63‑1 (MQ=255)
                                                       cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa        <  1:3618446/71‑1 (MQ=255)
                                                       cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa        <  1:2792971/71‑1 (MQ=255)
                                                        acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa       <  1:1450063/71‑1 (MQ=255)
                                                        acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa       <  1:3900925/71‑1 (MQ=255)
                                                        acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa       <  1:2680586/71‑1 (MQ=255)
                                                        acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa       <  1:40117/71‑1 (MQ=255)
                                                         cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtgt                 >  1:1537788/1‑60 (MQ=255)
                                                         cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa      <  1:484202/71‑1 (MQ=255)
                                                         cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa      <  1:2915695/71‑1 (MQ=255)
                                                          gATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                        >  1:124831/1‑52 (MQ=255)
                                                           aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc    <  1:3423260/71‑1 (MQ=255)
                                                           aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc    <  1:3099023/71‑1 (MQ=255)
                                                            tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccg   <  1:3864657/71‑1 (MQ=255)
                                                              gCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccga  <  1:1753339/70‑1 (MQ=255)
                                                                gTCATCCACCAGTAACCCTATCGCCAGCACCATCTCGAACATAGTCAAc                     <  1:3740599/49‑1 (MQ=255)
                                                                |                                                                   
CTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGA  >  minE/2585989‑2586120

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: