I am having an issue with the trim.seqs command in an ITS workflow. I do not have an oligos file to work with, because the people at my sequencing facility have already trimmed the primers and barcodes from my sequences. So what I’m trying to do is use the trim.seqs function with the maxlength as my parameter. At first, I set the maxlength to 250. The command ran, but when I tried to run my summary table I got an error due to missing values. I then ran trim.seqs again, but this time with a maxlength of 300 basepairs. This time the summary.seqs command for my output worked, but it DRASTICALLY reduced my number of sequences to a point that I think it wrong.
Here is my summary.seqs for the output of my screen.seqs function:
And here is the summary.seqs for the output of my trim.seqs after I have set the maxlength to 300 basepairs:
I would dearly appreciate any advice about how to proceed with this trim.seqs command.
All the best