Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,067 A→G 11.1% L404L (TTA→CTA) ‡ mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0670AG11.1% 317.7 / 24.8 117L404L (TTA→CTA) ‡mdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base A (43/61);  new base G (8/5);  total (51/66)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

Only 100 of 122 total aligned reads displayed.

GGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCA  >  minE/2585997‑2586137
                                                                      |                                                                      
ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa                                                                        >  1:3302608/1‑71 (MQ=255)
    ggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                                                                     >  1:1364286/3‑70 (MQ=255)
     tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                                   <  1:290303/71‑1 (MQ=255)
     tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                                   <  1:3595756/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:1274685/1‑70 (MQ=255)
          tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                               >  1:3234193/1‑70 (MQ=255)
           aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGcc                                                              <  1:3554745/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)
                  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:1467828/1‑56 (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: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)
                            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)
                            gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                              >  1:944151/1‑69 (MQ=255)
                            gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                              >  1:1872219/1‑69 (MQ=255)
                              tCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                                          >  1:187113/1‑71 (MQ=255)
                               cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                               <  1:928251/65‑1 (MQ=255)
                               cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                               <  1:413724/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:2210775/71‑1 (MQ=255)
                                   gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                     <  1:2033255/71‑1 (MQ=255)
                                   gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                     <  1:2441673/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)
                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                             <  1:836523/56‑1 (MQ=255)
                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                             <  1:1956328/56‑1 (MQ=255)
                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                              <  1:916688/71‑1 (MQ=255)
                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                              >  1:3757272/1‑71 (MQ=255)
                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                              >  1:3481322/1‑71 (MQ=255)
                                            caccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                                    <  1:1421996/63‑1 (MQ=255)
                                               cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                         <  1:2792971/71‑1 (MQ=255)
                                               cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                         <  1:3618446/71‑1 (MQ=255)
                                                acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                        <  1:40117/71‑1 (MQ=255)
                                                acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                        <  1:2680586/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)
                                                   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)
                                                         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:1729638/1‑70 (MQ=255)
                                                           aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc              <  1:1155688/70‑1 (MQ=255)
                                                           aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc              >  1:1610756/1‑70 (MQ=255)
                                                              caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc           <  1:923957/70‑1 (MQ=255)
                                                              caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc           <  1:1942710/70‑1 (MQ=255)
                                                               accaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgc               <  1:625019/65‑1 (MQ=255)
                                                               accaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACaa         <  1:673907/71‑1 (MQ=255)
                                                                ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                                         <  1:2271318/38‑1 (MQ=255)
                                                                ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc           >  1:2024760/1‑68 (MQ=255)
                                                                ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAg        <  1:694486/71‑1 (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:3929530/71‑1 (MQ=255)
                                                                 caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa       <  1:852538/71‑1 (MQ=255)
                                                                 caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa       <  1:2857097/71‑1 (MQ=255)
                                                                 caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa       <  1:2619898/71‑1 (MQ=255)
                                                                  aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTg                                >  1:664354/1‑45 (MQ=255)
                                                                  aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTg                                >  1:146125/1‑45 (MQ=255)
                                                                      aCCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGc   <  1:356836/70‑1 (MQ=255)
                                                                      aCCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCa  <  1:3838294/71‑1 (MQ=255)
                                                                      |                                                                      
GGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCA  >  minE/2585997‑2586137

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: