Annotate Ref/Alt for a list of rsIDs?
Hi,
I have a large list of rsIDs from dbSNP. I would like to get the ref/alt for each of the snps. What would be the best way of handling this? It's a VCF containing position and dbSNP id, but no ref/alt columns.
Thanks.
• 4,275 views
•
link
updated
by
GenoMax
•
written
by
blazer9131 •
•
written
by
blazer9131 •
1 answer
bcftools annotate -a dbsnp.vcf -c ID -o output.vcf input.vcf.compressed.indixed.vcf.gz should work
• 0 views
•
link
updated
by
GenoMax
•
written
by
binodregmi30 •
•
written
by
binodregmi30 •
Log in to answer this question.
More posts like this
-
Converting custom VCF to standard VCF
written by avelarbio46 •
-
Filtering VCF variants based on txt file with variants CHROM, POS and ID
written by avelarbio46 •
-
Best way to list variants in UK Biobank bed files
written by davidjohngreen •
-
Liftover dbSNP IDs
written by Arko •
-
Converting Affymtrix SNP data to plink and VCF formats
written by robjohn70000
-
Help understanding CHROM and POS fields in VCF
written by bdolin
-
best way of filtering a VCF file using a list of SNP IDs and ref/alt alleles
written by auraf85 •
-
Best way to query VCF for specific variants
written by Sean
-
Get gene names from rs SNP ids
written by Endre Bakken Stovner
-
Computing B-Allele Frequency From Bam Or Vcf Related To Exome Samples
written by Stephwen
How large is large? Give us a number.
Also, I assume single species?
It's roughly 600k IDs. And yes, human.
My first thought would be to download the latest version of dbSNP in VCF format and then use Python, grep, or BCFtools to pull out the information that you need.
dbSNP in VCF can be downloaded from here: https://www.ncbi.nlm.nih.gov/variation/docs/human_variation_vcf/
I reckon that it would take a good few hours to pull out all info.
Kevin