How to obtain constraint per base?

Good afternoon and thanks for the amazing database,

Using V4.1 data, I would like to estimate the constraint at each base for many (~400 000) 1000bp-sized regions of the genome.
Here, the constraint at a position is the proportion of GnomAD genomes that exhibit a variant at the position.
I know there exists files that indicate constraint for 1000bp bins, but I couldn’t find files that indicate constraint for each base.

For now, I thought to use the GnomAD API to request data for each region, but it is pretty expensive, both in time and in internet use. Moreover, I imagine I could obtain the information using the VCF files.

Do you have a better solution ?

Thanks,

Elliot