How to convert ss to rsid (SNPs)
1
0
Entering edit mode
3.3 years ago

Hi everybody! I have a long list of submitted SNPs (ss) and I would like to convert them to rsid nomenclature. I know I can do this one by one using dbSNP but how can I do this in batch mode? Thanks!!
snps • 1.1k views
ADD COMMENT •
0
Entering edit mode

Have you tried eutills/edirect?

ADD REPLY •
0
Entering edit mode

I'm gonna take a look into it, thanks!

ADD REPLY •
0
Entering edit mode

Can you provide a couple of examples?

ADD REPLY •
0
Entering edit mode

I have a list like this comprising 43.000 ss:

ss68767636
ss68769665
ss68769692
ss68769749
ss68771111
ss68772343
ss68774024
ss68776006
ss68779227

If I search ss68767636 on dbSNP I get -> rs1381888

ADD REPLY •
0
Entering edit mode

If you use the method method below be sure to sign up for and use NCBI_API_KEY. Add in a sleep 2 directive since you are processing thousands of queries.

ADD REPLY •
0
Entering edit mode

It worked! Thanks a lot!!

ADD REPLY •
0
Entering edit mode

Please accept GenoMax's answer below.

Upvote|Bookmark|Accept

ADD REPLY •
2
Entering edit mode
3.3 years ago
GenoMax 141k

Using EntrezDirect:

$ more id
ss68767636
ss68769665
ss68769692
ss68769749
ss68771111
ss68772343
ss68774024
ss68776006
ss68779227

$ for i in `cat id`; do printf ${i}"\n"; esearch -db snp -query ${i} | esummary | xtract -pattern DocumentSummary -element SNP_ID; printf "\n";done
ss68767636
1381888
1381888

ss68769665
1796935
1796935
1796935

ss68769692

ss68769749
227104
227104
227104

ss68771111

ss68772343
6683621
6683621
6683621
6683621

ss68774024
4658257
4658257
4658257
4658257

ss68776006

ss68779227
ADD COMMENT •

Login before adding your answer.

Similar Posts
Loading Similar Posts
Traffic: 2287 users visited in the last hour
Content Search
Users
Tags
Badges
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6