def get_regions_overlap_with_tfs(regions, encode_pr): 
	"""Returns boolean matrix indicating whether regions overlaps with an ENCODE TF."""
	joined = regions.pr.join(encode_pr).as_df()
	is_overlap = joined.groupby(["peak_id", "Name"]).size().unstack(-1, fill_value=0) > 0
	is_overlap = is_overlap.reindex(regions.index, fill_value=False)
	return is_overlap