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