library(BSgenome.Hsapiens.UCSC.hg19)
library(SNPlocs.Hsapiens.dbSNP.20120608)
snps.file <- system.file("extdata", "pca.enhancer.snps", package = "motifbreakR")
snps <- as.character(read.table(snps.file)[,1])
snps.mb <- snps.from.rsid(snps,
dbSNP = SNPlocs.Hsapiens.dbSNP.20120608,
search.genome = BSgenome.Hsapiens.UCSC.hg19)
Run the code above in your browser using DataLab