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 | 84.6% | 69.9 / 9.1 | 26 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (3/1); new base C (10/12); total (13/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.93e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGT > NC_002947/1126503‑1126786 | gggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGCCCTCCTGCCAGGTTGTTGGGTGGCCATTGAGCCCTGTAGGCGCGGGGTTACCCGCCCATGCGGACGTATAGCCACAAAAG‑CCTTCGC‑CGG‑GA‑AA‑cccccg > 7:135495/1‑145 (MQ=255) cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCGCGAATGCGGCGGGTTAGCCCCCATCG‑CCTTCGC‑GGGTTA‑AC‑CCGCTCCGAcaca > 4:187324/1‑147 (MQ=255) atCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAAAGCGGAGGGATATCCACCCTCGCCATCCCC‑GGT‑AA‑AA‑CCGCGACAAAACAGAGc > 8:243363/1‑149 (MQ=255) aaCAGGTCGAAGATGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAATGCGGGGGGATAGGCACCATCG‑CCTTAGC‑GGG‑TA‑AACCCGCGCCTACCCAGCGccccccc > 4:101233/1‑148 (MQ=255) tCGAAGTCGCTGAGTGACTTGTGTACCTAGTGCCATGTTGTTGGGTTTACATTGAGCCCTTTAGGGGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCaa < 3:194489/148‑1 (MQ=25) tCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACACGCCAATTAAGAGGAATAGACACCATCC‑CATTACC‑CCG‑TA‑AACCCCCTCCAACAAAAACCCCCCCACCCaa > 4:194489/1‑148 (MQ=255) aaGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGc < 5:280408/138‑1 (MQ=255) aGACGCTCAGGGACTTGTGGACCTCCTGCCCGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAAGGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAg < 8:135495/149‑1 (MQ=38) agggccATGGGGACCTCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCggg < 3:187324/144‑1 (MQ=38) tGACATGTGGACCTCCTCCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTCCCCGCGAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa < 1:295389/149‑1 (MQ=38) tcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAAACCCGAATTCGGAGGTATAGCCACCATAG‑CATTAGCGGGG‑TA‑AC‑CCGCTCCTAAAAACAGCCCGCCACCCAACCCCCGGGGCGGGGAGGGCTATCAGGCg > 7:178688/1‑149 (MQ=255) cctGCCAGGTGGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCGCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGt < 7:243363/149‑1 (MQ=255) tGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGATGGCTATCAGGCGTcc > 3:79051/1‑149 (MQ=255) cccTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCAGCC‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt < 7:210540/127‑1 (MQ=255) cTGTAGGCGCGGGTTAACCCGCGAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTc < 7:62915/120‑1 (MQ=255) tttACCCGCGAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg < 8:178688/149‑1 (MQ=255) gcgAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACAcc > 5:258941/1‑134 (MQ=255) gcgAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACAcc < 6:258941/134‑1 (MQ=255) ggAGGTATAGCCACAATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATAAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCACGTTGCGGTAGCCGAACTCg > 4:62124/1‑147 (MQ=255) ataGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc < 4:259484/83‑1 (MQ=255) ataGCCACCATCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc > 3:259484/1‑83 (MQ=255) ccaccaTCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTg > 8:75712/1‑85 (MQ=255) ccaccaTCG‑CATTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTg < 7:75712/85‑1 (MQ=255) aTTCGC‑GGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCa > 6:108409/1‑96 (MQ=255) cgc‑gGG‑TATAC‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGATGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc < 4:79051/148‑1 (MQ=255) gc‑gGG‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc > 7:97970/1‑147 (MQ=255) gg‑TA‑AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGt > 3:177705/1‑148 (MQ=255) | GGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCG‑CATTCGC‑GGG‑TA‑AA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGT > NC_002947/1126503‑1126786 |
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: ‑ |