Genome-wide heterozygosity

Heys,

I want to calculate genome-wide heterozygosity for a couple of samples I have. I have been looking in internet but I just find theoretical definitions of how you should do it, rather than practical examples. Is there any software or custom batch script already optimized? I have bam or vcf files!

Thanks in advance!