Hi, I’m working with trainset9_032012 but it doesn’t work as i expected because the result is aproximatelly 0.01% of my starting sequences.
At the beggining i had 41614 sequences and latter, i used the command “classify.seqs(fasta=stability.trim.contigs.good.unique.precluster.pick.fasta, count=stability.trim.contigs.good.unique.precluster.denovo.vsearch.pick.count_table, reference=trainset9_032012.pds.fasta, taxonomy=trainset9_032012.pds.tax, cutoff=80)” and the result were 501 sequences. Does anybody knows about a trainset that can i use?
Thanks.
Hi - yeah the trainset9 you are using is the RDP’s 16S rRNA gene sequence database, not the ITS. There is a UNITE database out there that you can get for fungal ITS sequences.
Pat