In this work, the Boubaker wavelet collocation method is used to numerically investigate a fractional-order fowlpox disease model. The intricate dynamics of transmission between birds and mosquitoes are captured by the analysis of the fowlpox model, which was developed using fractional-order derivatives. The model's solutions are approximated using the Boubaker wavelet technique, which is well-known for its effectiveness and precision in resolving fractional ordinary differential equations (FODEs). Using the Boubaker wavelets, we built the integrated operational matrices. To effectively solve fractional-order systems, the Boubaker wavelet collocation method (BWCM) is employed. Numerical results produced using the Runge-Kutta and NDSolve methods are compared to show how successful the proposed method is at solving fractional-order disease models. In contrast to traditional numerical methods, the BWCM offers better accuracy and efficiency while consuming minimal computer power. To determine how fractional-order parameters and other important factors affect the disease's dynamics, sensitivity studies have been carried out. The study reveals new insight on the dynamics of fowlpox transmission and emphasizes the value of fractional-order models and sophisticated numerical methods in epidemiological studies.