Hi, I finally solved my problem, the above mentioned file, which I have from https://console.cloud.google.com/storage/browser/genomics-public-data/resources/broad/hg38/v0 has a problem with the chromosome lengths in the vcf header. When I run CalculateGenotypePosteriors with -supporting 1000G.phase3.integrated.sites_only.no_MATCHED_REV.hg38.vcf.gz it gives me an error:

A USER ERROR has occurred: Input files reference and features have incompatible contigs: Found contigs with the same name but different lengths: contig reference = chr15 / 101991189 contig features = chr15 / 90338345

Inspecton of the file reveals that chr15 length has been replace by chr16 length and that chr17 length has been used twice for chr16 and chr17 length

contig=<ID=chr15,assembly=GCF_000001405.26,length=90338345>

contig=<ID=chr16,assembly=GCF_000001405.26,length=83257441>

contig=<ID=chr17,assembly=GCF_000001405.26,length=83257441>

I replace the lengths of chr15 and chr16 with the correct ones (which I have from the vcf file of my aligned data), resulting in these lines

contig=<ID=chr15,assembly=GCF_000001405.26,length=101991189>

contig=<ID=chr16,assembly=GCF_000001405.26,length=90338345>

contig=<ID=chr17,assembly=GCF_000001405.26,length=83257441>

Now I can run CalculateGenotypePosteriors without problems.

To add to that: What I did not realize was that the file 1000G.phase3.integrated.sites_only.no_MATCHED_REV.hg38.vcf.gz only has variants until chr15, the rest is missing. So while the program can be run with the above modifications of the file, it would probably be affected by the missing data > chr15... So it is basically back to square one....