Get the names of the peaks close to a gene promoter

Get the names of the peaks close to a gene promoter

Usage

get_promoter_peaks(
  peaks,
  gene,
  max_dist_to_promoter_peak = 500,
  tss_intervals = "intervs.global.tss"
)

Arguments

peaks

a PeakIntervals object

gene

a gene name, e.g. "CD4"

max_dist_to_promoter_peak

how far from gene's TSS to search for promoter-proximal peaks. Default: 500

tss_intervals

name of the intervals set containing the TSSs

Value

names of the peaks close to gene's promoter. If no peaks are found, NULL is returned.

Examples

if (FALSE) {
get_promoter_peaks(mc_atac@peaks, "CD4")
}