Entering edit mode
6.9 years ago
benanfel
▴
10
hey guys i need your help about something,i need to write a programm using those instructions:
input:
1Two protein sequences
2. Algorithm chosen: Needleman-Wunsch or Smith-Waterman (both algorithms must be implemented and it is up to the user to choose the execution of one of the two)
3. Substitution matrix: PAM250 or BLOSUM62
4. Constant penalty: delta
output:
All optimal alignments with the score obtained:
In order to obtain the optimal alignments and the corresponding score, the pairwise2 module can be used as follows (local or global alignment):
Module Import:
From Bio import pairwise2
Call functions:
Alignments = pairwise2.align.globalds (seq1, seq2, matrix, delta, delta)
Alignments = pairwise2.align.localds (seq1, seq2, matrix, delta, delta)