Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 1,126,647 | +C | intergenic (‑101/+48) | tdcG‑II ← / ← gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,126,645 | 1 | . | C | 84.0% | 68.1 / 8.4 | 25 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (4/0); new base C (5/16); total (9/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.96e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAA‑CCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGG > NC_002947/1126500‑1126781 | gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTTACCGCCAATGGGGACGTTTAGCCCCCATCGCCTTAGC‑GGGT‑aaaac > 4:110642/1‑146 (MQ=255) atgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTTAACGCGAATTCGGTGGGATATGCACCATCGCATTAGC‑CGGTAAAA‑CCGCTCCTac > 7:131238/1‑148 (MQ=255) ccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCCGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGTTTTACCCGCTAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGTT‑AAACCCGCTCCTAcaca < 6:96690/148‑1 (MQ=255) aaCCGGGCGCAGACGCTCTGGGACATGTGGACGTCCTGCCAGGTTGTGGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATCGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACaga < 8:110386/141‑1 (MQ=25) aaCAGGGCGAAGACGCTCAGGGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACaga < 7:74399/141‑1 (MQ=255) gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTTAACCCGCAAGGGGAGGTATATGCCCCAAAGCAATAGC‑CGGTAAAA‑CCCCGGCTACACAACGCCCCC‑CCCcccacccacggg > 4:127285/1‑144 (MQ=255) tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTTACCC‑CCAATGCGGGGGGTTAGCCCACATCGCCCTAGCGGGGT‑AAC‑CCCCGGCGACACCCAGCCCCC‑CCCCCCACCCACGTGGc > 5:89783/1‑148 (MQ=255) tCAGTGACATGTGGACCTCCTGCCAGGTTGTGGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGc < 5:39766/148‑1 (MQ=255) tcctGCGAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTCCCCGCGAATGCGGCGGTATAGCCACCATCGCATTCGG‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc < 2:102024/148‑1 (MQ=38) gCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTTACCGCGAATTAGGTGGTATAGCCCCCATAGCATTAGC‑CGGTAAAA‑CCGCGACAAAACAAAGCCCCCACCCCCAAG‑CACGTGGCGGGGTGGGGGATCAGGCGt > 4:82830/1‑146 (MQ=255) cagggggttgGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTg < 3:127285/142‑1 (MQ=255) gttgttGGGTGGACATTGAGCCCTGTCGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATa < 3:110642/148‑1 (MQ=255) gttgtgGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAAGGCGGAGGTATAGCCACCATCGCATTCGC‑GGTT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGCGGAGGGCTATCAGGCGTCCTGATa < 6:139316/148‑1 (MQ=39) tgGGGGCACATTGAGCCCGGTAGGCGC‑GGGTTTACCCGCGAAGGCGGAGGTATCGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGGGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct < 7:55380/148‑1 (MQ=37) catgccgCCCTGACGGCGG‑GGTTTTACCCGCTAAGGCGGAGGAATCCCCACCATCCCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc < 5:8321/142‑1 (MQ=255) cATTGCGCCCTGTCGGCGC‑GGGTTTACCCGCGCAGGCGGAGGTATCCCCACCACCGCAGTCGC‑GGTA‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCGCTCGATCGAc < 8:62825/148‑1 (MQ=255) aTTGAGCCCTGTAGGCGA‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GTGT‑AAACCCGCTCCTAAACAGAGCCCGC‑CACCAAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACg < 8:29988/148‑1 (MQ=255) ttccccGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg < 5:122028/145‑1 (MQ=255) ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg < 4:157611/148‑1 (MQ=255) ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGCTCGACGGGCAGGCGCACACCAGGTTGCGGTCg > 1:94329/1‑148 (MQ=255) ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGCCGGGCAGGCGCACACCAGGTTGCGGTCg > 6:78442/1‑148 (MQ=255) ggAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt > 7:154926/1‑148 (MQ=255) gTATAGCCACCATCGCATTCGC‑GGTT‑AAACCCGCTCCTACACAGAGCCGGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACTCCAGGTTGCGGTCGGCGAACACGTTGt < 8:131238/147‑1 (MQ=255) cgCCACCATCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc < 2:115337/147‑1 (MQ=255) caTCGCATTCGC‑GGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCg > 1:6188/1‑148 (MQ=255) gc‑gGGT‑AAACCCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcgg > 2:178817/1‑146 (MQ=255) | GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGT‑AAA‑CCGCTCCTACACAGAGCCCGC‑CACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGG > NC_002947/1126500‑1126781 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |