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 | 90.5% | 129.6 / 7.1 | 42 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (4/0); new base C (16/22); total (20/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.33e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAA‑C‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG > NC_002947/1126497‑1126785 | cTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑aaa < 3:177295/149‑1 (MQ=255) ggAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGG‑AGGGATAGCCACCATAG‑CAATCGC‑CGGT‑AAA‑c‑c > 6:116207/1‑149 (MQ=255) ccgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTcc > 1:41732/1‑148 (MQ=255) aGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAAATCGGTAGGT‑TAGCCCCCAAAG‑CCATAGCGGGGT‑AAA‑C‑CGCGGCGAACCACAGCCCGcccccc > 2:142187/1‑149 (MQ=255) gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGT‑AGTTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGc < 4:132194/148‑1 (MQ=255) gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGc < 5:116207/148‑1 (MQ=255) cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg < 2:41732/149‑1 (MQ=255) cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGATTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg < 5:341481/149‑1 (MQ=255) aCCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg < 5:302058/149‑1 (MQ=255) ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAATGCGG‑GGGTATAGCCACCAAAGCCATAGCC‑GGGT‑AAA‑C‑CCCAACAAAACACAGCCCCCCCCCCCACCCCCGGGGCGGGGGGGGCGATAAGGc > 4:348485/1‑149 (MQ=255) ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCCAATTAGG‑AGGTATATCCACCAAAC‑CAATAGC‑CGGTAAAA‑C‑CCCTCCAACAAAGACCCCCCCCCCCCAGCCACGTGGCGGGGGGGGCTATCATGc > 8:177917/1‑149 (MQ=255) ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc < 3:348485/149‑1 (MQ=255) ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc < 7:177917/149‑1 (MQ=255) ctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCTCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTc < 8:177944/149‑1 (MQ=255) ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGTTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc < 1:142187/147‑1 (MQ=255) cAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCTGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa < 5:347170/149‑1 (MQ=255) gttgttGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGACGCCTGATAg < 1:183176/149‑1 (MQ=38) gttgGGTGGACATTCCGCCCTGTAGGCGCGGGTTTACCCGCTAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGct < 1:179941/148‑1 (MQ=38) gtgGGTGGACATTGAGCCCTGTCGGCGCGGGTTTACCCGCGAATGCGG‑CGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct < 5:181434/148‑1 (MQ=255) tgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct < 8:16430/148‑1 (MQ=255) ggggTGGAATTTTAGCCCGTTAGGCGCGGTTTTACCCGCTAAGGCGT‑CTGTATAGCCACCATCG‑AATTCGC‑GGTT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGAGCTATCAGGACGCCTGATAGccct < 5:284370/147‑1 (MQ=255) gggTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc < 7:330498/148‑1 (MQ=255) gTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt > 7:202099/1‑113 (MQ=255) gTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt < 8:202099/113‑1 (MQ=255) tttACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAA‑CACGCTCCTACACAGAGCAAGCAAAACAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAAGGGCAGGCGCACACCAGGTTGCg > 5:208577/1‑149 (MQ=38) aCCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc > 3:26294/1‑149 (MQ=255) cgcgAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc > 4:54616/1‑98 (MQ=255) cgcgAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc < 3:54616/98‑1 (MQ=255) gcgAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg < 3:345572/95‑1 (MQ=255) gcgAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg > 4:345572/1‑95 (MQ=255) gg‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt > 5:10955/1‑149 (MQ=255) gg‑AGGTATAGCCACCATCG‑CAATCCC‑GGGT‑AAACC‑CGCTCCTCCACAGAGCCCGCCCCCCAACCCAGCTGGCGGGGAGGGCTATCCGGCGTCCAGATAGCTCTCGATCGACGGGCCGGCGCACACCACGATGCGCTCGCCGAACACgat > 5:410822/1‑147 (MQ=255) ataGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc > 2:361229/1‑149 (MQ=255) taGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc > 1:228900/1‑149 (MQ=255) taGCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc > 8:360966/1‑149 (MQ=255) gCCACCATCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg > 5:181574/1‑149 (MQ=255) ccaccaTCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg > 3:2015/1‑149 (MQ=255) accaTCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac > 4:15607/1‑132 (MQ=255) accaTCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac < 3:15607/132‑1 (MQ=255) ccaTCG‑CATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc > 8:29254/1‑148 (MQ=255) cATTCGC‑GGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc > 7:22014/1‑149 (MQ=255) cgc‑gGGT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc < 4:26294/148‑1 (MQ=255) gc‑gGGA‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCa < 1:203271/149‑1 (MQ=255) ggT‑AAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg > 4:320047/1‑147 (MQ=255) | CTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGG‑AGGTATAGCCACCATCG‑CATTCGC‑GGGT‑AAA‑C‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG > NC_002947/1126497‑1126785 |
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: ‑ |