[Biojava-l] How to parse pairwise alignment result string
David Zhao
wzhao6898 at gmail.com
Thu Jul 9 20:22:47 UTC 2009
Hi there,
I've successfully aligned 2 peptide sequence using Needleman-Wunsch algorithm:
SequenceAlignment aligner.pairwiseAlignment(
query, // first sequence
target // second one
);
Now, the only output I can get is this string returned by
aligner.getAlignmentString();
showing below:
Time (ms): 63
Length: 25
Score: 180
Query: query, Length: 25
Target: target, Length: 25
Query: 1 KLFVGGIKEDTEEHHLRDYFEEYGK 25
| ||||||||||||||||||| |||
Target: 1 KIFVGGIKEDTEEHHLRDYFEQYGK 25
How can create an gapped alignment object in biojava from this,
so I can retrieve
score, gap information, etc. from the object?
Thanks in advance!
David
More information about the Biojava-l
mailing list