Retrieve gene sequences from NCBI by accession number.
ncbi_byid(ids, format = NULL, verbose = TRUE)
(character) GenBank ids to search for. One or more. Required.
(character) Return type, e.g., "fasta". NOW IGNORED.
"fasta"
(logical) If TRUE (default), informative messages printed.
TRUE
Data.frame of the form:
taxon - taxonomic name (may include some junk, but hard to parse off)
gene_desc - gene description
gi_no - GI number
acc_no - accession number
length - sequence length
sequence - sequence character string
If bad ids are included with good ones, the bad ones are silently dropped. If all ids are bad you'll get a stop with error message.
ncbi_search, ncbi_getbyname
ncbi_search
ncbi_getbyname
# NOT RUN { # A single gene ncbi_byid(ids="360040093") # Many genes (with different accession numbers) ncbi_byid(ids=c("360040093","347448433")) # }
Run the code above in your browser using DataLab