Extracting exon neighbourhoods from the human reference genome

Hi all,

Starting from the GRCh38/hg38 human reference genome, I would like to scan through all the exons with an aim to extract all the neighbourhoods exon–intro–exon–intro–exon

I would later like to do some operations e.g. calculating lengths and doing some sequence predictions

Anyone familiar what would be the best starting point to get these exon neighbourhoods extracted?

Thanks,
Sergio