This work aims to explore the solution of a nonlinear fractional integro-differential equation in the complex domain through the utilization of both analytical and numerical approaches. The demonstration of the existence and uniqueness of a solution is established under certain appropriate conditions with the use of Banach fixed point theorems. To date, no research effort has been undertaken to look into the solution of this integro equation, particularly due to its fractional order specification within the complex plane. The validation of the proposed methodology was performed by utilizing a novel strategy that involves implementing the Rationalized Haar wavelet numerical method with the application of the Bernoulli polynomial technique. The primary reason for choosing the proposed technique lies in its ability to transform the solution of the given nonlinear fractional integro-differential equation into a representation that corresponds to a linear system of algebraic equations. Furthermore, we conduct a comparative analysis between the outcomes obtained from the suggested method and those derived from the rationalized Haar wavelet method without employing any shared mathematical methodologies. In order to evaluate the precision and effectiveness of the proposed method, a series of numerical examples have been developed.