Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,062 C→G 13.3% L405L (CTG→CTC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0620CG13.3% 293.8 / 27.4 105L405L (CTG→CTCmdtFmultidrug 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 = 9.74e-01

Only 100 of 123 total aligned reads displayed.

TCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGA  >  minE/2585995‑2586132
                                                                   |                                                                      
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)
         gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                            <  1:1975963/71‑1 (MQ=255)
         gCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGAtg                                                            >  1:3295347/1‑70 (MQ=255)
            tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                          >  1:1274685/1‑70 (MQ=255)
            tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                          >  1:3234193/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)
                  ccGCAATGACACTCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                                     <  1:1871226/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:3748009/1‑56 (MQ=255)
                        tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                            >  1:1467828/1‑56 (MQ=255)
                        tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc                                              >  1:416224/1‑70 (MQ=255)
                         gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                                    <  1:786780/63‑1 (MQ=255)
                          acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                           >  1:1833409/4‑71 (MQ=2)
                          acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                           >  1:3790593/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)
                              gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                     <  1:148099/57‑1 (MQ=255)
                              gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                     <  1:1200652/57‑1 (MQ=255)
                              gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                                                                     <  1:1123023/41‑1 (MQ=255)
                              gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                         >  1:944151/1‑69 (MQ=255)
                              gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                         >  1:1872219/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)
                                  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:417058/1‑71 (MQ=255)
                                    cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                                 <  1:3045170/71‑1 (MQ=255)
                                     gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                >  1:1140588/1‑71 (MQ=255)
                                     gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                <  1:2033255/71‑1 (MQ=255)
                                     gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                <  1:2210775/71‑1 (MQ=255)
                                      ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                               >  1:1738287/1‑71 (MQ=255)
                                       tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtg                              >  1:142341/1‑71 (MQ=255)
                                       tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                               <  1:526951/70‑1 (MQ=255)
                                            accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                        <  1:836523/56‑1 (MQ=255)
                                            accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                        <  1:1956328/56‑1 (MQ=255)
                                            accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                         >  1:3481322/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:40117/71‑1 (MQ=255)
                                                  acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                   <  1:1450063/71‑1 (MQ=255)
                                                  acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                   <  1:3900925/71‑1 (MQ=255)
                                                   cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                  <  1:484202/71‑1 (MQ=255)
                                                   cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                  <  1:2915695/71‑1 (MQ=255)
                                                     aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc                <  1:3099023/71‑1 (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)
                                                             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:923957/70‑1 (MQ=255)
                                                                 accaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgc          <  1:625019/65‑1 (MQ=255)
                                                                 accaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACaa    <  1:673907/71‑1 (MQ=255)
                                                                  ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc      >  1:2024760/1‑68 (MQ=255)
                                                                   caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa                          <  1:738470/47‑1 (MQ=255)
                                                                   caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                  >  1:522119/1‑55 (MQ=255)
                                                                   caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                  >  1:417796/1‑55 (MQ=255)
                                                                   caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa  <  1:2857097/71‑1 (MQ=255)
                                                                   caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa  <  1:852538/71‑1 (MQ=255)
                                                                   caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa  <  1:2619898/71‑1 (MQ=255)
                                                                   caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa  <  1:1782093/71‑1 (MQ=255)
                                                                   caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa  <  1:3929530/71‑1 (MQ=255)
                                                                   |                                                                      
TCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGA  >  minE/2585995‑2586132

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: