Mothur

Mothur pipeline to group sequences 99% similar

unique.seqs(fasta=felix.fasta)

filter.seqs(fasta=felix.unique.fasta, trump=.)

dist.seqs(fasta=felix.unique.filter.fasta)

read.dist(column=felix.unique.filter.dist, name=felix.names)

cluster(cutoff=0.01)

get.oturep(column=felix.unique.filter.dist, name=felix.names, fasta=felix.fasta, list=felix.unique.filter.fn.list)