We study linear perturbations about non rotating black hole solutions in scalar-tensor theories, more specifically Horndeski theories. We consider two particular theories that admit known hairy black hole solutions. The first one, Einstein-scalar-Gauss-Bonnet theory, contains a Gauss-Bonnet term coupled to a scalar field, and its black hole solution is given as a perturbative expansion in a small parameter that measures the deviation from general relativity. The second one, known as 4-dimensional-Einstein-Gauss-Bonnet theory, can be seen as a compactification of higher-dimensional Lovelock theories and admits an exact black hole solution. We study both axial and polar perturbations about these solutions and write their equations of motion as a first-order (radial) system of differential equations, which enables us to study the asymptotic behaviours of the perturbations at infinity and at the horizon following an algorithm we developed recently. For the axial perturbations, we also obtain effective Schrödinger-like equations with explicit expressions for the potentials and the propagation speeds. We see that while the Einstein-scalar-Gauss-Bonnet solution has well-behaved perturbations, the solution of the 4-dimensional-Einstein-Gauss-Bonnet theory exhibits unusual asymptotic behaviour of its perturbations near its horizon and at infinity, which makes the definition of ingoing and outgoing modes impossible. This indicates that the dynamics of these perturbations strongly differs from the general relativity case and seems pathological.