[BioC] readGAlignmentPairsFromBam with yieldSize & which

[Stefano Calza]

Hi!

I am experiencing this weird (to me...) behaviour:

bfl = BamFile(some/bam/file,yieldSize=100000L)

param=ScanBamParam(which=RangesList(chr22=IRanges(1,51304566)))

## 1
readGAlignmentsFromBam(bfl)

## 2

readGAlignmentsFromBam(bfl,param=param)


Now 1) is much faster than 2). Indeed 1) takes ~ 6 secs while 2) at 
least a couple of minutes till I killed it...

Am I doing something wrong or missing something here?

Stefano