Regionfilter not filtering correct for one base regions
Regions of size 1 (closed interval [x,x]) are not included as regions because of a query filter for start < end
. This should be start <= end
.
Consequence: Variants that fall in intervals like this will be erronously filtered out because the region is not considered. Example region: NM_145061.5, cds_start=21729289, exon_end=21729290 should create a region [21729289, 21729289], but does not.
It will also consistenly fail to include splice regions or UTR regions if filter_config is e.g. {"splice_region": [-1,1], "utr_region": [-1,1]}
. There should be tests for these cases.
Edited by Øyvind Evju