Feature suggestion as discussed via email with credits to @asntech .
It seems from looking at the source code that the genomic regions overlap is calculated using pybedtools intersect default settings which counts as intersection any regions that overlap by at least 1 bp. Could please consider
adding a note to the Docs, and/or
consider adding an option to, for instance, consider an overlap only regions matching 100 %? The perfect solution would be to an option, where user will be able to list any of the arguments available for bedtool’s intersect command. Something like this: --bedtool-options f=0.8,r=true,etc
The default is perfectly sensible, but it might lead to unexpected results for those not aware of it.
Feature suggestion as discussed via email with credits to @asntech .
It seems from looking at the source code that the genomic regions overlap is calculated using
pybedtools intersect
default settings which counts as intersection any regions that overlap by at least 1 bp. Could please considerThe default is perfectly sensible, but it might lead to unexpected results for those not aware of it.