noteworthy_seqs: Find noteworthy sequences in the dataset (aligned sequences)
In michalstolarczyk/BALCONY: Better ALignment CONsensus analYsis

Description Usage Arguments Value Examples

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.

1	noteworthy_seqs(percentage, alignment)

`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

`best_consensus`	Sequence closest to the consensus
`worst_consensus`	Sequence most different to the consensus
`most_common`	Most common sequence in the alignment

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