Skip to content

Fixing traceback of global alignments - #1

Open
tbertelsen wants to merge 1 commit into
ahmedmoustafa:masterfrom
tbertelsen:master
Open

Fixing traceback of global alignments#1
tbertelsen wants to merge 1 commit into
ahmedmoustafa:masterfrom
tbertelsen:master

Conversation

@tbertelsen

Copy link
Copy Markdown

The traceback for global alignments is not created correctly.

Test case

(source):

>F12Fcsw_257171
TTGGGCCGTGTCTCAGTCCCAATGTGGCCGATCACCCTCTCAGGTCGGCTACCGATCGTCGGCTTGGTGGGCCGTTACCT
CACCAACTACCTAATCGGACGCGAGCCCACCCCAAACCGATAATTCTTTTACCCCAGAACCATGTGATCCCGTGGTCTTA
TGCGGTATTAGTACACCTTTCGGTGTGTTATTCCCTCGTCTGGGAAAGGGTTAGTCTCACGCGTTACTCCACCCGTCCCG
CCCGCCTAAAACAAAGCTCTAA
>M13Fcsw_294419
TTGGGCCGTGTCTCAGTCCCAATGTGGCCGATCACCCTCTCAGGTCGGCTACCGATCGTCGGCTTGGTGGGCCGTTACCT
CACCAACTACCTAATCGGACGCGAGCCCACCCCAAACCGATAAATCTTTTACCTCAGAACCATGTGATCCCGTGGTCTTA
TGCGGTATTAGTACACCTTTCGGTGTGTTATTCCCCTGTCTGGGAAAGGTTGCTCACGCGTTACTCACCCGTCCGCCGCT
AAAACAGCT

Alignment before change:

Alignment.getScore() = 175
Alignment.calculateScore() = 131

println(new String(alignment.getSequence1) + "\n" new String(alignment.getMarkupLine) + "\n" new String(alignment.getSequence2))

TTGGGCCGTGTCTCAGTCCCAATGTGGCCGATCACCCTCTCAGGTCGGCTACCGATCGTCGGCTTGGTGGGCCGTTACCTCACCAACTACCTAATCGGACGCGAGCCCACCCCAAACCGATAATTCTTTTACCCCAGAACCATGTGATCCCGTGGTCTTATGCGGTATTAGTACACCTTTCGGTGTGTTATTCCCTCGTCTGGGAAAGGGTTAGTCTCACGCGTTACTCCACCCGTCCCGCCCGCCTAAAACAAAGCTCTAA
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||.|||||||||.|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||..||||||||||||.|.............|...||..|||.|.|.....|....             
TTGGGCCGTGTCTCAGTCCCAATGTGGCCGATCACCCTCTCAGGTCGGCTACCGATCGTCGGCTTGGTGGGCCGTTACCTCACCAACTACCTAATCGGACGCGAGCCCACCCCAAACCGATAAATCTTTTACCTCAGAACCATGTGATCCCGTGGTCTTATGCGGTATTAGTACACCTTTCGGTGTGTTATTCCCCTGTCTGGGAAAGGTTGCTCACGCGTTACTCACCCGTCCGCCGCTAAAACAGCT-------------

Alignment after change:

Alignment.getScore() = 175
Alignment.calculateScore() = 175

println(new String(alignment.getSequence1) + "\n" new String(alignment.getMarkupLine) + "\n" new String(alignment.getSequence2))

TTGGGCCGTGTCTCAGTCCCAATGTGGCCGATCACCCTCTCAGGTCGGCTACCGATCGTCGGCTTGGTGGGCCGTTACCTCACCAACTACCTAATCGGACGCGAGCCCACCCCAAACCGATAATTCTTTTACCCCAGAACCATGTGATCCCGTGGTCTTATGCGGTATTAGTACACCTTTCGGTGTGTTATTCCCTCGTCTGGGAAAGGGTTAGTCTCACGCGTTACTCCACCCGTCCCGCCCGCCTAAAACAAAGCTCTAA
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||.|||||||||.|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||..||||||||||||.|   .||||||||||||||      .||||.|||||...||.|.||||    
TTGGGCCGTGTCTCAGTCCCAATGTGGCCGATCACCCTCTCAGGTCGGCTACCGATCGTCGGCTTGGTGGGCCGTTACCTCACCAACTACCTAATCGGACGCGAGCCCACCCCAAACCGATAAATCTTTTACCTCAGAACCATGTGATCCCGTGGTCTTATGCGGTATTAGTACACCTTTCGGTGTGTTATTCCCCTGTCTGGGAAAGGTT---GCTCACGCGTTACTC------ACCCGTCCGCCGCTAAAACAGCT----

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant