Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,872,673 C→A 33.1% N432K (AAC→AAA)  ygcE → predicted kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,872,6730CA33.1% 135.2 / 73.9 116N432K (AAC→AAA) ygcEpredicted kinase
Reads supporting (aligned to +/- strand):  ref base C (20/56);  new base A (31/7);  total (52/64)
Fisher's exact test for biased strand distribution p-value = 2.33e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 116 total aligned reads displayed.

GGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAT  >  minE/1872628‑1872739
                                             |                                                                  
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAAGATGGTTCGCGTGAAAGATATCtt                                             <  1:2001960/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:1014479/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:271771/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:2377137/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:1180047/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:926472/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:356783/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:2260819/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:2026466/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:1981500/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:470986/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:1667164/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:1634574/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:1519356/69‑1 (MQ=255)
ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGATCGCGTGAAAGATATCtt                                             <  1:2061144/69‑1 (MQ=255)
 gTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                                             <  1:2844952/68‑1 (MQ=255)
    ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                                         >  1:2147206/1‑69 (MQ=255)
    ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                                         >  1:908499/1‑69 (MQ=255)
    ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                                         >  1:1947710/1‑69 (MQ=255)
      tCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATAcc                                       <  1:1783625/69‑1 (MQ=255)
        tGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGa                                     <  1:301596/69‑1 (MQ=255)
        tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCg                                      >  1:2233681/1‑68 (MQ=255)
        tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCg                                      >  1:230725/1‑68 (MQ=255)
        tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                     >  1:1663494/1‑69 (MQ=255)
        tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                     >  1:3238365/1‑69 (MQ=255)
        tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                     >  1:62038/1‑69 (MQ=255)
        tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                     >  1:3326747/1‑69 (MQ=255)
        tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                     >  1:3371099/1‑69 (MQ=255)
        tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                     >  1:847866/1‑69 (MQ=255)
        tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                     >  1:2553891/1‑69 (MQ=255)
        tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                     >  1:1335231/1‑69 (MQ=255)
        tGGGGCTATACCCGCATTACACAACGGCTGTCGATAATATGGTTCGCGTGAAAGATATCTTCATACCGa                                     >  1:2431865/1‑69 (MQ=255)
        tGGGGCGATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa                                     >  1:3020/1‑69 (MQ=255)
          gggCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGa                                     >  1:2714913/1‑67 (MQ=255)
               ataCCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTgaga                               >  1:1582180/1‑68 (MQ=255)
                taCCCGGATTACGCAACGGCTGTCGATAACATGGTTCGGGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:2464983/69‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:1512264/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:363335/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:3371270/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:2655540/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:3327540/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:1143153/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:2761953/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:2844417/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:168967/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:2139977/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:388596/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:726139/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:1425078/67‑1 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:3381351/1‑69 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:1360916/1‑69 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:268767/1‑69 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:2826826/1‑69 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:1524093/1‑69 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:1481573/1‑69 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:2086497/1‑69 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:3274752/1‑69 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:221201/1‑69 (MQ=255)
                  cccGGATTACGCAACGGCTGTCGATAACATGGATCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:2132856/1‑69 (MQ=255)
                  ccaggATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa                           >  1:2277287/4‑69 (MQ=255)
                  cacGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc                             <  1:931381/65‑1 (MQ=255)
                    cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                         <  1:2248126/69‑1 (MQ=255)
                    cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                         <  1:1616353/69‑1 (MQ=255)
                    cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                         <  1:1395448/69‑1 (MQ=255)
                    cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                         <  1:434827/69‑1 (MQ=255)
                    cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                         <  1:1833466/69‑1 (MQ=255)
                    cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                         <  1:2583056/69‑1 (MQ=255)
                    cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                         <  1:200735/69‑1 (MQ=255)
                    cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg                         <  1:1271559/69‑1 (MQ=255)
                               aCGGCTGTCGATAAAATGGTTCGCGTGAAATATATCTTCATACCGATTGAGAGCAATGCCAAACGCTac              <  1:1996874/69‑1 (MQ=255)
                               aCGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTac              <  1:267360/69‑1 (MQ=255)
                               aCGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTac              <  1:1354842/69‑1 (MQ=255)
                               aCGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTac              <  1:2418151/69‑1 (MQ=255)
                               aCGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTac              <  1:3087565/69‑1 (MQ=255)
                               aCGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTac              <  1:1481915/69‑1 (MQ=255)
                                 ggCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTacga            <  1:1059939/69‑1 (MQ=255)
                                 ggCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTacga            <  1:2289315/69‑1 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCt                >  1:1707379/1‑61 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGActcg        >  1:2638911/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:3016647/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:1248946/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:869596/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:857445/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:164707/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:636503/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:632220/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:1703275/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:523043/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:429925/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:429752/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:1894122/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:2365126/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:243214/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:2509417/1‑69 (MQ=255)
                                     gTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:2990893/1‑69 (MQ=255)
                                     gTCAATAAAATGGTTCGCGTGAAAGATATCTTCATACCGATTGAGAGCAATGCCAAACGCTACGAcgcg        >  1:122519/1‑69 (MQ=255)
                                           aaCATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAt  <  1:177963/69‑1 (MQ=255)
                                           aaCATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAt  <  1:1358426/69‑1 (MQ=255)
                                           aaCATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAt  <  1:3074233/69‑1 (MQ=255)
                                           aaCATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAt  <  1:1151296/69‑1 (MQ=255)
                                             |                                                                  
GGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATGCCAAACGCTACGACGCGATGAAT  >  minE/1872628‑1872739

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: ‑