This is my first time with Mothur! I’m analyzing some 16S rRNA samples from an anoxic basin and I need some help.
My samples contain a mix of archaeal and bacterial sequences, and I want to classify and split them using the SILVA database as a reference. I tried classify.seq and remove.seq with the archaeal database and all went fine. I tried to do the same with the bacterial one, but the program gave me an error, telling me that it couldn’t match some sequences in the reference database (silva.gold.ng.fasta) with those in the taxonomic reference (silva.bacteria.silva.tax). I checked both files and found that the identifiers in the taxonomy files are different from those in the fasta file:
Taxonomy: AB000389.1 Bacteria;Proteobacteria;Gammaproteobacteria;Alteromonadales;Pseudoalteromonadaceae;Pseudoalteromonas;
FASTA: 7000004128189528 [sequence]
How can I fix this error?