## Not run:
# # Extract embedded classifications from UNITE FASTA file offline
# file_path <- system.file("extdata", "unite_general_release.fasta", package = "metacoder")
# sequences <- ape::read.FASTA(file_path)
# unite_ex_data_3 <- extract_taxonomy(sequences,
# regex = "^(.*)\\|(.*)\\|(.*)\\|.*\\|(.*)$",
# key = c(name = "obs_info", seq_id = "obs_info",
# other_id = "obs_info", "class_name"),
# database = "none")
# # Look up taxonomic data online using sequence ID
# unite_ex_data <- extract_taxonomy(sequences,
# regex = "^(.*)\\|(.*)\\|(.*)\\|.*\\|(.*)$",
# key = c(name = "name", seq_id = "obs_id",
# other_id = "obs_info", tax_string = "obs_info"))
# ## End(Not run)
Run the code above in your browser using DataCamp Workspace