How can I generate a fasta file with all sequences assigned to a given otu?
The file I have is an *.an.list, generated at the stage:
mothur>cluster.split(fasta=final.fasta, taxonomy=final.taxonomy, name=final.names, taxlevel=6, processors=2)
We have a large set of 250nt long sequences of a variable region of the 16S gene in the first otu0001 at the 0.01 difference level,
and we would like to run unique.seqs on the obtained fasta file to check the precise taxonomic affiliation best hit of the different unique sequences in the otu.
Thanks in advance,