getAlignedReads {R453Plus1Toolbox} | R Documentation |
For a given AVASet, this function imports all aligned reads belonging to all (or some selected) amplicons of all samples.
getAlignedReads(object, amplicons, dir)
object |
An instance of the |
amplicons |
An (optional) character vector of amplicon names as
mentioned in the amplicon feature data (see |
dir |
Usually, the method tries to retrieve the path to the AVA
project from the given |
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.
One DNAStringSet
that contains all aligned reads for all samples (eventually restricted to some given amplicons).
Christoph Bartenhagen
# 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)