globalAligmentProteinSequence: Finds the optimal global nucleotide sequence aligment with...

In frosinastojanovska/Bioinformatics: A collection of functions for Bioinformatics

Description

Finds the optimal global nucleotide sequence aligment with Needleman-Wunsch algorithm with BLOSUM50 as score matrix.

Usage

globalAligmentProteinSequence(sequence1, sequence2, gapPenalty, blosumName)

Arguments

sequence1	First vector with characters indicating some protein sequence.
sequence2	Second vector with characters indicating some protein sequence.
gapPenalty	A number indicating the gap penalty value.
blosumName	Name of the BLOSUM matrix.

Value

List containing the score, the resulting aligned sequences and the matrix of aligment.

Examples

globalAligmentProteinSequence(c("P","A","W","H","E","A","E"), c("H","E","A","G","A","W","G","H","E","E"), 8, "BLOSUM50")

frosinastojanovska/Bioinformatics documentation built on May 16, 2019, 3:32 p.m.