Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 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 | ||
* | AE015451 | 1,126,645 | 1 | . | C | 90.5% | 63.9 / 1.6 | 21 | intergenic (‑99/+50) | tdcG‑II/gcvP‑I | L‑serine dehydratase/glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (9/10); total (9/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGT > AE015451/1126411‑1126847 | cGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑c < 2:215145/236‑1 (MQ=255) ccGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCaa < 1:126675/235‑1 (MQ=255) ttCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑cccgctc > 1:355476/1‑197 (MQ=255) cGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg < 1:321232/234‑1 (MQ=255) aGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg > 1:390727/1‑221 (MQ=255) gtgcATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCAAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGCTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTaca < 1:444373/184‑1 (MQ=255) ggAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg < 2:335559/194‑1 (MQ=255) gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGc > 1:491647/1‑165 (MQ=255) aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑cccg < 2:274727/111‑3 (MQ=255) aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑cccg > 1:274731/1‑109 (MQ=255) ggACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatc > 2:442519/1‑172 (MQ=255) cctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGATTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg > 1:346143/1‑235 (MQ=255) aTTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGCGGGCTATCGGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg > 1:382265/1‑168 (MQ=255) cccTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg < 2:311205/104‑1 (MQ=255) cccTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAgg < 1:205493/118‑1 (MQ=255) gcgcGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCt < 1:269405/236‑1 (MQ=255) cGGAGGTATAGCCACCATCGGATTCGCGGGTAACC‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgg < 2:397042/75‑1 (MQ=255) caccaTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATgggc > 1:323115/1‑224 (MQ=255) ccaTCGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTTGCGGCTGTAg > 1:126242/1‑215 (MQ=255) cGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgat < 2:240563/90‑1 (MQ=255) cGCATTCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgat > 1:240567/1‑90 (MQ=255) ttCGCGGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTc < 2:390723/174‑1 (MQ=255) cgtgggTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGc < 1:312765/171‑1 (MQ=255) gcgGGTAAGC‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc < 1:143432/80‑1 (MQ=255) cgGGTAA‑ACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTc > 1:321865/1‑170 (MQ=255) | CGCTGGCGGTGCGGGCCAGCAGGCCGTCACGGCGCAGCCCTTCGGCAAACCGCGCAGCGGCGCGCATCGGGCCGACCGTGTGGGAGCTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAA‑A‑CCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCAGGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGT > AE015451/1126411‑1126847 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |