What's the diff between seqfilter & seqtrim?

1 次查看(过去 30 天)
Both of the function seems to remove sequence with bad quality, doesn't make them overlaped?

采纳的回答

Tim DeFreitas
Tim DeFreitas 2020-5-13
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.

更多回答(1 个)

Rik
Rik 2020-5-13
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.

类别

Help CenterFile Exchange 中查找有关 Bioinformatics Pipeline 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by