I imagine this is an easy question, but I haven't been able to find the answer! I have a bed file list of all HG19 genes. I also have BigWig files of MACS2 peak calls from ChIP-seq data. I would like to iterate through my list of genes, and assign the "score" value of the bed file as the signal from the BigWig at the bed file's indicated "start" value.
I have used DeepTools to achieve a similar goal, but the computeMatrix output, "sorted/filtered regions," does not have a true "score".