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 | 95.7% | 72.6 / ‑2.8 | 23 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (9/13); total (10/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.35e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAA‑CC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGC > NC_002947/1126503‑1126779 | gggCCGCTGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGACCACCATCGCATTCGCGGGTAAA‑cccgc < 1:22265/149‑4 (MQ=38) gccgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTcc > 7:91160/1‑149 (MQ=255) gTCGAAGACGCTCAGTGACCTGTGCCCCTCCTGCCAGGTTGTTGGGTGGACATTGCGCCCTGTCGGCGCGGGTTAACCCGCTATTGCGGAGGTATAC‑CCACAATCGCATTCGCGGTTAAACCC‑GCTCCTACACAGAGCCCGCCACCCaa < 8:198101/149‑1 (MQ=255) gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCaa < 2:87583/149‑1 (MQ=255) gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAA‑c > 6:144946/1‑118 (MQ=255) gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAA‑c < 5:144946/118‑1 (MQ=255) gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGACCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGcc < 7:300340/149‑1 (MQ=255) tGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa < 4:34007/149‑1 (MQ=255) gACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATc < 4:157601/148‑1 (MQ=255) aCCTCCTGCCAGGTTGTTGGTTGTACATTGAGCCCTGTAGTCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCTCGTGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg < 4:272384/149‑1 (MQ=38) ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc < 1:128022/149‑1 (MQ=255) ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTg < 8:91160/149‑1 (MQ=255) gttgttGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCaa < 7:300874/112‑1 (MQ=255) ccGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGgcgc > 8:313885/1‑149 (MQ=255) ccGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGAGCTATCAGGACGCCTGATAg > 7:309799/1‑109 (MQ=37) ccGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGAGCTATCAGGACGCCTGATAg < 8:309799/109‑1 (MQ=37) ccGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAA‑CCAGCTCCTACACAGAGCCCGCAACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc > 6:357667/1‑149 (MQ=255) cgcgAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc > 5:40171/1‑149 (MQ=255) gcgAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc > 3:159807/1‑148 (MQ=255) aaTGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac < 7:18876/149‑1 (MQ=255) ggAGGTATAG‑CCACCATCGCATTAGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt > 2:221817/1‑149 (MQ=255) aG‑CCACCATCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc > 4:227540/1‑149 (MQ=255) caccaTCGCATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc < 7:85190/149‑1 (MQ=255) cATTCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc > 4:209390/1‑149 (MQ=255) ttCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc < 4:78383/140‑1 (MQ=255) ttCGCGGGTAAACCC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc > 3:78383/1‑140 (MQ=255) | GGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAG‑CCACCATCGCATTCGCGGGTAAA‑CC‑GCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGC > NC_002947/1126503‑1126779 |
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: ‑ |