getAlignedReads: Import reads from an Amplicon Variant Analyzer project

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

Description

For a given AVASet, this function imports all aligned reads belonging to all (or some selected) amplicons of all samples.

Usage

1
getAlignedReads(object, amplicons, dir)

Arguments

object

An instance of the link{AVASet-class}.

amplicons

An (optional) character vector of amplicon names as mentioned in the amplicon feature data (see fDataAmp).

dir

Usually, the method tries to retrieve the path to the AVA project from the given link{AVASet} object. However, if it fails to find the directory, dir can be used to set the root directory of the AVA project.

Details

This function reports all reads for all samples together. If you want to get the reads for some samples individually, try subsetting your AVASet as in the examples below.

Value

One DNAStringSet that contains all aligned reads for all samples (eventually restricted to some given amplicons).

Author(s)

Christoph Bartenhagen

See Also

AVASet-class, fDataAmp

Examples

1
2
3
4
5
6
7
# load an AVA dataset containing 6 samples, 4 amplicons and 259 variants
data(avaSetExample)

# import all reads for amplicon "TET2_E11.04" of the first sample
avaProjectDir = system.file("extdata", "AVASet", package = "R453Plus1Toolbox")
alnReads = getAlignedReads(avaSetExample[, 1], dir=avaProjectDir, amplicons="TET2_E11.04")
show(alnReads)


Search within the R453Plus1Toolbox 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.