Comparing the genomes from closely related species provides a powerful tool to identify functional elements in a reference genome. Many methods have been developed to identify conserved sequences across species, however, existing methods only model conservation as a decrease in the rate of mutation and have ignored selection acting on the pattern of mutations.