Hi, say I'm working with some RNAseq reads
What I'm usually doing is taking coverage() and giving it a BAM file and a GRange for the param and I assume this will calculate coverage for the segment that I specify in the GRange.
Would it be possible to get the reads for the segment with maybe scanBam even the ones that aren't completely inside the scanned segment. Take a look at their positions and remove/shift a few based on position. Then run coverage() based on the altered set of reads? Preferably without generating new large bam files ?
(although if this is the only way its acceptable although I'm pretty sure its not necessary)
Thanks.