I am currently working with the Phyloseq package in R to analyze 16s sequences. I sequenced the samples with Illumina's MiSeq platform.
I want to see if I have pathogens in my samples and want to BLAST all the sequences from the Bacteroidetes OTUs but am having trouble figuring out how to determine which OTUs were bacteroidetes. Does anyone know the R commands for phyloseq that would allow me to do this so that I can go back to my QIIME outputs and get the sequences from those OTUs?
I want to see if I have pathogens in my samples and want to BLAST all the sequences from the Bacteroidetes OTUs but am having trouble figuring out how to determine which OTUs were bacteroidetes. Does anyone know the R commands for phyloseq that would allow me to do this so that I can go back to my QIIME outputs and get the sequences from those OTUs?