I’m trying to work out why some of our sequences aren’t classifying, it’d be really useful if Mothur could (optionally) write out .fasta files of classified or unclassified reads.
You can run get.lineage and set “unclassified” as your taxon.
Thank you Pat, that worked and has saved me writing a script. 