How to download all available sequences of a gene from all bacteria using R
I need to download all/many sequences of a specific bacterial gene from Genbank nuccore database from entries limited to complete genome sequences. I prefer using R. Querying: 'Bacteria[ORNG] AND gyrB[GENE] AND complete genome[TI] ' in web interface results in >10k hits. I do not want to download whole genome sequences but only extracted gyrB sequences to make a local database. I tried
library(rentrez):
db = "nuccore"
query = "Bacteria[ORGN] AND gyrB[GENE] AND complete[TI]"
found = entrez_search(db, query, config = NULL, retmode = "xml", use_history = FALSE, retmax = 90000)
but this fetch ids for whole genome sequences. Is it possible to get fasta sequences for gryB genes or at least gyrB coordinates however I'm not into downloading whole genome sequences of thousands of genomes.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
efetch from NCBI E-utilities returns "curl error s 400 & 500" and takes a very long time
written by Eugene •
-
How to download large protein data from NCBI?
written by enriqp02 •
-
Perl script to download gene sequences
written by dllopezr •
-
How to avoid http failure when using rentrez to fetch records from a long list of species names
written by lvogel •
-
E-utilities for obtain gene sequences from the gene database
written by dllopezr •
-
rentrez fetch error - unable to obtain query
written by klrdna •
-
Entrez.efetch HTTP Error 400: Bad request
written by wanderingstefan •
-
Download complete bacterial genomes and associated plasmid sequences from NCBI
written by wanderingstefan •
-
Fetching multiple FASTA from NCBI
written by mikael.lenz.strube •
-
NCBI - Gene File to GenBank file in R
written by jaschelite •
You can get this data from Ensembl bacteria using the Ensembl Genomes perl API or maybe using the R package biomartr.
It would be a great option but I found that BioMart is not currently available for Ensembl Bacteria. https://support.bioconductor.org/p/82585/