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