I lunched a trim.seqs command, and Mothur returned an error telling me that my sequences are in my fasta file and not in my quality file. The previous command (sort.seqs) did work properly and found the same amount of sequences in both files.
To check, I tried to lunch a summary.qual command but mothur the error “M02442-223 has a quality scores of 80, expecting values to be less than 40.” for each sequences until the command crashes.
I never had any problem with the quality files on this dataset before…
Any ideas on what went wrong and how to fix it?
Previous commands were fastq.info, make.contigs, screen.seqs and sort.seqs. All apparently working withour errors.