Clear Filters
Clear Filters

FASTA AND BLAST

6 views (last 30 days)
huda nawaf
huda nawaf on 14 Nov 2011
hi, are there fasta and blast functions in matlab?
thanks

Answers (3)

Walter Roberson
Walter Roberson on 14 Nov 2011
  3 Comments
huda nawaf
huda nawaf on 15 Nov 2011
thanks
but i did not find the heuristic software BLAST for in matlab.
can help me to find it if their .
Walter Roberson
Walter Roberson on 15 Nov 2011
Well, that is odd, considering that the blast* functions in that toolbox specifically include the reference,
Altschul, S.F., Gish, W., Miller, W., Myers, E.W., and Lipman, D.J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410.

Sign in to comment.


Lucio Cetto
Lucio Cetto on 19 Nov 2011
None of these alignment algorithms are included in Matlab.
Do not mix the "fasta" alignment algorithm with the "fasta" file format. Matlab provides fastaread and fastawrite for I/O FASTA formatted files.
In the case of the "blast" algorithm Matlab does provides a wrapper that allows you to pass the information to a locally installed blast (blastlocal), and provides a read function (blastlocalread) that will let you put the results of the alignments back in Matlab, but you will need to install the blast algorithm (provided by NCBI) by yourself. Also Matlab provides a function that allows you to submit queries to the NCBI's blast server, balstncbi (i.e. aligning against their databases). The function blastread will read the web data back into Matlab.
  3 Comments
Lucio Cetto
Lucio Cetto on 21 Nov 2011
This is an interesting question that I have come accross several times.
Matlab allows you to run an alignment (using a Needleman&Wunshtype algorithm) with symbols other than aminoacids or nucleotides, for intance look at this <http://www.mathworks.com/matlabcentral/newsreader/view_thread/281532#743065 thread>. Note that this is undocummented.
Also, if you want to align instead numeric values (or even vectors) you can look at the |samplealign| function. It also uses a dynamic programming algorithm. Run the examples in the m-help, this is a function hard to understand but very handy once you understand how it works.
Note that in both of the solutions I give, none of them use an heuristic approach (such as in the BLAST or FASTA algorithms).
HTH
huda nawaf
huda nawaf on 29 Nov 2011
thanks Lucio
I did not see this last comment , walter draw my attention to this comment.
i will see that link and that function

Sign in to comment.


Jonathan
Jonathan on 14 Nov 2011
The following commands and output says there are not such functions.
>> help fasta
fasta not found.
Use the Help browser Search tab to search the documentation, or
type "help help" for help command options, such as help for methods.
>> help blast
blast not found.
Use the Help browser Search tab to search the documentation, or
type "help help" for help command options, such as help for methods.
If you can tell us the purpose of the functions, perhaps we can answer if there are equivalent functions in Matlab.
  4 Comments
Walter Roberson
Walter Roberson on 26 Dec 2011
Are you still trying to work with non-biological data? If so then be aware that the Smith-Waterman algorithm is "Identification of Common Molecular Subsequences" by "T. F. SMITH AND M. S. WATERMAN", and is thus inherently about biological data.
It is possible that the techniques used by Smith-Waterman could be adapted for non-biological data, but then it would not be the Smith-Waterman algorithm.
I note that Wikipedia http://en.wikipedia.org/wiki/Smith%E2%80%93Waterman_algorithm claims that the Smith-Waterman algorithm in its original form is not often used now, partly because it is slow.
Diego
Diego on 28 Dec 2011
Did you already look at the swalign.m file located inside the bioinfo/bioinfo folder?.

Sign in to comment.

Tags

Community Treasure Hunt

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

Start Hunting!