This checks if the _original_ phyloseq object's columns represent taxa. If true (taxa_are_rows == FALSE), then the new OTU table is transposed before returning. This code makes the asummption that newsub was tranposed to have rows as taxa.

Or maybe I am missing something.