What's the diff between seqfilter & seqtrim?
1 view (last 30 days)
Show older comments
Both of the function seems to remove sequence with bad quality, doesn't make them overlaped?
0 Comments
Accepted Answer
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.
0 Comments
More Answers (1)
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.
0 Comments
See Also
Categories
Find more on Bioinformatics Pipeline in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!