Rsamtools
and extract the reads for each chromosome.
bam2gr(bamFile, chr=NULL, PE=FALSE, verbose=FALSE)
logical
. This should be set to TRUE for paired-end sequencing data.logical
. Print additional information about the data.GRanges
of all the reads for each chromosome.
segmentPICS