skbio.alignment.
global_pairwise_align_protein
(seq1, seq2, gap_open_penalty=11, gap_extend_penalty=1, substitution_matrix=None, penalize_terminal_gaps=False)[source]¶Globally align pair of protein seqs or alignments with Needleman-Wunsch
State: Experimental as of 0.4.0.
Parameters: | seq1 : Protein or TabularMSA[Protein]
seq2 : Protein or TabularMSA[Protein]
gap_open_penalty : int or float, optional
gap_extend_penalty : int or float, optional
substitution_matrix: 2D dict (or similar), optional
penalize_terminal_gaps: bool, optional
|
---|---|
Returns: | tuple
|
See also
local_pairwise_align
, local_pairwise_align_protein
, local_pairwise_align_nucleotide
, skbio.alignment.local_pairwise_align_ssw
, global_pairwise_align
, global_pairwise_align_nucelotide
Notes
Default gap_open_penalty
and gap_extend_penalty
parameters are
derived from the NCBI BLAST Server [R205206].
The BLOSUM (blocks substitution matrices) amino acid substitution matrices were originally defined in [R206206].
This function can be use to align either a pair of sequences, a pair of alignments, or a sequence and an alignment.
References
[R205206] | (1, 2) http://blast.ncbi.nlm.nih.gov/Blast.cgi |
[R206206] | (1, 2) Amino acid substitution matrices from protein blocks. S Henikoff and J G Henikoff. Proc Natl Acad Sci U S A. Nov 15, 1992; 89(22): 10915-10919. |