ggmsa (version 1.0.2)

seqdiff: seqdiff

Description

calculate difference of two aligned sequences

Usage

seqdiff(fasta, reference = 1)

Arguments

fasta

fasta file

reference

which sequence serve as reference, 1 or 2

Value

SeqDiff object

Examples

# NOT RUN {
fas <- list.files(system.file("extdata", "GVariation", package="ggmsa"),
                  pattern="fas", full.names=TRUE)
seqdiff(fas[1], reference=1)
# }