Methods for computing genome-wide four-gamete compatible intervals describing local phylogenetic structure