[Biojava-l] dp pairwisealignment demo

Stephen Montgomery smontgom at bcgsc.bc.ca
Thu Jun 12 12:31:28 EDT 2003


hi -
i've been trying out the PairwiseAlignment demo to align to piece of 
genomic dna (after changing the alphabet from protein to dna).
i would like to view the optimum path generated by the viterbi algorithm 
but i seem to be a bit stumped as to how to extract this from the 
statepath.  it is extracted in the viterbialign demo but under a 
different hmm scenario (not a pairwise alignment).

i tried this (among other things, to no avail)

SymbolList sl = result.symbolListForLabel(StatePath.SEQUENCE);
Iterator i = sl.iterator();
while (i.hasNext()) {
Symbol sym = (Symbol) i.next();
System.out.println(sym.getName());

}

//iterate through symbol list and printing symbol names give this output

gap
(thymine cytosine)
([] [])
(adenine cytosine)
([] [])
(guanine adenine)
([] [])
(guanine guanine)
([] [])
(guanine cytosine)
([] [])
(guanine adenine)

this seems to correlate with the input sequence files.

thanks for the help.
cheers,
stephen




More information about the Biojava-l mailing list