## 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")
# ## End(Not run)
Run the code above in your browser using DataLab