This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SNP ID

Hello Community,

I'm Shamsur. I was dealing with variant calling.

Is it possible to get snp id's while generating vcf file, like this image vcf file for bangladeshi goat breed

I got this from my professor. It represent bangladeshi goat breed snp. Here this vcf file has ID for snp's and also dbSNP_143 was used to get these ID's.

Now here is what I got while doing varient calling. In my case ID has a dot(.). I did not get any ID for this vcf file. Someone please help me, how can I get vcf file with ID's

enter image description here

Thank you.

snp id dbsnp_143
• 2,111 views • link
written by shamsurrahaman132 • shamsurrahaman132

If I'm not wrong, in GATK, you can pass the -D argument in HaplotypeCaller and pass your DBSNP file as input. But, you can always annotate your vcf file after generating it!

• 0 views • link
written by brunobsouzaa brunobsouzaa

Thank you. @brunobsouzaa.

I did look at GATK documentation and your suggation worked.

• 0 views • link
written by shamsurrahaman132 • shamsurrahaman132

1 answer

Most variant callers don't look up variant ID. You would need to run your VCF through a variant annotator, eg VEP to add ID and other useful annotation.

• 0 views • link
written by Emily Emily

Thank you @Emily_Ensembl

This one is new for me. I will check it out.

• 0 views • link
written by shamsurrahaman132 • shamsurrahaman132

Log in to answer this question.

More posts like this