Skew Schur function

edit

I have worked with the Schur function using the following command for example

sage: Sym = SymmetricFunctions(QQ)
sage: P = Sym.p()  # to expand symmetric functions in power symmetric basis
sage: s = Sym.schur()
sage: P(s([2,2]))
1/12*p[1, 1, 1, 1] + 1/4*p[2, 2] + (-1/3)*p[3, 1]

Recently I need to work with the skew Schur function. I want to expand the shifted schur function $$ s_{\lambda}(x_1 +y , x_2 +y ,\ldots , x_n +y) $$ I am guessing that it can be the following identity $$ s_{\lambda}(x_1 +y , x_2 +y ,\ldots , x_n +y) = \sum_{\mu} s_{\lambda \setminus mu}(y) s_{\mu}(x_1 , x_2, \ldots, x_n) $$ For $y=1$ my idenity can be proven as I am sure there is an expression. The expression $s_{\lambda \setminus mu}(y)$ is skew Schur function. I am wondering if there are packages in sagemath where I can work with skew schur function and verify my identity is correct or not.