# NOT RUN {
file_path <- system.file("extdata", "unite_general_release.fasta", package = "metacoder")
sequences <- ape::read.FASTA(file_path)
unite_ex_data_1 <- extract_taxonomy(sequences[!grepl(pattern = "\\|UDB", names(sequences))],
regex = "^(.*)\\|(.*)\\|(.*)\\|.*\\|(.*)$",
key = c(name = "obs_info", "obs_id",
other_id = "obs_info", tax_string = "obs_info"),
database = "ncbi")
# }
# NOT RUN {
# }
Run the code above in your browser using DataCamp Workspace