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 | 92.1% | 114.8 / 1.8 | 38 | 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 . (3/0); new base T (19/16); total (22/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.49e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
AAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑A‑TTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAAT > NC_000913_3_hsa_tpiA/389771‑390060 | aaGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttatt‑a > 1:598475/1‑149 (MQ=255) gCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAACCTTTATTTTTGGCTGTTGGTTatt‑c‑tt > 2:599898/1‑149 (MQ=255) tttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTCATCTTTATTTATAGCTTTTGGTTATT‑A‑TTTT‑TTTGGGttttcgt > 1:47309/1‑145 (MQ=255) gAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGCCACACCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTgcgc > 1:190390/1‑149 (MQ=255) caaaGTTGTCGAGATGTTGATAAGAATGGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGCAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTGGGTTgctc < 2:83086/148‑3 (MQ=255) atttttCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCgctgc < 1:3877/144‑1 (MQ=255) ttGTCGCGATGGTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATTTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCta < 2:24319/149‑1 (MQ=255) cgcgATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATaaa < 1:639190/149‑1 (MQ=255) gcgATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTCTAAAt > 2:268323/1‑149 (MQ=255) ttGACAAGAAGAAAATTGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGaata < 1:258321/149‑1 (MQ=255) gACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATCTGTTTAATCTTTCTTTCAAGCCTCGGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGaataat > 1:646992/1‑149 (MQ=255) gACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTTTTTTTTTGTTGCGCTGCTATCAGTTTcataat > 2:4698/1‑149 (MQ=255) gACAAGAAGAGAATCGAAGAGAGGCCAGGACATCTGTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTGATT‑ATTTTT‑TTGGAGTTTGGGTGCGCTGCTATAAATTGaataat < 2:679703/149‑1 (MQ=255) aCAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAAtt < 1:403440/149‑1 (MQ=255) aagaagAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTaa > 1:273281/1‑149 (MQ=255) gaagaGAATGGAAGTGAGGCCAGGAGATCTTTTTAATAGAAACATCAATAATTTTATATTCACTGTAAATATTTTTAATCTTTATTTATAGCTGTTGTTTATT‑ATTTTT‑TTGGAGTTCGGTTGCGCTGCTATAAATTGAATAATTaaaa < 2:582990/149‑1 (MQ=255) ggAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTAACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGc > 1:407375/1‑149 (MQ=255) agagagGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGATTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCatt > 2:391992/1‑148 (MQ=255) agagGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTTGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAATTTGAATAATTAAAACTTTGTTGCATCtt > 1:252989/1‑149 (MQ=255) agagGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGAATAtt > 1:188425/1‑149 (MQ=255) gCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATa > 1:687638/1‑149 (MQ=255) gCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTGGGTTATTCT‑TTTT‑TTGGCGTCTTCTCGTGCTGATATAAATTGAATGAGTAAAACTTTGTTGCATATTCTTa > 2:413691/1‑149 (MQ=255) gACAGCTTTTTAATAGAAACAGCAATAATTTGATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑A‑TTTTGTTGTAGGTTCGTTAGAGTGAGAAGAAGTGTGTAATTCAACCTTTGTTGCATATTCATATGata > 2:486130/1‑149 (MQ=255) tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATCGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGTTGCTATAAATTGAATAAGTAAAACTTTGTTGCATATTCTTAGTGTATTtctc > 1:613611/1‑149 (MQ=255) tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATc > 1:624736/1‑149 (MQ=255) tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATc < 1:537788/149‑1 (MQ=255) aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt > 1:169187/1‑149 (MQ=255) aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCGATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgttt > 1:592304/1‑149 (MQ=255) tAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGTTTTTTTTTTGTTATGCTATAAATTAGATAATCAAAACTTTGTTTCATATTAATATGATATTTATCTgttgtt > 2:29800/1‑144 (MQ=255) aaaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGTATTTTGGTTGCGCTGCTATAAATTGAATAATTAAAATTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt > 2:503632/1‑149 (MQ=255) aaaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt < 1:316219/149‑1 (MQ=255) agaaataatTTTATACTCAATGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtttgt < 2:434695/146‑1 (MQ=255) aataatTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCTCTGCTATAAATTGAATAATTAAAACTTTTTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGat > 1:235825/1‑149 (MQ=255) tttATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatttat < 1:274497/149‑1 (MQ=255) atatTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTtattat < 1:130950/149‑1 (MQ=255) ttCACTGAAAAGATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCg < 1:134665/149‑1 (MQ=255) aCTGAAAATATTTTTAATCTTTTTTTATAGCTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATGTTTATCTGGTTTGTGTTTGTGATTTTTATTATGCGTAt < 1:595582/149‑1 (MQ=255) atggttatttatAGCTGTTGGTTAGT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAAt < 2:93773/145‑1 (MQ=255) cTGTTGGTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGt > 1:51086/1‑149 (MQ=255) gTTATT‑ATTTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCACAATGtgggact > 1:428630/1‑143 (MQ=255) | AAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATT‑A‑TTTT‑TTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAAT > NC_000913_3_hsa_tpiA/389771‑390060 |
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: ‑ |