Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,872,656 G→A 19.2% A427T (GCA→ACA)  ygcE → predicted kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,872,6560GA19.2% 143.1 / 21.3 93A427T (GCA→ACA) ygcEpredicted kinase
Reads supporting (aligned to +/- strand):  ref base G (21/54);  new base A (12/6);  total (33/60)
Fisher's exact test for biased strand distribution p-value = 4.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATG  >  minE/1872590‑1872716
                                                                  |                                                            
aaCGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCa                                                            >  1:347793/1‑69 (MQ=255)
              tcTGGGGGCTGCGATCAATACGGCGGTAGGTCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:2899630/68‑1 (MQ=255)
              tcTGGGGGCTGCGATCAATACGGCGGTAGGTCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:1199489/68‑1 (MQ=255)
              tcTGGGGGCTGCGATCAATACGGCGGTAGGTCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:2001243/68‑1 (MQ=255)
              tcTGGGGGCTGCGATCAATACGGCGGTAGGTCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:1900192/68‑1 (MQ=255)
              tcTGGGGGCTGCGATCAATACGGCGGTAGGGCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:2063511/68‑1 (MQ=255)
              tcTGGGGGCTGCGACCAATACGGCGGTAGGTCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:3303850/68‑1 (MQ=255)
              ccTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATaa                                              <  1:2189897/69‑1 (MQ=255)
                     gCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGAt                                                <  1:111536/60‑1 (MQ=255)
                     gCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGAt                                                <  1:2413390/60‑1 (MQ=255)
                     gCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGAt                                                <  1:1637856/60‑1 (MQ=255)
                     gCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGAt                                                <  1:1410913/60‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAAGATGGTTCGCGTGAAAGATATCtt                      <  1:2001960/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:926472/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:271771/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2708246/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1014479/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2377137/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2920826/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2260819/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1180047/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1254766/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2026466/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1981500/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1402168/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:356783/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1634574/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1569386/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1667164/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1519356/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:470986/69‑1 (MQ=255)
                                      ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGATCGCGTGAAAGATATCtt                      <  1:2061144/69‑1 (MQ=255)
                                       gTAGGTCTGTGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2347966/68‑1 (MQ=255)
                                       gTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2844952/68‑1 (MQ=255)
                                          ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                  >  1:1947710/1‑69 (MQ=255)
                                          ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                  >  1:908499/1‑69 (MQ=255)
                                          ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                  >  1:2231693/1‑69 (MQ=255)
                                          ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                  >  1:2147206/1‑69 (MQ=255)
                                          ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                  >  1:2515628/1‑69 (MQ=255)
                                            tCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATAcc                <  1:1783625/69‑1 (MQ=255)
                                              tGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGa              <  1:301596/69‑1 (MQ=255)
                                              tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCg               >  1:230725/1‑68 (MQ=255)
                                              tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCg               >  1:2233681/1‑68 (MQ=255)
                                              tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:3326747/1‑69 (MQ=255)
                                              tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:3238365/1‑69 (MQ=255)
                                              tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:2553891/1‑69 (MQ=255)
                                              tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:3371099/1‑69 (MQ=255)
                                              tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:1663494/1‑69 (MQ=255)
                                              tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:1335231/1‑69 (MQ=255)
                                              tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:62038/1‑69 (MQ=255)
                                              tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:847866/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:388596/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:363335/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:984664/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:3327540/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:2844417/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:3371270/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:726139/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:1143153/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:2655540/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:1425078/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:1512264/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:2139977/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:168967/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:2761953/67‑1 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:3381351/1‑69 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:1360916/1‑69 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:1481573/1‑69 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:1524093/1‑69 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:1675949/1‑69 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:2826826/1‑69 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:2086497/1‑69 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:3274752/1‑69 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:268767/1‑69 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:221201/1‑69 (MQ=255)
                                                        cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:205960/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:434827/69‑1 (MQ=255)
                                                          cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:1271559/69‑1 (MQ=255)
                                                          cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:1395448/69‑1 (MQ=255)
                                                          cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:2583056/69‑1 (MQ=255)
                                                          cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:2248126/69‑1 (MQ=255)
                                                          cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:1616353/69‑1 (MQ=255)
                                                          cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:349137/69‑1 (MQ=255)
                                                          cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:1833466/69‑1 (MQ=255)
                                                          cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:200735/69‑1 (MQ=255)
                                                                  |                                                            
AACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATG  >  minE/1872590‑1872716

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: