Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913_3_hsa_tpiA | 389,919:1 | +T | 100% | intergenic (‑192/‑332) | hemB ← / → yaiT | 5‑aminolevulinate dehydratase (porphobilinogen synthase)/pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_hsa_tpiA | 389,919 | 1 | . | T | 83.3% | 87.6 / 13.8 | 36 | intergenic (‑192/‑332) | hemB/yaiT | 5‑aminolevulinate dehydratase (porphobilinogen synthase)/pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein |
Reads supporting (aligned to +/- strand): ref base . (2/4); new base T (16/14); total (18/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
GTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑A‑T‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTT > NC_000913_3_hsa_tpiA/389778‑390055 | gTAGGGATGTTTTATGAAAGTTGTCGCGATGTTTACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAAGAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑A‑t‑t‑ttttt < 2:361346/149‑2 (MQ=255) ttttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTgg > 2:242151/1‑149 (MQ=255) ttttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGCTATT‑ATT‑T‑TTTT‑GGGGTTTgt > 2:10265/1‑148 (MQ=255) cGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTgcg > 1:103646/1‑149 (MQ=255) aaaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGTTTGCgct > 2:13011/1‑149 (MQ=255) gACAAGAAGAGAATGGAAGAGAGGCTAGGACAGCTTTTTAATATAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGaatcac < 2:173356/149‑4 (MQ=255) agaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTaaa > 2:517736/1‑149 (MQ=255) gagaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACtt > 1:63292/1‑149 (MQ=255) aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTttgt > 2:236034/1‑149 (MQ=255) gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑A‑T‑T‑TTTTCGGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCa < 1:581908/149‑1 (MQ=255) ccAGGACATATATTTAATAGAAAAATCAATAATTTTAGATTCACTGAAAATATTTTTAATCTTCATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCaccc < 2:611149/149‑4 (MQ=255) aGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGa > 1:440644/1‑149 (MQ=255) aGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGa > 1:97803/1‑149 (MQ=255) ttAATAGTAACAGCAATAATTATATATGCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑TGTTATT‑ATT‑T‑TTTT‑GGAGTTTGCTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTgg < 2:309983/149‑1 (MQ=255) tAATAGAGACATCAATAATTTTATATTCACTGAAATTATTTTTAATCTTTATTTATAGCTGTT‑TGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCACTGCTATAAATTTAATCAATCACAAAGAT‑CGGTATATTCATAGGATATTTATCTGgt < 2:305717/149‑1 (MQ=255) aaTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTGTAGCTGTG‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAGACTTTGTTGCATATTCATAGGATATTTCTGTGgtt > 1:22429/1‑149 (MQ=255) aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTAGATCTTT‑T‑T‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt > 1:125384/1‑149 (MQ=255) aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt > 1:434651/1‑149 (MQ=255) aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GCTTCTT‑A‑T‑TCTTTT‑GGAGTTTGGTTGCGCGGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt > 1:190389/1‑149 (MQ=255) aaaCAGCAATAATTTTATATTGCCTGAAAATATTTTTAATCTTTATTTAGAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt < 1:87218/149‑1 (MQ=255) aaaaGCAATAATTTTATATTCACTGAAAATATTTTTAATCTCTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtt < 2:413953/146‑1 (MQ=255) aatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatt < 2:440644/149‑1 (MQ=255) ttttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattta < 1:543929/149‑1 (MQ=255) tttATATTAACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGATGCTATAAATTGAAT‑AATTAAAACTTTATTGCATTATCATAGGTTATTTGTCTTGTTTTTGTATGTGttatttat > 2:390900/1‑149 (MQ=255) aaaaTATTTTTAATCTTTATTTATAGCTGTG‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAAGTTTTTTTCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATTCGTATGCtt > 2:566488/1‑149 (MQ=255) aaaaTAGTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGGTGCGATAAATTTAAT‑AATTTAAAGGTTGTTGCATATTTGTAGGATGTTTATGTGGTTTGTGTTTGTGATATTTATTATGCGTATGCtt < 2:386063/149‑1 (MQ=255) tatTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCaaa < 1:242151/149‑1 (MQ=255) tattCTTTGTTTATAGGTGTT‑GTTTATT‑ATT‑T‑TTTT‑TGAGTTTGTTCGTGTTGCTATAAATTGAAA‑TATTAAAACTTCGTTGCATAGTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa < 2:392286/146‑1 (MQ=255) tAATTTGTATTTAGAGCTGTT‑GGTTATTGT‑T‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTTAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGCTATTTATTATGCGTATGCTTCAAAAACaaa < 2:392381/149‑1 (MQ=255) aaTCATTATTTATAGCTGTT‑TGTTATT‑A‑T‑T‑TTTTGGGAGTTGGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATTTTCATAGGATATTTAGCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaaa < 2:597122/149‑1 (MQ=255) tCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAATCAAACtt > 2:191827/1‑149 (MQ=255) tCTTTATTTATAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAAtt > 1:428884/1‑149 (MQ=255) tttatttaTAGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTAt < 1:517736/149‑1 (MQ=255) tatttatGGCGGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt < 1:177306/149‑1 (MQ=255) aGCTGTT‑GGTTATT‑ATT‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa > 2:644609/1‑135 (MQ=255) aGCTGTT‑GGTTATT‑A‑TCT‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaa < 1:644609/135‑1 (MQ=255) tGTT‑GGTTATT‑ATT‑T‑TTTT‑GCAGTTTGGTTGCGGTGCTATCCTTTGAAT‑ATTTAAAACTTTACTACATCACAAATCGTGATTCACTCAATTGGTTTTCTTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGtt < 1:318232/149‑1 (MQ=255) | GTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTT‑GGTTATT‑A‑T‑T‑TTTT‑GGAGTTTGGTTGCGCTGCTATAAATTGAAT‑AATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTT > NC_000913_3_hsa_tpiA/389778‑390055 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |