I've got a list of rs SNPs I'd like to enter into the David functional annotation tool. Since it does not support rsids I need to get the refseq gene id (or something similar) of the genes these SNPs overlap with first.
Is there a simple way of getting a gene ID for a SNP? Solution in BioPython or R is fine.
Since I have 22k snps I need an automatic way. And if a text file that maps these values exist, a link to it would be enough.
Ps. preferably looking for a solution that does not use the position of the SNPs; I would be able to solve the problem this way myself.