How to execute smith waterman algorithm in matlab?
13 次查看(过去 30 天)
显示 更早的评论
I want to get an optimal sequence for given set of sequences. So, I am using Smith Water-man algorithm. How to execute that algorithm in Matlab?
0 个评论
采纳的回答
Walter Roberson
2015-10-16
See also http://www.mathworks.com/matlabcentral/fileexchange/45831-matlab-audio-analysis-library/content/library/smithWaterman.m if you want code for the algorithm itself.
更多回答(1 个)
Vanessa Srebny
2019-4-29
Hej,
thanks for the post, I was trying to use your function to split my fasta file, however it only returns a file fasta.txt.1 nothing more. It is a big file more than 200000 characters it should be more than 1 resulting file right? I would like to use the Smith -Waterman function to compare a whole genome sequence to a protein sequence. When without splitting there are to many characters for matlab.
I attached my to 2 sequences:
Thank you very much :)
Vanessa
1 个评论
Walter Roberson
2019-4-30
? Which function were you using to split your fasta file?
When I use fastaread() to read each of your files, and use swalign() to align them, I am able to do that provided that I do not ask for the undocumented fourth output of swalign. I can also nwalign() them, which gives a different answer (with larger gaps).
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Genomics and Next Generation Sequencing 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!