How to pre-process Next Generation Sequencing data using MATLAB?
显示 更早的评论
I am trying to use MATLAB for pre-processing of NGS data. Can anyone suggest a comprehensive code for this procedure. I have tried codes suggested in this page but the codes can only be used for a limited number of tasks. For example I don't know how to filter (or mask) reads shorter than 10 nucleotides or how to treat paired-end reads. Moreover, how can I filter reads that have more than two N nucleotides? can anyone suggest a comprehensive reference for these tasks?
采纳的回答
Luuk van Oosten
2016-9-22
Dear Ehsan,
The page you refer to is a good start to get familiar with processing NGS data, but there are (a lot more!!) functions in the BioInformatics toolbox that will help you with preprocessing. Now to your specific questions:
(1) how to filter (or mask) reads shorter than 10 nucleotides:
your_filtered_data = seqfilter(yourdata.fastq, 'Method','MinLength','Threshold',10)
(2) how to treat paired-end reads
You are in luck, as there exist this thingy called 'seqsplitpe', which allows you to split merged paired-end sequences into separate files (if that is something you want).
(3) how can I filter reads that have more than two N nucleotides
This is probably a combination of (a) importing your sequences and then (b) searching your sequences for your specific repeat of >N nucleotides. I believe there are no pre-fabricated functions in MATLAB for this, but there are numerous functions which allow you to analyze sequences in the Bioinformatics Toolbox.
Best regards
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Bioinformatics Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
