What genome are you working with ? For any genome, if you already have genes in fasta, there must be genomic coordinates for those genes in either GTF/GFF or bed format. If you have GTF, convert it to a bed file of exons or gene start end and then calculate coverage using bedtools. Either you can calculate per base pair or per region ( for each gene)