This is a test version of Biostars. For the public version, visit https://www.biostars.org.
customProDB: issues getting non-model organism dbsnp

I'm trying to use PrepareAnnotationEnsembl() function in customProDB. Based on Ensembl's site the current dbSNP version for Macaca mulatta is snp131. However, if I try to enter that into the function I get an error.

ensembl = useMart(biomart="ensembl")
ensembl = useMart("ENSEMBL_MART_ENSEMBL", dataset = "mmulatta_gene_ensembl", host="may2012.archive.ensembl.org")
PrepareAnnotationEnsembl(mart=ensembl, annotation_path=annotation_path_mm, splice_matrix = FALSE, dbsnp="snp131", COSMIC=TRUE)

I have also tried just using the number e.g. dbsnp="131". Both commands return the following error:

Error in pmatch(dbsnp, dbsnps) : object 'dbsnps' not found

Without adding the dnsnp parameter, it retrieves the annotation data. Also, in the source of the function, it is set up to query human and mouse dbSNP data. Could it be that non-model organisms aren't supported for dbSNP reference data retrieval? Any help on this would be greatly appreciated.

-Kristen

ensembl var snp customprodb

0 answers

No answers yet.

Log in to answer this question.