I started mothur by following the SOP with the example data downloaded from the wiki page. I run exactly the same commands (except te parameter processors):
mothur>make.file(inputdir=., type=fastq, prefix=stability)
mothur>make.contigs(file=stability.files)
mothur>summary.seqs(fasta=stability.trim.contigs.fasta)
At this point I found an inconsistency in the # Seqs: in my case it was 152307 compared to the 152360 shown in the example of the wiki. The numbers should be exactly the same, right? I continued anyway
mothur>screen.seqs(fasta=stability.trim.contigs.fasta, group=stability.contigs.groups,
maxambig=0, maxlengt h=275)
mothur>summary.seqs(fasta=stability.trim.contigs.good.fasta)
Here I get another inconsistency in numbers (in the Minimum):
Start End NBases Ambigs Polymer NumSeqs
Minimum: 1 29 29 0 3 1
2.5%-tile: 1 252 252 0 3 3221
25%-tile: 1 252 252 0 4 32208
Median: 1 252 252 0 4 64416
75%-tile: 1 253 253 0 5 96624
97.5%-tile: 1 253 253 0 6 125611
Maximum: 1 270 270 0 12 128831
Mean: 1 252 252 0 4
number of Seqs: 128831
unique.seqs(fasta=stability.trim.contigs.good.fasta)
count.seqs(name=stability.trim.contigs.good.names, group=stability.contigs.good.groups)
[ERROR]: is not in your groupfile, please correct.
I´m running mothur 1.42.0 in a virtual machine with ubuntu 14.04