clustal: Run Multiple Sequence Alignment (ClustalW) on a Set of...

In mhahsler/rMSA: Interface for Popular Multiple Sequence Alignment Tools

Run Multiple Sequence Alignment (ClustalW) on a Set of Sequences

Description

Executes Clustal on a set of sequences to obtain a multiple sequence alignment.

Usage

clustal(x, param)
clustal_help()

Arguments

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 clustal_help()).

Details

For installation details see: https://github.com/mhahsler/rMSA/blob/master/INSTALL

Value

An object of class DNAMultipleAlignment (see BioStrings).

Author(s)

Michael Hahsler

References

Larkin M., et al. Clustal W and Clustal X version 2.0, Bioinformatics 2007 23(21):2947-29

Examples

## 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)

mhahsler/rMSA documentation built on May 24, 2024, 3:36 p.m.