# NOT RUN {
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(seq_name = "obs_info", seq_id = "obs_info",
other_id = "obs_info", "class"),
class_regex = "^(.*)__(.*)$",
class_key = c(unite_rank = "taxon_info", "name"),
class_sep = ";")
# }
# NOT RUN {
# }
Run the code above in your browser using DataCamp Workspace