make.contigs(file=nemabiome.files, processors=2)
screen.seqs(fasta=nemabiome.trim.contigs.fasta, minlength=200, maxlength=450, maxambig=0, group=nemabiome.contigs.groups, processors=2)
align.seqs(candidate=nemabiome.trim.contigs.good.fasta, template=nematode_ITS2_database_version1_3.fasta)