Option to log fraction of ctx_db / regions overlap

[alexlenail]

I'm using the precomputed SCREEN regions ctx_db, but I don't know what fraction of my regions overlap with the regions in that database. This is the code that overlaps a database with a set of regions:

pycistarget/pycistarget/utils.py

Lines 282 to 288 in 7daf370

	join_pr = target_pr.join(query_pr, report_overlap = True)
	join_pr.Overlap_query = join_pr.Overlap/(join_pr.End_b - join_pr.Start_b)
	join_pr.Overlap_target = join_pr.Overlap/(join_pr.End - join_pr.Start)
	join_pr = join_pr[(join_pr.Overlap_query > fraction_overlap) | (join_pr.Overlap_target > fraction_overlap)]
	target_regions = [str(chrom) + ":" + str(start) + '-' + str(end) for chrom, start, end in zip(list(join_pr.Chromosome), list(join_pr.Start), list(join_pr.End))]
	query_regions = [str(chrom) + ":" + str(start) + '-' + str(end) for chrom, start, end in zip(list(join_pr.Chromosome), list(join_pr.Start_b), list(join_pr.End_b))]
	target_to_query = pd.DataFrame([target_regions, query_regions], index=['Target', 'Query']).T

There doesn't appear to be an option to log statistics from this. There should be! That would help me decide whether I need to create my own database.