In this work, an algorithm for finding numerical solutions of linear fractional delay-integro-differential equations (LFDIDEs) of variable-order (VO) is introduced. The operational matrices are used as discretization technique based on shifted Chebyshev polynomials (SCPs) of the first kind with the spectral collocation method. The proposed VO-LFDIDEs have multiterms of integer, fractional-order derivatives for delayed or nondelayed and mixed Volterra-Fredholm integral terms. The introduced model is a more general form of linear fractional VO pantograph, neutral, and mixed Fredholm–Volterra integro-differential equations with delay arguments. Caputo’s VO fractional derivative operator is used to generate the matrices of the derivative. Operational matrices are presented for all terms. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments. Also, some examples are included to improve the validity and applicability of the techniques. Finally, comparisons between the proposed method and other methods were used to solve this kind of equation.