This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Map rs to gene symbol in python script

I have snp ids (rs18008665, rs3212227...) Now, I want to map these snp ids to gene symbols using Python script. I have found similar post, but as I have recently started to learn programing and with no bioinformatics knowledge, I´m unable to create the python script which works. Any script example to use?

snp python script gene
• 2,282 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by leireahedo02 • leireahedo02

1 answer

using bash

$ wget  -O  - "http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=snp&id=rs18008665&retmode=xml"  |\
xmllint --format --xpath  '//*[name(.)="FxnSet"]/@symbol' - |\
tr " " "\n" | sort | uniq | cut -d '"' -f 2
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thank you Pierre. Easy to understand and concise solution, great help

• 0 views • link
written by leireahedo02 • leireahedo02

Log in to answer this question.

More posts like this