copy_detect {copyseparator}

copy_detect

Description

Separates two or more gene copies from a single subset of short reads.

Usage

copy_detect(filename, copy_number, verbose = 1)

Arguments

filename

A fasta file contains short reads from a single subset generated by "subset_downsize".

copy_number

An integer (e.g. 2,3, or 4) giving the anticipated number of gene copies in the input file.

verbose

Turn on (verbose=1; default) or turn off (verbose=0) the output.

Value

A fasta alignment of the anticipated number of gene copies.

Examples

## Not run: 
copy_detect("inst/extdata/toysubset.fasta",2,1)

## End(Not run)