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 | 93.0% | 131.5 / 3.5 | 43 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (3/0); new base C (17/23); total (20/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.24e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA > NC_002947/1126504‑1126784 | ggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGGTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAA‑cccgctc < 5:337776/149‑6 (MQ=255) cgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTa < 4:355954/149‑1 (MQ=255) atgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTaca < 8:105931/149‑1 (MQ=255) cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTTCCCGCGAAGGCGGGGGGATAGGCACCCGCGCCTTCGCGGGG‑TAAC‑CCGCGCCGGcaca > 6:297034/1‑147 (MQ=255) aaGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCCG‑ATGCGGGGGTATAGGCACCCAAGCCATAGC‑GGGTTAAA‑CCCCGCCGAAACAGAGccccccccccccccca > 6:231222/1‑143 (MQ=255) aaGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTTACCCCGAATTGGGAGGTATTGCCACCATCGCCTTAGC‑GGGTTAAC‑CCCCGCCGACACACAGCCCCCCCCCCAACCCa > 8:245173/1‑149 (MQ=255) aGTGACATGTGGACCTCCTGCCAGTTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGTATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCggg < 1:313236/149‑1 (MQ=255) gACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTAAACGCGAAAGCGGAGGTATAGCCACCAAAGCATTCGC‑CGG‑TAAACCCCCACCGACACAGAGCCCCCCACCCCACCCCGGTGGCGGGGAGGGCTATc > 2:20900/1‑148 (MQ=255) cctcctGCCAGGTGGTTGGGTGCACATTCAGCCCTGTAGGCGCGGG‑TTAACCCGCGAATGCGGAGGTATACCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAgg < 7:147208/149‑1 (MQ=38) ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTAGC‑GGG‑TAAACCCGCTCCTAAACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc > 2:313486/1‑149 (MQ=255) cctGCCAGGTTGTGGGGGGGACATTCAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGt < 7:332185/149‑1 (MQ=38) ctGCCTGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTc < 1:51154/149‑1 (MQ=255) gCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc < 1:118248/148‑1 (MQ=255) cccAGGTTGTTGGGGGCACATTGAGCCCTGTCGGCGCGGG‑TTACCCCGCTAATGCGGAGGTATCGCCACCATCGCATTCGG‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCt < 7:339802/148‑1 (MQ=255) ggTTGGTGGGTCCACATTCCCCCCTCTAGGCGCGGT‑TTTCCCCGCTAATCCGGCGTAAAACCCACCACCGCATTCGC‑GGT‑AAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATa < 3:205042/149‑1 (MQ=255) gAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGTAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGc < 5:231222/148‑1 (MQ=255) ccTGTCGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg < 3:330614/142‑1 (MQ=255) cTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcg < 7:123264/148‑1 (MQ=255) aGGAGCGGG‑TTTACCCGCGAATGCGATGGTGGCTATACCTCCGCATTCGC‑GGG‑TAAACCCGCGCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCacac > 5:331974/1‑149 (MQ=18) tttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGT‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg < 7:245173/149‑1 (MQ=255) ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc > 7:66867/1‑149 (MQ=255) ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc < 1:20900/149‑1 (MQ=255) cgcgAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc < 8:66867/149‑1 (MQ=255) cgcgAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc > 5:216777/1‑149 (MQ=255) gcgAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGCGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCCGGCGCCCCCCAGGTTGCGGTCGcc > 8:140968/1‑148 (MQ=39) ggAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCa < 5:413956/116‑1 (MQ=255) ggAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCa > 6:413956/1‑116 (MQ=255) ggAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg < 1:260397/132‑1 (MQ=255) ggAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg > 2:260397/1‑132 (MQ=255) ggTATAGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc > 7:335638/1‑149 (MQ=255) aGCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc > 7:21538/1‑149 (MQ=255) gCCACCATCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc < 6:216777/148‑1 (MQ=255) ccaccaTCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg > 2:348149/1‑149 (MQ=255) ccaccaTCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg > 3:270736/1‑149 (MQ=255) caccaTCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg > 6:346812/1‑69 (MQ=255) caccaTCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg < 5:346812/69‑1 (MQ=255) caTCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGATCCCGCCACCCAAGCCAGGTGGCCGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac > 4:188063/1‑130 (MQ=255) caTCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac < 3:188063/130‑1 (MQ=255) aTCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAc < 4:270736/149‑1 (MQ=255) tCGCATTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAcc > 4:238087/1‑149 (MQ=255) aTTCGC‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg > 6:28937/1‑149 (MQ=255) gc‑gGG‑TAAACCCGCGCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCa < 6:331974/149‑1 (MQ=255) g‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc < 2:203717/140‑1 (MQ=255) g‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc > 1:203717/1‑140 (MQ=255) | GGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGG‑TTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGC‑GGG‑TAAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA > NC_002947/1126504‑1126784 |
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: ‑ |