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 | 89.5% | 52.9 / 1.1 | 19 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (2/0); new base C (5/12); total (7/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.23e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGG‑TAAA‑C‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCC > NC_002947/1126529‑1126772 | gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGGCACCATCGCATTAGCGGG‑GAAACC‑CGCTCCTACACAGAGCCCGCCCCCCCAg > 6:69721/1‑147 (MQ=255) gTGACTTGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGT‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCggg < 2:116811/148‑1 (MQ=255) cATGTGGACCTCCTGCCAGGTTGTGGGGGGTACTTTGAGCCCTGTAGGCGCGGGTTTACCCGCTAATGCTGAGGTATAGCCACCATCGCATTCGCGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg < 8:103556/147‑1 (MQ=37) gtgGACCTCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCt < 8:82699/148‑1 (MQ=255) tcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAATGCGGAGGTATATCCACCATCGCATTAGCGGGTTAAC‑C‑CCCGCCGACACACACCCCCCCCCCCAACCCCCGTGGCGGGGCGGGCTATCAGGc > 3:33913/1‑148 (MQ=255) tggTTGTTGGGTGGACATTGAGCCCTTTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGGGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa < 8:83529/147‑1 (MQ=255) aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa < 3:69632/148‑1 (MQ=255) ttgGGGGGACATTGAGCCCTGTCGGCGCGGGTTTCCCCGCGAATGCGTCGGTATAGCCACCATCGCAGTCGCGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGct < 2:113621/147‑1 (MQ=37) ggggccATTGCGCGCTGTAGGCGCGGGGTTCCCCGCGAAGGCGGAGGTATACCCACCACCGCATTCGCGGT‑AAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTAGCAGGCGTCCTGATAGCTCtcga < 5:8497/143‑1 (MQ=255) gTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCGAAAGCGGAGGGATAACCACCATACCAATAGCCGG‑GAAA‑CACCCGCCTAAAAAAAGCCCCCCCCCCAACCCACGGGGCGGGGTGGGCTATAAGGCGGCCAGATaagtctag > 2:30063/1‑140 (MQ=255) cttggcgCCCGGTTGGGGCGGTTTTACCCGCGAATGCGCAGCTATAGCCACCAGCCCATTCGCGGG‑AAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc < 2:64160/142‑1 (MQ=255) acgccTGTAGGCGCGGGTTTAGCCGCGAATGCGGAGGTATAGCCCCCATCGCATTCGCGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACggg < 7:30017/143‑1 (MQ=38) ccTGTAGGCGCGGGTTTACCCCCGAATGCGGAGGGATAGCCACCCTCGCATTAGCGGG‑TAAACC‑CCCTCCTACACAGAGCCCCCCCCCCCACCCACGGGGCGGGGAGGGCTATCAGGCGGCCGGCTAGCGGGCGAACGACGGGCCGgc > 6:68641/1‑148 (MQ=255) taaTGCGGAGGTATAGCCACCATCGCATTCGCGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGa < 2:113758/147‑1 (MQ=255) ggAGGTATAGCCACCATCGCATTCGCGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt > 7:59357/1‑148 (MQ=255) ggAGGTATAGCCACCATCGCATTCGCGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt > 3:18908/1‑148 (MQ=255) ggAGGTATAGCCACCATCGCATTCGCGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt > 2:39484/1‑148 (MQ=255) gAGGTATAGCCACCATCGCATTCGCGGG‑TAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt < 4:33913/148‑1 (MQ=255) ccaCCGCATGCGGGGT‑AAAACC‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc < 5:102745/148‑1 (MQ=255) | GAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGG‑TAAA‑C‑CGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCC > NC_002947/1126529‑1126772 |
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: ‑ |