I have a basic question. I got 16S MiSeq paired microbiome dataset and I did taxonomy annotation using mothur MiSeq SOP. In my annotation process, I used silva.v4.fasta as reference in align.seqs() step and for classify.seqs() step I used trainset18_062020.rdp.fasta as reference. Did I make mistake doing like this?
For my research, I use the silva for the align step, but for classify I use silva, PR2 or whatever is better suited for the group I am studying (protists/prokaryotes/metazoans/etc).