BALCONY (version 0.2.10)

consensus: Consensus sequence determination

Description

Function calculates consensus sequence for given alignment with a threshold of user's choice.

Usage

consensus(alignment, threshold)

Arguments

alignment

output of of read.alignment function or grouped alignment created with: align_seq_mtx2grs and alignment2matrix

threshold

minimal fraction of amino acids on the certain position in all sequences of the alignment to be taken for consensus letter on this position; number in range 0-100.

Value

consensus_sequence

A character vector of length of the aligned sequence containing consesus sequence based on the input alignment

Details

If maximum fraction of any amino acid on the certain position is lower than a threshold then "*" is printed instead.

See Also

s2c

Examples

# NOT RUN {
data("alignment")
alignment = delete_isoforms(alignment)
threshold=80 # Set the consensus threshold
consensus_sequence=consensus(alignment, threshold)
# }