remove the singletons using remove.rare in mothur
Entering edit mode
6.2 years ago

Can you please advise if the following steps will do the job? Many thanks. If not, could you please also give us some suggestions (especially which commands should be used, and at which steps). Thanks.

I did these: Following mothur SOP,
1. mothur > sub.sample(, size=6105) Sampling 6105 from each group. 0.03

Output File Names:

  1. mothur > remove.rare(, nseqs=1) 0.03

Output File Names:

  1. mothur > count.groups( Total seqs: 207161.

Output File Names:

Get the minimal number of OTUs that were present in all of these samples, then go from there for : collect.single(, calc=chao-invsimpson, freq=100)

Does this insertion make senses? I mean, if I want to compare the results with (these procedures) and without the removing the singletons (by using the original mothur SOP).

Have a great day!

next-gen sequencing • 2.1k views

Login before adding your answer.

Traffic: 840 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6