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: |
|
---|---|
Returns: |
|
Return type: | 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 [1].
The BLOSUM (blocks substitution matrices) amino acid substitution matrices were originally defined in [2].
This function can be use to align either a pair of sequences, a pair of alignments, or a sequence and an alignment.
References
[1] | http://blast.ncbi.nlm.nih.gov/Blast.cgi |
[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. |