Question: Finding dbSNP rsIDs from a VCF file
gravatar for marcus.grewal
2.7 years ago by
marcus.grewal0 wrote:

I have a vcf file containing a couple thousand SNPs. None of them have a number in the ID column. I know that most of these SNPs do not have a dbSNP ID. However, I want to find the ID for the ones that do.

Is there any way to loop over the .vcf file and query a database to fill in the ID for the SNPs that are discovered?

snp dbsnp vcf • 1.8k views
ADD COMMENTlink modified 2.7 years ago by genecats.ucsc570 • written 2.7 years ago by marcus.grewal0

if you have dbSNP vcf and sample vcf, you can intersect both of them

ADD REPLYlink modified 2.7 years ago • written 2.7 years ago by cpad011212k
gravatar for genecats.ucsc
2.7 years ago by
genecats.ucsc570 wrote:

If your organism of interest is Human, then the Variant Annotation Integrator (VAI) at the UCSC Genome Browser can get these rsID's for you. Load your VCF file as a custom track, head to the VAI, make whatever selections you want, then make sure the "Include dbSNP rs# ID if one exists" box is checked, and click get output. You'll get some output with a bunch of info that maybe you aren't interested in, but any rsId's will be included. You can output these results to a file then filter for lines that have the rsId.

If you have any questions about using the VAI, or about the UCSC Genome Browser in general, feel free to send them to our mailing lists: for most questions, or if your question involves private data

ChrisL from the UCSC Genome Browser

ADD COMMENTlink written 2.7 years ago by genecats.ucsc570

Thanks for the informative answer!

ADD REPLYlink written 2.7 years ago by marcus.grewal0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1435 users visited in the last hour