import-bam-2: Import Bam-2
In pappewaio/AllelicImbalance: Investigates Allele Specific Expression

Description Usage Arguments Details Value Author(s) Examples

Imports bla bal bal a specified genomic region from a bam file using a GenomicRanges object as search area.

1	impBamGRL.old(UserDir, searchArea, verbose = TRUE)

`UserDir`	The relative or full path of folder containing bam files.
`searchArea`	A `GenomicRanges object` that contains the regions of interest
`verbose`	Setting `verbose=TRUE` gives details of procedure during function run.

These functions are right on tahea wrappers to import bam files into R and store them into either GRanges, GAlignments or GappedAlignmentpairs objects.

It is recommended to use the impBamGAL() which takes information of gaps into account. It is also possible to use the other variants as well, but then pre-filtering becomes important keps to understand because gapped, intron-spanning reads will cause problems. This is because the GRanges objects can not handle if gaps are present and will then give a wrong result when calculating the allele (SNP) count table.

impBamGRL returns a GRangesList object containing the RNA-seq reads in the region defined by the searchArea argument. impBamGAL returns a list with GAlignments objects containing the RNA-seq reads in the region defined by the searchArea argument. funImpBamGAPL returns a list with GappedAlignmentPairs object containing the RNA-seq reads in the region defined by the searchArea argument.

Jesper R. Gadin, Lasse Folkersen

#Declare searchArea
searchArea <- GRanges(seqnames=c('17'), ranges=IRanges(79478301,79478361))

#Relative or full path
pathToFiles <- system.file('extdata/ERP000101_subset', package='AllelicImbalance')