Obtaining The Snp Rs Number With The Chromosomal Position
3
6
Entering edit mode
9.0 years ago
Sheila ▴ 390

This question is similar to this one (Get rs number based on position).

I have a text file with SNPs in the chr:position format

 10:71086
10:72876
10:75794


I was wondering if there is an R package (or perhaps one in Bioconductor) that can take these as input values and provide the SNP rs number associate with them????

r snps bioconductor • 41k views
2
Entering edit mode
9.0 years ago

I think JC has already pointed out a nice resource. The other option would be

Get rs number based on position

4
Entering edit mode
9.0 years ago
JC 13k

It's not in R, but Kaviar (http://db.systemsbiology.net/kaviar/) can give you the RS information using the coordinates and can be accessed locally or remotely.

1
Entering edit mode
3.5 years ago
hsiaoyi0504 ▴ 60

I would like to share another approach. Using EDirect (https://dataguide.nlm.nih.gov/edirect/install.html).

Then you can get the result by running following command:

esearch -db snp -query chr_id[CHR] AND organism[ORGN] AND chr_pos[CPOS]

(modify chr_id, organism, chr_pos to fit your application).

For python, take a look of this function I wrote before https://github.com/hsiaoyi0504/bcinfo/blob/546c54c2f72febd361da59e18b9901591837f826/bcinfo/search.py#L5.

0
Entering edit mode

0
Entering edit mode

Here is another Python wrapper for easy handling of the searches, but without the need to install edirect: https://github.com/krassowski/easy-entrez#example-obtaining-the-snp-rs-id-number-from-chromosomal-position