Question on Hap.py

[elowy01]

Dear Hap.py authors,

First, thanks a lot for Hap.py, it is an excellent tool.

I've a question on the possibility of using Hap.py in order to assess the accuracy of the phasing (haplotype blocks) of a single sample phased VCF query using a truth phased VCF (for example from GIAB) . Is that possible?

I've read the manual and I understand that it can normalize the variants and calculate the precision and recall of your query VCF, but it is not so clear to me if Hap.py can be used to check if the VCF is correctly phased.

Thanks,

ernesto