View source: R/pairwise_alignment.R
pairwise_alignment_sequence_identity | R Documentation |
Calculate the percentage of pairwise sequence identity
pairwise_alignment_sequence_identity( seqs, aln_type = "global", pid_type = "PID1" )
seqs |
A named character vector to convert into a
|
aln_type |
A character vector of one containing the alignment type. Possible options are "global" (Needleman-Wunsch),"local" (Smith-Waterman) and "overlap". |
pid_type |
A character vector of one containing the definition of percent sequence identity. Possible options are "PID1", "PID2", "PID3" and "PID4". |
A long DataFrame with the results.
global
: align whole strings with end gap
penalties (Needleman-Wunsch).
local
: align string fragments (Smith-Waterman).
overlap
: align whole strings without end gap penalties.
PID1
: 100 * (identical positions) / (aligned positions + internal
gap positions).
PID2
: 100 * (identical positions) / (aligned positions).
PID3
: 100 * (identical positions) / (length shorter sequence).
PID4
: 100 * (identical positions) / (average length of the two
sequences).
data(phmmer_2abl) pairwise_alignment_sequence_identity( seqs = phmmer_2abl$hits.fullfasta[6:10], aln_type = "overlap", pid_type = "PID2" )
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.