Uncertainty always lives with us. We cannot take exact measurement of initial conditions or parameters values in a mathematical model. As humans, we are remaining alive in an environment where the uncertainties lie in the modelling of physical phenomena. There might be some incomplete information or estimation of the parameter or initial values. To handle uncertainty, we use fuzzy operators rather than classical operators. In this paper, we study a model of HIV-1 infection by taking uncertainty in the initial data under Caputo fractional operator. We explore the existence and uniqueness of the results through fixed-point theory. We study the Ulam–Hyres stability of the considered model. By using the fuzzy Laplace Adomian decomposition method, numerical results are obtained for specific fuzzy initial conditions. To better understand the behaviour of the fuzzy solution, we present the obtained numerical results graphically for various fractional orders where the uncertainty lies in [0, 1].