Hi Biostars I'm trying to find out to what level on a phylogenetic tree a regulatory sequence is conserved. I take known regulatory sites from E. coli and try to find where this region was gained in evolution. My question is how do I do it in an elegant way rather than just extracting the sequences, deciding which are similar enough to the one of E. coli and reconstructing the event on the tree to get the gain node. I found the program dless from phastCons that should do that but this program finds a conserved sequence that was gained not necessarily in my clade of interest. Do you have any idea? Thanks!