[Biojava-l] Pairwise alignment of protein sequences
David Zhao
wzhao6898 at gmail.com
Thu Jul 9 18:13:45 UTC 2009
David Zhao <wzhao6898 <at> gmail.com> writes:
I found the answer from a post here(http://portal.open-
bio.org/pipermail/biojava-l/2007-February.txt#). In short, I need to use
AlphabetManager.alphabetForName("PROTEIN-TERM"); instead of
AlphabetManager.alphabetForName("PROTEIN"); to parse the matrix file.
Another question though, now I have the alignment,
how do I retrieve score and
mismatch, gap information from it?
Time (ms): 62
Length: 25
Score: 180
Query: query, Length: 25
Target: target, Length: 25
Query: 1 KLFVGGIKEDTEEHHLRDYFEEYGK 25
| ||||||||||||||||||| |||
Target: 1 KIFVGGIKEDTEEHHLRDYFEQYGK 25
Thanks!
David
More information about the Biojava-l
mailing list