Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,612,061 A→G 36.0% D364D (GAT→GAC) ‡ eutA ← reactivating factor for ethanolamine ammonia lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,612,0610AG36.0% 51.9 / 52.1 86D364D (GAT→GAC) ‡eutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base A (24/31);  new base G (10/21);  total (34/52)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGG  >  minE/1611995‑1612124
                                                                  |                                                               
cAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGa                                                                 >  1:1072591/1‑67 (MQ=255)
cAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGa                                                                 >  1:3405068/1‑67 (MQ=255)
cAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGa                                                                 >  1:2727832/1‑67 (MQ=255)
cAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGa                                                                 >  1:2627489/1‑67 (MQ=255)
cAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGa                                                                 >  1:2402260/1‑67 (MQ=255)
cAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGa                                                                 >  1:2090626/1‑67 (MQ=255)
cAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGa                                                                 >  1:1973629/1‑67 (MQ=255)
cAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGa                                                                 >  1:1446484/1‑67 (MQ=255)
cAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGa                                                                 >  1:1328585/1‑67 (MQ=255)
       gcggcgTAACGCACAGGCGGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:2822798/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGGGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:2131727/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:1098722/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:2273625/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:2379495/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:457178/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:454207/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:255518/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:2666691/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:2448026/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:404641/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:39514/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:2858563/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:390895/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:3199334/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:3275275/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:3286778/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:767336/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:720355/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:980380/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:1511481/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:700784/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:605132/1‑69 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:2057277/69‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATc                                                         >  1:2353891/1‑68 (MQ=255)
       gcggcgTAACGCACAGGCAGAGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:3087461/1‑69 (MQ=255)
            gTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTg                                                   <  1:666643/69‑1 (MQ=255)
             tAACGCATAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:2404313/69‑1 (MQ=255)
             tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:9830/69‑1 (MQ=255)
             tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:163095/69‑1 (MQ=255)
             tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:1856232/69‑1 (MQ=255)
             tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:413125/69‑1 (MQ=255)
             tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:2579801/69‑1 (MQ=255)
             tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:1307234/69‑1 (MQ=255)
             tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:2143613/69‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:1520207/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:515934/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2179688/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:3243509/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:474538/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:3083994/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:3049725/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2968649/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:1462152/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:1265533/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:257595/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:1621681/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2809380/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:472156/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2766479/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2425130/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:840527/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2598080/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2587846/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:446669/67‑1 (MQ=255)
               aCGCACAGGCAGCGAGCCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2279345/67‑1 (MQ=255)
                          gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                     <  1:985722/69‑1 (MQ=255)
                          gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                     <  1:3387491/69‑1 (MQ=255)
                          gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                     <  1:2975654/69‑1 (MQ=255)
                          gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                     <  1:2817081/69‑1 (MQ=255)
                          gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                     <  1:2204798/69‑1 (MQ=255)
                          gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                     <  1:2489156/69‑1 (MQ=255)
                           cGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                     <  1:2284128/68‑1 (MQ=255)
                             aGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCTTGTTGc                                   <  1:4499/68‑1 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGcc                                  >  1:2777597/1‑68 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                 >  1:1389468/1‑69 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                 >  1:202669/1‑69 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                 >  1:695394/1‑69 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                 >  1:610392/1‑69 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                 >  1:219918/1‑69 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                 >  1:886534/1‑69 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                 >  1:1370616/1‑69 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                 >  1:1166626/1‑69 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGAATCAGCGCCTGTTGCCa                                 >  1:2949640/1‑69 (MQ=255)
                                     aGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGccccc                          >  1:3151785/1‑69 (MQ=255)
                                     aGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGcccc                           >  1:1386367/1‑68 (MQ=255)
                                            gCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                     <  1:2155431/51‑1 (MQ=255)
                                                         aGTTTTGGGATCAAGATCCAGCTGAATCCGCGCCTGTTGCCAGGCCCCCACCAGAt                   >  1:2575353/1‑56 (MQ=255)
                                                         aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt      >  1:1098557/1‑69 (MQ=255)
                                                         aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt      >  1:491111/1‑69 (MQ=255)
                                                         aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt      >  1:2494811/1‑69 (MQ=255)
                                                         aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt      >  1:2669320/1‑69 (MQ=255)
                                                         aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt      >  1:40910/1‑69 (MQ=255)
                                                         aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAAt      >  1:3398330/1‑69 (MQ=255)
                                                             ttGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCggg  <  1:1389875/69‑1 (MQ=255)
                                                             ttGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCggg  <  1:2420082/69‑1 (MQ=255)
                                                                  |                                                               
CAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGG  >  minE/1611995‑1612124

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: