Modeling regulatory sites with higher order position-dependent weight matrices