Entering edit mode
I used the pairwiseAlignment
from the Biostrings
package to get a pairwise alignment of two protein sequences and pid
function to get the percentage sequence identity. How can I get the percentage similarity and percentage gap of the alignment?
Also, what is the difference in using pairwise alignment of two proteins like this from Biostrings and blastp
from NCBI?