I was trying to run the trim command with fasta and qual files I have received from the sequencing service.
Mothur recognize the fasta file but cannot recognize the qual file. I get the following message for each sequence:
"We found more than 25% of the bases in sequence IAQWM6U02HEYNX to be ambiguous. Mothur is not setup to process protein sequences
I tried again to run the trim command - this time it recognized the file and did the trim.
However - when I tried to check the qual file with summary.seq command Mothur cannot read the file. I use to run the summary.seq command before I start the analysis just to see that there is a match between the fasta and the qual.
So, the command was: summary.seqs(fasta=XXX.qual) - and it gave me an error. In older versions of Mothur there was no problem to run this command
I HAVE SAME PROBLEM i used 1.48 mothur the command i used
make.contigs(file=stability. Files)
it should be give me two file fastaand qual but that is not happen !!
i try summary.qual(qfile=summary.qual(qfile=stool.trim.pick.qual) but i change qfile=stability this is my file name and it does not work what is your recommendation
Hi - can you post htis to a new thread? This one is 10 years old and a lot has changed over that time. FWIW - make.contigs will give you a conesnsus fasta file and count file. It should not give a quality score file.