Converts a FASTA formatted multiple alignment to a matrix.
Usage
msa2mat(msa)
Value
A matrix where each row is a vector of aligned bases/amino acids.
Arguments
msa
A fasta object with a multiple alignment, see msalign`.
Author
Lars Snipen.
Details
This function converts the fasta object msa, containing a multiple alignment,
to a matrix. This means each position in the alignment is a column in the matrix, and the
content of the Header column of msa is used as rownames of theh matrix.
Such a matrix is useful for conversion to a DNAbin object that is used by the ape
package for re-constructing phylogenetic trees.
msa.file <- file.path(path.package("microseq"),"extdata", "small.msa")
msa <- readFasta(msa.file)
msa.mat <- msa2mat(msa) # to use with ape::as.DNAbin(msa.mat)