# NOT RUN {
file_path <- system.file("extdata", "ncbi_basidiomycetes.fasta", package = "metacoder")
sequences <- ape::read.FASTA(file_path)
genbank_ex_data <- extract_taxonomy(sequences,
regex = "^.*\\|(.*)\\|.*\\|(.*)\\|(.*)$",
key = c(gi_no = "obs_info", "obs_id", desc = "obs_info"),
database = "ncbi")
# }
Run the code above in your browser using DataCamp Workspace