I couldn’t find stability.contigs.count_table to run commands like summary.seqs(fasta=stability.trim.contigs.fasta, count=stability.contigs.count_table). Please see attached the files make.contigs created. I didn’t have this problem earlier when I used Mothur. I tried various versions but none output count_table from make.contigs. Is it Ok to proceed?
I suspect you have an older version of mothur. If you open that mothur logfile that’s at the top you’ll find the version number towards the top. This was a change that was made in recent versions of mothur.
Hopefully one last question, I am running 12 samples to set up Mothur to run some clinical experiment samples. I followed the Miseq SOP. Could you tell me why the dist.seqs(fasta=final.fasta, cutoff=0.03) is taking ever to finish?Is it expected or am I doing something stupid? I am running it on a supercomputer using 900 GB of memory.
I suspect you have a ton of sequences. Are you sequencing the V4 region or another region? Regardless, you should probably checkout using cluster.split instead. You can see how to use this in the mothur MiSeq SOP.