Get rs ids using chromsome coordinates

0

Entering edit mode

3.7 years ago

Mahan ▴ 70

I have: CHR START END Effect_Allele Ref_Allele columns (hg38)

I want: CHR START END Effect_Allele Ref_Allele RS_ID

I would like to get the SNP rs id using chromosome coordinates, please let me know how I can do this?

Thanks in Advance

rs id Chromosome Coordinates Hg38 • 1.1k views

ADD COMMENT • link updated 3.6 years ago by Biostar 20 • written 3.7 years ago by Mahan ▴ 70

1

Entering edit mode

Download UCSC vcf (hg 38)
Intersect with bcftools or bed tools for ID column. Both bcftools and bedtools, support tab and vcf formats for intersection.

ADD REPLY • link 3.7 years ago by cpad0112 21k

0

Entering edit mode

You can download dbSNP in VCF format and then use bcftools to extract %CHROM\t%POS\t%REF\t%ALT\t%ID\n as a tab separated file from this VCF. Then, use R (data.table) to join your table to this table withIDs. That's how I'd do it.

ADD REPLY • link 3.7 years ago by Ram 43k

0

Entering edit mode

Hi, RamRS ~ Thanks for your suggestion. Do you mind sharing the detailed steps and R code?

ADD REPLY • link 3.7 years ago by Mahan ▴ 70

2

Entering edit mode

Yes, I do mind sharing that. Please invest some effort and let us know if you have any specific questions.