seqdiff: seqdiff

Description Usage Arguments Value Author(s) Examples

View source: R/seqdiff.R

Description

calculate difference of two aligned sequences

Usage

1
seqdiff(fasta, reference = 1)

Arguments

fasta

fasta file

reference

which sequence serve as reference, 1 or 2

Value

SeqDiff object

Author(s)

guangchuang yu

Examples

1
2
fas <- list.files(system.file("examples","GVariation", package="seqcombo"), pattern="fas", full.names=TRUE)
seqdiff(fas[1], reference=1)

Example output

sequence differences of Mont and CF_YL21 
1181 sites differ:
  A   C   G   T 
286 315 301 279 

seqcombo documentation built on Nov. 8, 2020, 6:27 p.m.