This function is used internally as part of the infrastructure to ensure that
component data types in a phyloseq-object have exactly the same taxa/species.
It relies heavily on the
Reduce function to determine the
strictly common species.
Returns a character vector of only those species that are present in
all species-describing components of
1 2 3
# ## data(GlobalPatterns) ## head(intersect_taxa(GlobalPatterns), 10)