Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,126,645:1 | +C | 100% | intergenic (‑99/+50) | 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 | 90.0% | 85.2 / 3.7 | 30 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (3/0); new base C (13/14); total (16/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.28e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTA > NC_002947/1126500‑1126787 | gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAA‑c < 5:335592/147‑1 (MQ=255) gatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTac < 3:361114/149‑1 (MQ=255) ttGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAA‑cccgc > 3:288627/1‑130 (MQ=255) ttGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAA‑cccgc < 4:288627/133‑4 (MQ=255) aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCCAAATG‑GGCGGGTTTGCCACCAACGCAATAGCGGGG‑A‑AAA‑CCCCGCCTACAAACACCCCCCCa > 6:241125/1‑149 (MQ=255) gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTAACCACCAAAGCGGGGGGATAGCCCCCAACGCAATAGC‑GGG‑G‑AAACCCGCGACATCACAGAGCCCCCCCCCCaa > 2:228767/1‑149 (MQ=255) gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTAACCGCGAATGCGGAGGTTTAGCCCCCATCGCATTCGC‑GGGTT‑AAC‑CCGCGGCTACACAGAGCCCGCCCCCCACGcc > 6:393858/1‑149 (MQ=255) aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg < 2:325236/148‑1 (MQ=255) aTGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTTACCGCCAATTCGGAGGTATAGCCACCATAGCATTAGC‑GGG‑TAAAA‑CCGCGCCTACACAGAGCCCCCCACCCAAGCCAGGGGGCGGGGAggg > 2:294266/1‑148 (MQ=255) tgtgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCAccc > 7:382729/1‑127 (MQ=255) tgtgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCAccc < 8:382729/127‑1 (MQ=255) ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTTACCCCGAAAGCGGAGGTATAGCCACCCTCCCCATAGC‑GGG‑T‑AAACCCGCGCCGACACAGAGCCCCCCCCCCAAGCCCGGGGGCGGGGAGGGCTATAAggg > 6:299453/1‑148 (MQ=255) ttgGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCACGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctc < 1:228767/148‑1 (MQ=255) gTGGACATTGAGCCCTGTAGGCGCGGG‑TTTTACCGCCAATGCGGAGGGATAGGCACCATCGCATTAGC‑GGG‑T‑AAACCCGCGCCGACACAGAGCCCCCCCCCCAACCCAGGTGGGGGGGGGGGCGATAAGGCGGGCGGATAGcgatagat > 4:130854/1‑142 (MQ=255) gcgGG‑TTTACCCGCAAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg < 1:294266/149‑1 (MQ=255) gg‑TTTACCCGCGAATGCGATGGTGGCTATACCTCCGCATTCGC‑GGG‑T‑AAACCCGCGCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg > 3:307456/1‑149 (MQ=21) g‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGc < 1:279266/149‑1 (MQ=255) tttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg > 3:377252/1‑149 (MQ=255) ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc > 8:135737/1‑149 (MQ=255) cgcgAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc > 5:405629/1‑149 (MQ=255) aTGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAaca < 6:405629/149‑1 (MQ=255) ggAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt > 6:397732/1‑97 (MQ=255) ggAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt < 5:397732/97‑1 (MQ=255) gAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTACTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTATCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTg > 7:9000/1‑149 (MQ=255) ggTATAGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc < 5:393858/149‑1 (MQ=255) tataGCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGa > 5:110952/1‑149 (MQ=255) gCCACCATCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg < 7:135737/149‑1 (MQ=255) accaTCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGAGCTATCAGGACGCCTGATAGCcctcccc < 6:132898/95‑4 (MQ=25) accaTCGCATTCGC‑GGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGAGCTATCAGGACGCCTGATAGCcctcccc > 5:132898/1‑92 (MQ=25) aacaTCGCATTCGC‑GGG‑T‑AAACCCGCTCCTAAATAGAGCCAGCCACCCAAGCCAGGTGGCTGGGAGGGCTATCAGGCGTCATGATAGCTCTCGATCGACGGGCAGGAGCACAACAGGTTGCGGTCGACGAACACGTTGTCGAcctgtcc > 2:69343/3‑144 (MQ=255) cctccGCATTCGC‑GGG‑T‑AAACCCGCGCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc < 1:252956/144‑1 (MQ=255) cgc‑gGG‑T‑AAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc < 3:140330/148‑1 (MQ=255) gg‑T‑AAACCCGCTCCGACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTa < 5:119645/149‑1 (MQ=255) | GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑T‑AAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTA > NC_002947/1126500‑1126787 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |