I have a bed file containing mappability values of reference genome. It starts from 10,000th position and goes on. The interval for each line varies, but there are no jumps and overlaps. For example, the first line is between 10,000 and 10,0039, the second line is 10,039 and 10,040. I want to change the interval to 1000 base pairs for each line and take the mean of mappability value for that interval.
Let's say: Input chr1 10,000 10,039 0.4 ... 10,039 ..... 0.2
Output chr1 10,000 11,000 0.2 (mean mappability value between 10,000 and 11,000 position)