« 回覆文章 #1 於: 五月 05, 2017, 01:55:37 am »
http://stackoverflow.com/questions/41490657/best-way-to-get-list-of-snps-by-gene-id/41496733library(biomaRt)
#load the human variation data
variation = useEnsembl(biomart="snp", dataset="hsapiens_snp")
#look up a single gene and get SNP data
getBM(attributes = c(
"ensembl_gene_stable_id",
'refsnp_id',
'chr_name',
'chrom_start',
'chrom_end',
'minor_allele',
'minor_allele_freq'),
filters = 'ensembl_gene',
values ="ENSG00000166813",
mart = variation
)