"read.taxonomy" = function(taxon.table, eda.dis, species = 1, family = 6, taxon.seq = 9)
{
	return(taxon.table[dimnames(eda.dis)[[1]], c(species, family, taxon.seq)])
}