consensus | R Documentation |
consensus of aligned sequences
consensus of aligned sequences
consensus(x, type = "DNA")
bs_consensus(x, type = "DNA", r = 1)
x |
BStringSet object |
type |
currently, only DNA supported |
r |
if any NT > r, it will be selected as representative base |
consensus sequence string
consensus sequence string
Guangchuang Yu
## Not run:
fa_file <- system.file("extdata/HA.fas", package="seqmagick")
x <- fa_read(fa_file)
## align first 5 sequences, use `bs_aln(x)` to align all sequences
aln <- bs_aln(x[1:5])
## or bs_consensus(aln)
consensus(aln)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.