Simple program to align sequences in Julia
This code provides two main functions, smithwaterman(seq1::String, seq2::String
, and needlemanwunsch(seq1::String, seq2::String)
(respectively local sequence alignment, global sequence alignment).
Allowed seq1,seq2
are just for DNA,RNA sequences (ACGT,acgt,ACGU,
acgu, or mixed capitalization). Score matrix is built in. So far the
program compute equivalent optimal solutions only for local alignments
(smithwaterman) but not yet for global alignment (needlemanwunsh).
The output is Array{(String,String,String),1}
where the first two
strings represent the alignment, and the third a vector of symbols:
+ = match, o = not match, - = insert
Just type using SeqAlign
at julia prompt
Protein sequences, custom score matrices.