## Not run:
#
# file_path <- system.file("extdata", "rdp_current_Archaea_unaligned.fa", package = "metacoder")
# sequences <- ape::read.FASTA(file_path)
# rdp_ex_data <- extract_taxonomy(sequences,
# regex = "^(.*?) (.*)\\tLineage=(.*)",
# key = c(id = "obs_info", description = "obs_info", "class"),
# class_regex = "(.+?);(.*?);",
# class_key = c("name", "taxon_info"))
# ## End(Not run)
Run the code above in your browser using DataLab