Bayesian sampling of genomic rearrangement scenarios via double cut and join