This study proposes analytical and computational methods for the solution of the sliding frictional contact problem of an anisotropic laterally graded layer loaded by an arbitrarily shaped rigid stamp. The plane-strain orthotropy prevails in the layer which is bonded to a rigid foundation. Each of four orthotropic stiffness coefficients is exponentially varied through the lateral direction of the elastic layer. The Fourier transformations of the field variables are employed in the formulation. The gradient of a displacement component on the surface is then converted to a singular integral equation of the second kind. The singular integral equation is solved by means of the Gauss-Jacobi quadrature integration techniques, a collocation method, and a recursive integration method for the Cauchy integral considering the flat and triangular stamp profiles. The finite element method solutions of the same contact problems are performed using the augmented Lagrange method which is implemented in virtue of ANSYS design parametric language. An iterative algorithm is additionally utilized for the (incomplete) triangular stamp problem to conveniently reach the solutions for predetermined contact lengths. The convergence and comparative analyses are carried out to elucidate the trustworthiness of the analytical and computational methods proposed. Moreover, the parametric analyses infer that the contact-induced damage risks can be effectively alleviated upon tuning the degree of orthotropy and the lateral heterogeneity of the elastic layer.