generate segregating sites, convert the output gene tree branch length to 4Ne