Thanks for your reply. More precisely, I am looking for a program that would penalize the gaps on the 5' side but not on the 3' side.
It sounds to me like you want two local alignments joined by a long insertion / linker.
Exactly!
If so, setting a gap penalty to zero may get you what you want with
For my project, I can't set the penalty to zero. In fact, I am looking for gapless alignments in the 3' side and 5' side, and an insert in between. Inserts are repetitive regions, and there is a high chance that there would be a match for them in the remaining parts of the template. Actually, when I used local alignment, I set the gap opening, extension penalty to 100 to avoid any gaps.
Smith-waterman and Needleman-wunch differ in how the score of each position of the alignment table is calculated. In Needleman-wunch, the score of each position of the table is calculated from the surrounding positions. In Smith-Waterman, the score of each position is calculated based on the surrounding positions + the first position in the table. In the book entitled "Bioinformatics Algorithms: An Active Learning Approach," the authors beautifully say local alignment is like having a free ride to the start of the alignment and from the end of alignment to the final position of the scoring table.
I need the free ride for only one side, not the other side. That's why I said I want a global alignment for the 3' side and local for the 5'.
Thanks in advance for your reply,