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 | 89.5% | 113.7 / 7.4 | 38 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (4/0); new base C (11/23); total (15/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.85e-02 | |||||||||||
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. |
TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG > NC_002947/1126498‑1126785 | tGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAA‑c < 8:58888/149‑1 (MQ=255) ggggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAA‑cccgc < 3:304449/149‑4 (MQ=255) gggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCCA‑ATTCGGATGTAAAGCCACCATAGCATTAC‑C‑GGG‑TAAA‑cccgc > 8:115245/1‑145 (MQ=255) cgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTa < 6:91738/149‑1 (MQ=255) gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGACATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGccc < 1:281150/147‑1 (MQ=255) gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGc < 3:62218/148‑1 (MQ=255) gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGC‑ATTGGGCGGGTTATGCCCCATCGCCTTAG‑CGGGG‑TAAC‑CCCCGCCGACACACACCCCCCCCCCCaa > 7:34085/1‑149 (MQ=255) tCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATATCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAg > 2:228254/1‑149 (MQ=255) aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGTGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg < 7:115245/147‑1 (MQ=255) aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGA‑AAGCGGAGGTTTAGCCACCCACGCATTAG‑C‑GGGTTAAC‑CCGCGCCGACAAACAGCCCGCCCCCCCAGCCACGTgggggg > 3:313192/1‑149 (MQ=255) aCATGGGGACCTCCTGCCAGGTTGTTGGGTGCACATTGACCCCTGTCGCCGCGGGTTAACCCGCGA‑ATGCGGAGGTATACCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAgg < 7:186081/149‑1 (MQ=37) cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg < 1:228254/149‑1 (MQ=255) cATGGGGACCCCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATACCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg < 7:42539/149‑1 (MQ=38) aTGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGC‑ATGCGGGGGTTTTGGCACCATCGCATTAG‑C‑GGGTTAAA‑CCGCGCCGACACAGAGCCCCCCCCCCCAACCCGGTGGGGGGGGGGGc > 7:262402/1‑149 (MQ=255) aTGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTAG‑C‑GGG‑TAAACCCGCGCCGACACAGAGCCCCCCCCCCAAGCCAGGTGGCGGGGAGGGc > 2:64994/1‑149 (MQ=38) tgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTAt < 2:42086/149‑1 (MQ=255) gACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGA‑ATGCGGGGGTTTAGGCACCAAAGCATTAGCC‑GGT‑TAAC‑CCGCGCCGAAACAGAGCCCCCCACCCAAGACACGTGGCGGGGGGGGCTAt > 2:282268/1‑147 (MQ=255) cAGGTTGTGGGGTGGAAATGGAGCCCTGTAGGCGGGGGTTTACCCGCTA‑ATGCGGAGGTATAGCCACCATCGCATTGG‑C‑GGT‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa < 1:308326/149‑1 (MQ=255) ttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct < 8:34085/149‑1 (MQ=255) aTTGCGCCCTGTAGGCGCGGGTTAACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg < 5:22389/149‑1 (MQ=255) cgtGTTGGCGGGGGTTACCCCGCTA‑TGGGGGTGGAATCCCCACCAGCGCATGCG‑G‑GGT‑AAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcg < 5:52408/147‑1 (MQ=255) tAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGcac < 1:282268/148‑1 (MQ=255) gcgA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGT‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCg < 5:228200/149‑1 (MQ=255) gcgA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc > 1:46792/1‑148 (MQ=255) ggAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt > 8:61980/1‑149 (MQ=255) ggAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt > 6:264815/1‑149 (MQ=255) aGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt < 4:313192/149‑1 (MQ=255) ggTATAGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc < 2:46792/149‑1 (MQ=255) taGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCAAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc > 4:24367/1‑149 (MQ=255) aGCCACCATCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc > 1:51488/1‑149 (MQ=255) gCCACCATCGCATTAG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcg > 2:313414/1‑111 (MQ=255) gCCACCATCGCATTAG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcg < 1:313414/111‑1 (MQ=255) accaTCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc > 8:214544/1‑149 (MQ=255) caTCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc > 7:73360/1‑123 (MQ=255) caTCGCATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc < 8:73360/123‑1 (MQ=255) cATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc < 2:51488/149‑1 (MQ=255) cATTCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc < 1:64994/149‑1 (MQ=255) ttCG‑C‑GGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcgg < 6:197191/149‑1 (MQ=255) c‑gGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc < 6:295003/146‑1 (MQ=255) c‑gGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc > 5:295003/1‑146 (MQ=255) c‑gGG‑TAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg < 2:217922/149‑1 (MQ=255) | TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGA‑ATGCGGAGGTATAGCCACCATCGCATTCG‑C‑GGG‑TAAA‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG > NC_002947/1126498‑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: ‑ |