Skin-friction drag is a major engineering concern, with wide-ranging consequences across many industries. Active flow-control techniques targeted at minimising skin friction have the potential to significantly enhance aerodynamic efficiency, reduce operating costs, and assist in meeting emission targets. However, they are difficult to design and optimise. Furthermore, any performance benefits must be balanced against the input power required to drive the control. Bayesian optimisation is a technique that is ideally suited to problems with a moderate number of input dimensions and where the objective function is expensive to evaluate, such as with high-fidelity computational fluid dynamics simulations. In light of this, this work investigates the potential of low-intensity wall-normal blowing as a skin-friction drag reduction strategy for turbulent boundary layers by combining a high-order flow solver () with a Bayesian optimisation framework. The optimisation campaign focuses on streamwise-varying wall-normal blowing, parameterised by a cubic spline. The inputs to be optimised are the amplitudes of the spline control points, whereas the objective function is the net-energy saving (NES), which accounts for both the skin-friction drag reduction and the input power required to drive the control (with the input power estimated from real-world data). The results of the optimisation campaign are mixed, with significant drag reduction reported but no improvement over the canonical case in terms of NES. Selected cases are chosen for further analysis and the drag reduction mechanisms and flow physics are highlighted. The results demonstrate that low-intensity wall-normal blowing is an effective strategy for skin-friction drag reduction and that Bayesian optimisation is an effective tool for optimising such strategies. Furthermore, the results show that even a minor improvement in the blowing efficiency of the device used in the present work will lead to meaningful NES.