What's the diff between seqfilter & seqtrim?

Both of the function seems to remove sequence with bad quality, doesn't make them overlaped?

 Accepted Answer

Rik is correct, seqfilter will filter out entire reads from a FASTQ file, whereas seqtrim modifies the sequences themselves. seqtrim is needed because reads typically have lower quality at the termini, which can lead to poor or incorrect alignment. However, even after trimming reads, many reads may not contain enough high-quality bases to be included in an analysis, and must be filtered out.
A typical QC pipeline will need to do both to produce high-quality reads before alignment. seqqcplot can be a helpful tool to analyze your reads to determine what preprocessing is appropriate for your application.

More Answers (1)

Rik
Rik on 13 May 2020
The way I read the documentation, seqfilter selects or rejects sequences and seqtrim edits the sequences themselves.
I have never worked with this toolbox before, so I might be wrong.

Categories

Products

Release

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!