extension to perform a double stage maximum parsimony -- peanlty
depended on changes of number of domain but also on the number of total
genes