Get read depth from Multiple bam files

Hello everyone,

I need some help, I have a Panel sequencing data (107 genes) from 18000 individuals, which is already aligned it to hg38 reference genome (1 bam file for each individual), As I have to merge this data with a WES data and I do not have the region information, I was just wondering if there is a way/tool to extract the regions from those 18000 bam files into a single file that can be further used to get those same regions out of the WES data.

I would really appreciate your help/suggestions.

Regards,