noteworthy_seqs {BALCONY}R Documentation

Find noteworthy sequences in the dataset (aligned sequences)

Description

This function detects noteworthy sequences (most common, closest to the consensus and most different from the consesus) to facilitate convenient detection of outlying sequences that might be excluded from the further analysis.

Usage

noteworthy_seqs(percentage, alignment)

Arguments

percentage

The identity of each sequence in the alignment to the consensus sequence. Output of the cons2seqs_ident function

alignment

Alignment loaded with read.alignment function

Value

best_consensus

Sequence closest to the consensus

worst_consensus

Sequence most different to the consensus

most_common

Most common sequence in the alignment

Author(s)

Alicja Pluciennik & Michal Stolarczyk

Examples

data("alignment")
consensus_seq = consensus(alignment, 30)
consensus_to_sequences_identity=cons2seqs_ident(alignment,consensus_seq)
noteworthy_seqs(consensus_to_sequences_identity, alignment)

[Package BALCONY version 0.2.10 Index]