This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to annotate novel SNPs using coding sequence position?

Hello!

I have a list of novel mutations with coding sequence positions (not in vcf format). I want to annotate them to know how many of them are deleterious.

As per my knowledge we require the genomic positions to use in for annotation in ANNOVAR (creating a vcf file or annovar format), however since I have coding sequence positions I am not sure how to proceed.

Any suggestions how to do it?

Thank you in advance.

snp
• 2,592 views • link
updated by Emily Emily • written by bioinfo89 • bioinfo89

3 answers

Are they in HGVS? You can use that as input for the VEP.

• 0 views • link
written by Emily Emily

I would recall your SNPs in a standard format - VCF. Then use Annovar, Ensembl VEP or snpEff to call. Otherwise you're going to have to script it with biopython or similar. That's going to be a world of pain considering it's essentially a solved problem.

• 0 views • link
written by colindaven colindaven

Hi colindaven,

I have a novel mutation list complied from literature review, so I cannot recall the SNPs. I wanted to know if I could use these mutation list with coding sequence positions (since genomic position was not available) to annotate.

Thanks!

• 0 views • link
written by bioinfo89 • bioinfo89

Sounds like you'll have to convert to a standard format. Try a minimal VCF version, it's probably easiest and best documented (in my experience). You probably don't need the INFO or GT columns. Good luck.

• 0 views • link
written by colindaven colindaven

Log in to answer this question.

More posts like this