Active flow-control techniques have shown promise for achieving high levels of drag reduction. However, these techniques are often complex and involve multiple tunable parameters, making it challenging to optimize their efficiency. Here, we present a Bayesian optimization (BO) approach based on Gaussian process regression (GPR) to optimize a wall-normal blowing and suction control scheme for a NACA 4412 wing profile at two angles of attack: 5 and 11 degrees, corresponding to cruise and high-lift scenarios, respectively. An automated framework is developed by linking the BO code to the CFD solver OpenFOAM. RANS simulations (validated against high-fidelity LES and experimental data) are used in order to evaluate the different flow cases. BO is shown to provide rapid convergence towards a global maximum, even when the complexity of the response function is increased by introducing a model for the cost of the flow control actuation. The importance of considering the actuation cost is highlighted: while some cases yield a net drag reduction (NDR), they may result in an overall power increase. Furthermore, optimizing for NDR or net power reduction (NPR) can lead to significantly different actuation strategies. Finally, by considering losses and efficiencies representative of real-world applications, still a significant NPR is achieved in the 11° case, while net power reduction is only marginally positive in the 5° case.