What's the diff between seqfilter & seqtrim?

1 view (last 30 days)
Both of the function seems to remove sequence with bad quality, doesn't make them overlaped?

Accepted Answer

Tim DeFreitas
Tim DeFreitas on 13 May 2020
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

Find more on Bioinformatics Pipeline in Help Center and File Exchange

Products


Release

R2019b

Community Treasure Hunt

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

Start Hunting!