DNA Sequence Alignment Scoring With Extended Gap and Similarity Matrix

DNA Sequence Alignment Scoring With Extended Gap and Similarity Matrix

0

Hey everyone! I am having a little difficulty trying to figure out how to align two sequences in Python using a Dynamic Programming algorithm. I’ve tried to take a look at the source code for jAligner, but I’m having issues understanding how it works.

Basically the problem is this:

Given two sequences, a scoring matrix given that gives scores between each nucleotide, a start gap penalty, and a continue gap penalty.
For example:

   A   G   C   T
A  1  -2  -2  -1
G -2   1  -1  -2
C -2  -1   1  -5
T -1  -2  -2   1

openGap = -3
continueGap = -1

How do we compute the alignment of two sequences A and B (with length N and M) such that the score is maximized?

I’ve taken a look at This Wikipedia Page, but it doesn’t mention anything about a continue/open gap, rather just one single gap score.


programming


dynamic


aligning


needleman-wunsch

• 19 views

Read more here: Source link