multipleSeqAlign: Compare multiple sequence alignments

In quick2063706271/Cov2Comparator: Analyzing SARS covid-2 genome across different geographic regions

Description

A function that calculates multiple sequence alignments

Usage

multipleSeqAlign(sequences, algorithm = "clustalw")

Arguments

sequences	A DNAStringset containing several sequences
algorithm	A string indicating the algorithm that user wants to use to calculate multiple sequence alignment

Value

Returns a MsaAAMultipleAlignment object which is the result of multiple sequence alignment

References

Charif D, Lobry J. 2007. “SeqinR 1.0-2: a contributed package to the R project for statistical computing devoted to biological sequences retrieval and analysis.” In Bastolla U, Porto M, Roman H, Vendruscolo M (eds.), Structural approaches to sequence evolution: Molecules, networks, populations, series Biological and Medical Physics, Biomedical Engineering, 207-232. Springer Verlag, New York.

Paradis E. & Schliep K. 2019. ape 5.0: an environment for modern phylogenetics and evolutionaryanalyses in R. Bioinformatics 35: 526-528.

U. Bodenhofer, E. Bonatesta, C. Horejs-Kainrath, and S. Hochreiter (2015) msa: an R package for multiple sequence alignment. Bioinformatics 31(24):3997- 9999. DOI: 10.1093/bioinformatics/btv176.

Examples

# Example 1
# Create a basic msa and then plot a tree from it
library(Biostrings)
set1 <- Biostrings::DNAStringSet("ATCGATCG")
set2 <- Biostrings::DNAStringSet("ATTTTTTT")
set3 <- Biostrings::DNAStringSet("ATCGATTT")
setTotal <- union(set1, set2)
setTotal <- union(setTotal, set3)
align <- multipleSeqAlign(setTotal)

quick2063706271/Cov2Comparator documentation built on Dec. 22, 2021, 10:59 a.m.