Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,059 C→G 13.3% V406V (GTG→GTC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0590CG13.3% 296.2 / 26.1 105V406V (GTG→GTCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base C (39/52);  new base G (8/6);  total (47/58)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

CTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAC  >  minE/2585989‑2586128
                                                                      |                                                                     
cttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcc                                                                       <  1:576957/71‑1 (MQ=255)
cttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcc                                                                       <  1:2360006/71‑1 (MQ=255)
  tccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                                     >  1:1285240/1‑71 (MQ=255)
   ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTTATccac                                                                     <  1:511122/70‑1 (MQ=255)
   ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                                     <  1:2412808/70‑1 (MQ=255)
       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)
           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)
                        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)
                            aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                           >  1:253318/1‑71 (MQ=255)
                            aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                           >  1:3364156/1‑71 (MQ=255)
                              tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                        >  1:3748009/1‑56 (MQ=255)
                               gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                                <  1:786780/63‑1 (MQ=255)
                                acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATGCCGaa                                       >  1:2101319/4‑71 (MQ=2)
                                acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                       >  1:1833409/4‑71 (MQ=2)
                                acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                       >  1:1284789/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:3711481/57‑1 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                 <  1:148099/57‑1 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                 <  1:1200652/57‑1 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                                                                 <  1:1123023/41‑1 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                     >  1:2023796/1‑69 (MQ=255)
                                    gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                     >  1:1872219/1‑69 (MQ=255)
                                     cTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                                            <  1:2353192/45‑1 (MQ=255)
                                      tCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                                 >  1:187113/1‑71 (MQ=255)
                                       cGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGt                                 >  1:3416428/1‑70 (MQ=2)
                                       cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                      <  1:413724/65‑1 (MQ=255)
                                       cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                      <  1:928251/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:1140588/1‑71 (MQ=255)
                                           gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                            <  1:2033255/71‑1 (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:2685255/1‑71 (MQ=255)
                                                  accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                     <  1:916688/71‑1 (MQ=255)
                                                  accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                     >  1:3757272/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:3900925/71‑1 (MQ=255)
                                                        acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa               <  1:2680586/71‑1 (MQ=255)
                                                        acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa               <  1:1450063/71‑1 (MQ=255)
                                                        acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa               <  1:40117/71‑1 (MQ=255)
                                                         cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtgt                         >  1:1537788/1‑60 (MQ=255)
                                                         cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa              <  1:2915695/71‑1 (MQ=255)
                                                         cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa              <  1:484202/71‑1 (MQ=255)
                                                          gATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                                >  1:124831/1‑52 (MQ=255)
                                                           aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc            <  1:3423260/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)
                                                                 tCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATg                    <  1:2577487/57‑1 (MQ=255)
                                                                 tCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgc      <  1:3925589/71‑1 (MQ=255)
                                                                  cATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                           >  1:3763548/1‑49 (MQ=255)
                                                                  cATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc     >  1:2434099/1‑71 (MQ=255)
                                                                   aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                    >  1:178616/1‑39 (MQ=255)
                                                                   aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAac             >  1:212543/1‑62 (MQ=255)
                                                                   aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc     >  1:1610756/1‑70 (MQ=255)
                                                                   aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc     <  1:1155688/70‑1 (MQ=255)
                                                                    tCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGa   >  1:3519300/1‑71 (MQ=255)
                                                                      caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc  <  1:1942710/70‑1 (MQ=255)
                                                                      caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc  <  1:923957/70‑1 (MQ=255)
                                                                      |                                                                     
CTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAC  >  minE/2585989‑2586128

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: