[Bioc-devel] GRangesList Conversion Fails For Unstranded Sequencing Data

Dario Strbenac dstr7320 at uni.sydney.edu.au
Tue Mar 14 07:00:10 CET 2017

Good day,

countOverlaps doesn't work for a GAlignmentPairs object with strandMode set to 0. This is because of an oversight in the grglist function. It has an if statement that checks whether the strand mode is 1 or 2. Then, it tries to subset the variable 'x_unlisted'. However, if strand mode is 0, neither of the conditional sections of code are executed and Error in .local(x, use.names, use.mcols, ...) : object 'x_unlisted' not found happens because the 'x_unlisted' variable has not been created. It's a surprise no one else has encountered this bug before.

Dario Strbenac
University of Sydney
Camperdown NSW 2050

More information about the Bioc-devel mailing list