that's what I was trying to get at talking about "changes of coordinates for trajectories". I'm being imprecise because I'm not sure how to formalise this. On the original space you cannot simply have a coordinate transformation $q\mapsto f(q)$, as you point out (I mean, not all phase space transformations correspond to this). It must be something somehow encoding the information in $\dot q\mapsto (...)$.
But then again, if you consider the equations on-shell, the momenta are directly related to the coordinates via $p=\dot q$, at least in simple cases, right? After doing a canonical transf…