mafft: Run Multiple Sequence Alignment (MAFFT) on a Set of Sequences

Description Usage Arguments Details Value Author(s) References Examples

Description

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

Usage

1
2
mafft(x, param="--auto")
mafft_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 (see output of mafft_help()).

Details

The MAFFT software has to be installed (see rMSA_INSTALL).

Value

An object of class DNAMultipleAlignment (see BioStrings).

Author(s)

Michael Hahsler

References

Katoh, Standley 2013 (Molecular Biology and Evolution 30:772-780) MAFFT multiple sequence alignment software version 7: improvements in performance and usability.

Examples

 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 <- mafft(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 <- mafft(rna)
al

### Proteins
aa <- readAAStringSet(system.file("examples/Protein_example.fasta",
	package="rMSA"))
aa

al <- mafft(aa)
al

## End(Not run)

mhahsler/rMSA documentation built on May 22, 2019, 8:55 p.m.