Question

customProDB: issues getting dbSNP data for human (hg19)

0

Entering edit mode

kristenbeck527 • 0

@kristenbeck527-8349

Last seen 9.4 years ago

United States

I am trying to use the variant annotation capabilities in customProDB. I would like to make a variant annotated fasta file, but I am running into issues retrieving the dbsnp data for hg19.

PrepareAnnotationRefseq(genome='hg19', CDSfasta, pepfasta, 
     annotation_path=annotation_path_hs, splice_matrix = FALSE, 
     dbsnp = "snp130", COSMIC = TRUE)

I also tried "snp131" as recommended in the vignette and "snp135" as recommended in the help page of the function. All snp variations returned the error below.

Error in PrepareAnnotationRefseq(genome = "hg19", CDSfasta, pepfasta, : invalid dbsnp name for specified genome

Any help would be appreciated. Thank you!

customprodb dbsnp refseq varianttools • 1.8k views

ADD COMMENT • link updated 9.4 years ago by xiaojing.wang ▴ 50 • written 9.4 years ago by kristenbeck527 • 0

score 2 · Accepted Answer · 2015-07-08

Try using ‘snp138’, ‘snp141’ or ‘snp142’ instead. You may use the code below to check the current dbSNP version provided by the UCSC table browser. I’ll update the document accordingly. library(rtracklayer) session <- browserSession() genome(session) <- 'hg19' tmp <- trackNames(session) tmp[grep('snp', tmp)] xiaojing From: kristenbeck527 [bioc] [mailto:noreply@bioconductor.org] Sent: Wednesday, July 08, 2015 5:10 PM To: Wang, Xiaojing Subject: [bioc] customProDB: issues getting dbSNP data for human (hg19) Activity on a post you are following on support.bioconductor.org<https: support.bioconductor.org=""> User kristenbeck527<https: support.bioconductor.org="" u="" 8349=""/> wrote Question: customProDB: issues getting dbSNP data for human (hg19)<https: support.bioconductor.org="" p="" 69670=""/>: I am trying to use the variant annotation capabilities in customProDB. I would like to make a variant annotated fasta file, but I am running into issues retrieving the dbsnp data for hg19. PrepareAnnotationRefseq(genome='hg19', CDSfasta, pepfasta, annotation_path=annotation_path_hs, splice_matrix = FALSE, dbsnp = "snp130", COSMIC = TRUE) I also tried "snp131" as recommended in the vignette and "snp135" as recommended in the help page of the function. All snp variations returned the error below. Error in PrepareAnnotationRefseq(genome = "hg19", CDSfasta, pepfasta, : invalid dbsnp name for specified genome Any help would be appreciated. Thank you! ________________________________ Post tags: customprodb, dbsnp, refseq, varianttools You may reply via email or visit customProDB: issues getting dbSNP data for human (hg19)