Re: [Bioc-devel] GenomicAlignments: using asMates=TRUE and yieldSize with paired-end BAM files

2014-03-27 Thread Valerie Obenchain
Hi Mike, This is fixed in Rsamtools 1.15.35. The bug was related to when the mate-pairing was performed wrt meeting the 'yieldSize' requirement. Thanks for sending the file and reproducible example. The file has ~115 million records: fl -

Re: [Bioc-devel] GenomicAlignments: using asMates=TRUE and yieldSize with paired-end BAM files

2014-03-27 Thread Valerie Obenchain
I should also mention that when both 'yieldSize' in the BamFile and 'which' in ScanBamParam are set the 'which' gets priority. The point of 'yieldSize' is to provide a reasonable chunk for processing the file. When 'which' is provided it's assumed that range is of reasonable chunk size so the

Re: [Bioc-devel] GenomicAlignments: using asMates=TRUE and yieldSize with paired-end BAM files

2014-03-27 Thread Hervé Pagès
Hi Val, On 03/27/2014 09:13 AM, Valerie Obenchain wrote: I should also mention that when both 'yieldSize' in the BamFile and 'which' in ScanBamParam are set the 'which' gets priority. The point of 'yieldSize' is to provide a reasonable chunk for processing the file. When 'which' is provided