Question: getGenome {biomartr} not working for viruses?
0
gravatar for peter.durr
13 months ago by
peter.durr0
peter.durr0 wrote:

Hi everyone

I am trying to retrieve various reference genomes of viruses held in Genbank, or more specifically NCBI's RefSeq database.

I am exploring the new package biomartr 0.7.0 - which does just this - i.e. retrieves reference genomes

whilst this seems to work well for most taxons - it is not performing for viruses

here is code in which I compared a bacterial pathogen (Mycobacterium tuberculosis) with a virus (rabies virus)

is.genome.available(organism = "Mycobacterium tuberculosis", db = "refseq")

is.genome.available(organism = "Rabies lyssavirus", db = "refseq")

both of these came back with the response "TRUE"
 

However, when I try fetching these two reference sequences using the following code:

Mtb.genome.refseq <- getGenome(db = "refseq", organism = "Mycobacterium tuberculosis",  path   = file.path("_ncbi_downloads","genomes"))

RABV.genome.refseq <- getGenome(db = "refseq", organism = "Rabies lyssavirus", path   = file.path("_ncbi_downloads","genomes"))

The M. tunberculosis fetch is successful, but the rabies virus fetch fails with the following message:

----------> No reference genome or representative genome was found for 'Rabies lyssavirus'. Thus, download for this species has been omitted. Have you tried to specify 'reference = FALSE' ?

adding the argument reference = FALSE does not solve the problem.

Thanks

PS. I am running the program in RStudio 1.1.447, R 3.5.0, Windows 10.

 

ADD COMMENTlink modified 13 months ago • written 13 months ago by peter.durr0
Answer: getGenome {biomartr} not working for viruses?
0
gravatar for Martin Morgan
13 months ago by
Martin Morgan ♦♦ 23k
United States
Martin Morgan ♦♦ 23k wrote:

biomartr is a CRAN package, so you should ask on it's support site.

ADD COMMENTlink written 13 months ago by Martin Morgan ♦♦ 23k

Thanks Martin

I had just assumed that biomatr was a Bioconductor package - apologies

I have moved my question to the developers github Issues site

Regards

ADD REPLYlink written 13 months ago by peter.durr0
Answer: getGenome {biomartr} not working for viruses?
0
gravatar for peter.durr
13 months ago by
peter.durr0
peter.durr0 wrote:

As Martin Morgan correctly pointed out, biomartr is a CRAN pacakge

I have moved my question to the developer's github Issues site

 

 

ADD COMMENTlink written 13 months ago by peter.durr0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 260 users visited in the last hour