Description Usage Arguments Details Value Author(s) References Examples
Executes Clustal on a set of sequences to obtain a multiple sequence alignment.
1 2 | clustal(x, param)
clustal_help()
|
x |
an object of class XStringSet (e.g., DNAStringSet) with the sequences to be aligned. |
param |
character string with the command line parameters for clustal (see output of |
The clustalw software has to be installed (see rMSA_INSTALL).
An object of class DNAMultipleAlignment (see BioStrings).
Michael Hahsler
Larkin M., et al. Clustal W and Clustal X version 2.0, Bioinformatics 2007 23(21):2947-29
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 | ## Not run:
### DNA
dna <- readDNAStringSet(system.file("examples/DNA_example.fasta",
package="rMSA"))
dna
al <- clustal(dna)
al
### inspect alignment
detail(al)
### plot a sequence logo for the first 20 positions
plot(al, 1, 20)
### RNA
rna <- readRNAStringSet(system.file("examples/RNA_example.fasta",
package="rMSA"))
rna
al <- clustal(rna)
al
### Proteins
aa <- readAAStringSet(system.file("examples/Protein_example.fasta",
package="rMSA"))
aa
al <- clustal(aa)
al
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.