Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,612,062 T→C 36.0% D364G (GAT→GGT) ‡ 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,0620TC36.0% 51.3 / 51.8 86D364G (GAT→GGT) ‡eutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base T (24/31);  new base C (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 = 4.12e-01

GCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGG  >  minE/1612002‑1612124
                                                            |                                                              
gcggcgTAACGCACAGGCGGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:2822798/1‑69 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGGGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:2131727/69‑1 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:2666691/69‑1 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:605132/1‑69 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:2273625/1‑69 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:2379495/69‑1 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:2448026/1‑69 (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:1098722/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:3286778/1‑69 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:3199334/1‑69 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:3275275/1‑69 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:720355/69‑1 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:767336/69‑1 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:2057277/69‑1 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:980380/69‑1 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        >  1:700784/1‑69 (MQ=255)
gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATcc                                                        <  1:1511481/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:1856232/69‑1 (MQ=255)
      tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:163095/69‑1 (MQ=255)
      tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:1307234/69‑1 (MQ=255)
      tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:9830/69‑1 (MQ=255)
      tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:2579801/69‑1 (MQ=255)
      tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:2143613/69‑1 (MQ=255)
      tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGa                                                  <  1:413125/69‑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:1265533/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:840527/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:3243509/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2809380/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2587846/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2766479/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:1462152/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2598080/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:257595/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:446669/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:1520207/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:1621681/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2425130/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:472156/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:474538/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:515934/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2179688/67‑1 (MQ=255)
        aCGCACAGGCAGCGAGCCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGGCAAAGATCCAGTTGa                                                  <  1:2279345/67‑1 (MQ=255)
                   gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                     <  1:3387491/69‑1 (MQ=255)
                   gCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGtt                                     <  1:985722/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: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:1389468/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)
                                                            |                                                              
GCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGG  >  minE/1612002‑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: