The haploblocks calculation based on existing phased and unphased variants is commendable but slightly unsafe and perhaps not necessary. Let's investigate if we have enough info, like haploblock ids, within the VCF or possibly by annotating with a separate BED, to skip the calculation.