sosorum variant calling and theta / h