Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | AE015451 | 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 | ||
* | AE015451 | 1,126,645 | 1 | . | C | 81.0% | 58.4 / 9.8 | 21 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (0/4); new base C (8/9); total (8/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.31e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01 |
GCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGC > AE015451/1126444‑1126834 | gcAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAgg < 1:220849/237‑1 (MQ=255) gCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTaca < 1:118231/211‑1 (MQ=255) cccTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGc < 1:298336/237‑1 (MQ=255) tCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg < 1:228386/235‑1 (MQ=255) gtgtGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTCCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑c < 1:83323/159‑1 (MQ=255) ggggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTCTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg > 1:171744/1‑185 (MQ=255) ggggCCGATGCCGATCTTGAACAGGGCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg < 2:171743/185‑1 (MQ=255) cgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAgg < 1:277614/194‑1 (MQ=255) gatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAACC‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTa < 1:17528/188‑1 (MQ=255) gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTGACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGCAA‑ACCCGCTCCTACACAGAGCACGCCACCCAAGCCAGGTGGCgg > 1:179397/1‑158 (MQ=255) aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGa > 1:199851/1‑237 (MQ=255) aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACATAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTCTCAGGCGTCCTGAt > 2:319040/1‑149 (MQ=255) ggTTGTTGGGTGGACATTAAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc < 2:27707/162‑1 (MQ=255) aCAGTGAGCGCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCa < 1:341345/237‑1 (MQ=255) ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAacac > 1:220062/1‑160 (MQ=255) ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGc > 2:216225/1‑237 (MQ=255) cccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAGC‑CCGCTCCTACGCAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGt < 1:45999/186‑1 (MQ=255) aTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGCGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGa > 2:268297/1‑127 (MQ=255) aTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTa < 1:216513/68‑1 (MQ=255) aTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTa > 2:216512/1‑68 (MQ=255) tCGCATTCGCGGGTAAGC‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGcac < 1:274943/107‑1 (MQ=255) cGCATTCGCGGGTAAGC‑CCGCTCCTACACAGCGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg < 2:220061/115‑1 (MQ=255) | GCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGC > AE015451/1126444‑1126834 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |