Massive read annotation

[rebrendi]

Hello,

I want to annotate all my reads (~1 billion) with respect to genomic features such as transcription start sites. I just want my old BED file with mapped reads to be modified in such a way that a new column would be added, containing the distance of a given read from the nearest genomic feature. Any ideas which program or a publicly available script can do this?

Thanks!