Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,555,181 C→T 50.0% R120R (CGC→CGT dsdA → D‑serine ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,555,1810CT50.0% ‑5.2 / 38.8 44R120R (CGC→CGTdsdAD‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (16/6);  new base T (15/7);  total (31/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGTTACTGAAAAAAGATAGCCATTTGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGAAGCGGGGTTGCTGACG  >  minE/1555126‑1555244
                                                       |                                                               
aaCTGTTACTGAAAAAAGATAGCCATTTGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTa                                                      >  1:3296450/1‑67 (MQ=255)
aaCTGTTACTGAAAAAAGATAGCCATTTGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTa                                                      >  1:1629463/1‑67 (MQ=255)
            aaaaaTACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGt                                                >  1:757434/1‑61 (MQ=255)
            aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:1146509/1‑67 (MQ=255)
            aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:846293/1‑67 (MQ=255)
            aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:3312805/1‑67 (MQ=255)
            aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:2680598/1‑67 (MQ=255)
            aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:1856553/1‑67 (MQ=255)
            aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:1700188/1‑67 (MQ=255)
            aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:1131745/1‑67 (MQ=255)
            aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:1526967/1‑67 (MQ=255)
            aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGcc                                         >  1:2861561/1‑67 (MQ=255)
            aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGcc                                         >  1:2705092/1‑67 (MQ=255)
             aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:999877/68‑3 (MQ=255)
             aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:1277331/68‑3 (MQ=255)
             aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:2445274/68‑3 (MQ=255)
             aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:346316/68‑3 (MQ=255)
             aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:3353819/68‑3 (MQ=255)
             aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:1606249/68‑3 (MQ=255)
             aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:3098251/68‑3 (MQ=255)
                aGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTg                                            >  1:2240034/1‑61 (MQ=255)
                aGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGcccatg                                     >  1:1939818/1‑66 (MQ=255)
                aGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGcccatg                                     >  1:1901940/1‑66 (MQ=255)
                aGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGcccatg                                     >  1:3225326/1‑66 (MQ=255)
                           tGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGAtt                                                        >  1:1925736/1‑38 (MQ=255)
                           tGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACt                           >  1:1973960/1‑67 (MQ=255)
                           tGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTgg                         >  1:293777/1‑69 (MQ=255)
                           tGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTgg                         >  1:813951/1‑69 (MQ=255)
                           tGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTgg                         >  1:2421764/1‑69 (MQ=255)
                                   tCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGaa                 >  1:262062/1‑69 (MQ=255)
                                          ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCAGACGCAGaaaa                              >  1:2143537/1‑49 (MQ=255)
                                          ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:2529459/49‑1 (MQ=255)
                                          ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:18646/49‑1 (MQ=255)
                                          ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              >  1:662595/1‑49 (MQ=255)
                                          ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              >  1:366951/1‑49 (MQ=255)
                                          ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:2548866/49‑1 (MQ=255)
                                          ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:2719207/49‑1 (MQ=255)
                                          ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              >  1:2826272/1‑49 (MQ=255)
                                          ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              >  1:2925413/1‑49 (MQ=255)
                                             tAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:1373960/46‑1 (MQ=255)
                                               aaaGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGAAGCGGGGTTGCTg     >  1:552345/1‑69 (MQ=255)
                                               aaaGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGAAGCGGGGTTGCTg     >  1:3016701/1‑69 (MQ=255)
                                                 aGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:91914/42‑1 (MQ=255)
                                                  ggaCGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGAAGCGGGGTTGCTGACg  >  1:1840782/3‑69 (MQ=255)
                                                       |                                                               
AACTGTTACTGAAAAAAGATAGCCATTTGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGAAGCGGGGTTGCTGACG  >  minE/1555126‑1555244

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