View source: R/parse_amptk_taxonomy_table.R
parse_taxonomy_amptk_batch | R Documentation |
Parse multiple taxonomy strings from AMPtk.
parse_taxonomy_amptk_batch(x, withID = TRUE, multithread = FALSE)
x |
Vector of character strings from AMPtk taxonomy |
withID |
Logical; set to TRUE (default) if charter strings contain sequence IDs |
multithread |
Logical; run the function in parallel |
This function splits AMPtk (ex-UFITS) result of taxonomic annotation by taxonomic rank and combines it into one table. The first part of taxonomy string is assumed to be an ID of the best hit (e.g., EF040844) or a taxonomic classificator (UTAX or SINTAX) followed by a semicolon.
Data frame with OTUs or species as rows and their taxonomic ranks as columns.
AMPtk: Amplicon ToolKit for NGS data (formally UFITS): https://github.com/nextgenusfs/amptk
parse_taxonomy_amptk
parse_taxonomy_qiime
## Taxonomy with sequence ID
tax <- c(
"UDB011192;k:Fungi,p:Basidiomycota,c:Agaricomycetes,o:Russulales,f:Russulaceae,g:Russula,s:Russula chloroides",
"SINTAX;k:Fungi,p:Chytridiomycota,c:Chytridiomycetes,o:Rhizophydiales",
"EF040844;k:Fungi,p:Basidiomycota,c:Agaricomycetes,o:Atheliales,f:Atheliaceae,g:Piloderma",
"EU240039;k:Fungi,p:Zygomycota,o:Mortierellales,f:Mortierellaceae,g:Mortierella",
"UTAX;k:Fungi,p:Basidiomycota,c:Agaricomycetes,o:Atheliales,f:Atheliaceae,g:Amphinema,s:Amphinema byssoides",
"SINTAX;k:Fungi,p:Chytridiomycota,c:Chytridiomycetes,o:Spizellomycetales"
)
parse_taxonomy_amptk_batch(tax)
## Taxonomy without sequence ID
tax2 <- c(
"k:Fungi,p:Basidiomycota,c:Agaricomycetes,o:Russulales,f:Russulaceae,g:Russula,s:Russula chloroides",
"k:Fungi,p:Chytridiomycota,c:Chytridiomycetes,o:Rhizophydiales",
"k:Fungi,p:Basidiomycota,c:Agaricomycetes,o:Atheliales,f:Atheliaceae,g:Piloderma",
"k:Fungi,p:Zygomycota,o:Mortierellales,f:Mortierellaceae,g:Mortierella",
"k:Fungi,p:Basidiomycota,c:Agaricomycetes,o:Atheliales,f:Atheliaceae,g:Amphinema,s:Amphinema byssoides",
"k:Fungi,p:Chytridiomycota,c:Chytridiomycetes,o:Spizellomycetales"
)
parse_taxonomy_amptk_batch(tax2, withID = FALSE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.