Add species-level annotation to a taxonomic table.

Share:

Description

addSpecies wraps the assignSpecies function to assign genus-species binomials to the input sequences by exact matching against a reference fasta. Those binomials are then merged with the input taxonomic table with species annotations appended as an additional column to the input table. Only species identifications where the genera in the input table and the binomial classification are consistent are included in the return table.

Usage

1
addSpecies(taxtab, refFasta, allowMultiple = FALSE, verbose = FALSE)

Arguments

taxtab

(Required). A taxonomic table, the output of assignTaxonomy.

refFasta

(Required). The path to the reference fasta file, or an R connection. Can be compressed. This reference fasta file should be formatted so that the id lines correspond to the genus-species binomial of the associated sequence:

>SeqID genus species ACGAATGTGAAGTAA......

allowMultiple

(Optional). Default FALSE. Defines the behavior when multiple exact matches against different species are returned. By default only unambiguous identifications are return. If TRUE, a concatenated string of all exactly matched species is returned. If an integer is provided, multiple identifications up to that many are returned as a concatenated string.

verbose

(Optional). Default FALSE. If TRUE, print status to standard output.

Value

A character matrix one column larger than input. Rows correspond to sequences, and columns to the taxonomic levels. NA indicates that the sequence was not classified at that level.

See Also

assignTaxonomy, assignSpecies

Examples

1
2
3
4
5
## Not run: 
 taxtab <- assignTaxonomy(dadaF, "rdp_train_set_14.fa.gz")
 taxtab <- addSpecies(taxtab, "rdp_species_assignment_14.fa.gz")

## End(Not run)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.