CreateChimeras: Create Artificial Chimeras

Description Usage Arguments Details Value Author(s) See Also Examples

Description

Creates artificial random chimeras from a set of sequences.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
CreateChimeras(myDNAStringSet,
               numChimeras = 10,
               numParts = 2,
               minLength = 80,
               maxLength = Inf,
               minChimericRegionLength = 30,
               randomLengths = TRUE,
               includeParents = TRUE,
               processors = 1,
               verbose = TRUE)

Arguments

myDNAStringSet

A DNAStringSet object with aligned sequences.

numChimeras

Number of chimeras desired.

numParts

Number of chimeric parts from which to form a single chimeric sequence.

minLength

Minimum length of the complete chimeric sequence.

maxLength

Maximum length of the complete chimeric sequence.

minChimericRegionLength

Minimum length of the chimeric region of each sequence part.

randomLengths

Logical specifying whether to create random length chimeras in addition to random breakpoints.

includeParents

Whether to include the parents of each chimera in the output.

processors

The number of processors to use, or NULL to automatically detect and use all available processors.

verbose

Logical indicating whether to display progress.

Details

Forms a set of random chimeras from the input set of (typically good quality) sequences. The chimeras are created by merging random sequences at random breakpoints. These chimeras can be used for testing the accuracy of the FindChimeras or other chimera finding functions.

Value

A DNAStringSet object containing chimeras. The names of the chimeras are specified as "parent #1 name [chimeric region] (distance from parent to chimera), ...".

If includeParents = TRUE then the parents of the chimeras are included at the end of the result. The parents are trimmed to the same length as the chimera if randomLengths = TRUE. The names of the parents are specified as "parent #1 name [region] (distance to parent #2, ...)".

Author(s)

Erik Wright DECIPHER@cae.wisc.edu

See Also

FindChimeras, Seqs2DB

Examples

1
2
3
4
db <- system.file("extdata", "Bacteria_175seqs.sqlite", package="DECIPHER")
dna <- SearchDB(db)
chims <- CreateChimeras(dna)
BrowseSeqs(chims)


Search within the DECIPHER package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.