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 | 97.3% | 123.6 / ‑3.3 | 37 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (17/19); total (18/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGA > NC_002947/1126498‑1126774 | tGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAA‑c < 1:94885/149‑1 (MQ=255) gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAA‑c < 1:108413/147‑1 (MQ=255) ggggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAAA‑GCGGATTACGATGTATAGCAGCCATCGCATTCGC‑GGGTAAA‑cccgg > 8:347132/1‑146 (MQ=255) ggggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAA‑cccgc < 8:73881/149‑4 (MQ=255) gggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCTGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAA‑cccgc < 2:287352/148‑4 (MQ=255) cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTAcaca > 5:383109/1‑147 (MQ=255) cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACC‑GCGAATGGGGGGGTATAGCCACCATAGCATTCGC‑GGGTAAACCCGCGCCGAcaca > 6:198653/1‑147 (MQ=255) gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACagag < 6:383109/143‑1 (MQ=255) aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGGGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTa > 5:202034/1‑135 (MQ=255) aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTa < 6:202034/135‑1 (MQ=255) aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACC‑GCGAATGCGGAGGTATAGCCACCATCGCATTAGC‑GGGTAAACCCGCGCCTACACAGAGCCCGccc > 5:94172/1‑148 (MQ=255) gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑CCGAATTCCGAGGTATAGCCACCATACCAATAGC‑CGGTAAACCCCCACCTACACAAAGCCCGccccccca > 2:303497/1‑147 (MQ=255) tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACC‑GCGAATGGGGGGGGTTAGGCACCAACGCATTAGCGGGGTAAC‑CCCCACCGACACAAAGCCCCCCCCCCCAACCACGTgggg > 2:338257/1‑147 (MQ=255) aTGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGc < 6:181599/149‑1 (MQ=255) tgGCCCTCCTGCCTGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGTGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTAt < 5:338016/149‑1 (MQ=38) tcctGCCAGGTTGTTGGGGGGACTTTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTg < 7:216204/129‑1 (MQ=255) ccAGGTTGTTGGGTGGACATTGAGCCCTGTCGGCGCGGGGTTACCC‑GCGAATGCGGCGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTcc < 1:368411/147‑1 (MQ=38) ggTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑gg < 3:255686/125‑1 (MQ=255) aTTGAGCCCTGTCGGCG‑GGGCTTACCCGGCCAATGCGGCGGTATACCCACCAACGAATTCGC‑GGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACgg < 3:10372/149‑1 (MQ=255) aTTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACgg < 5:198653/149‑1 (MQ=255) aTTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACgg < 3:153358/149‑1 (MQ=255) ggCGGGGTTTTCCCC‑GCTAATGCGTAGGAAAACCCCCCATCCCATGCGG‑GGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACAcc < 5:252635/149‑1 (MQ=255) gcgcGGTTTTACCC‑GCGAATGCTGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATATCT‑CTCGATCGACGTGCAGGCGCACAcc < 7:70980/148‑1 (MQ=38) gcgGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACACCAgg > 4:81341/1‑149 (MQ=255) ttACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTc > 5:72942/1‑106 (MQ=255) ttACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTc < 6:72942/106‑1 (MQ=255) cc‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGcaca < 1:146930/134‑1 (MQ=255) cc‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGcaca > 2:146930/1‑134 (MQ=255) gcgAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGCGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCATGCTAGCT‑CTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc > 8:297014/1‑148 (MQ=255) gcgAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGgcgc > 2:188413/1‑129 (MQ=255) gcgAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGgcgc < 1:188413/129‑1 (MQ=255) ggAGGTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt > 7:334555/1‑149 (MQ=255) ggTATAGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGGTCCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt < 1:303497/149‑1 (MQ=39) ataGCCACCATCGCGTTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc > 7:342881/1‑149 (MQ=255) taGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAgg > 8:70589/1‑111 (MQ=255) taGCCACCATCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAgg < 7:70589/111‑1 (MQ=255) ccaccaTCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGTGGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc > 1:121017/1‑149 (MQ=255) accaTCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc > 8:195732/1‑149 (MQ=255) caTCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGCCGGGCAGGCGCACACCAGGTTGCGGTCGCCGCACACGTTGTCGACCCGGCCGa > 2:373611/1‑149 (MQ=255) caTCGCATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGa < 7:195732/149‑1 (MQ=255) cATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt > 3:15595/1‑130 (MQ=255) cATTCGC‑GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt < 4:15595/130‑1 (MQ=255) | TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCC‑GCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGGTAAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGC‑‑GGGGAGGGCTATCAGGCGTCCTGATAGCT‑CTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGA > NC_002947/1126498‑1126774 |
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: ‑ |