Can anyone help me evaluate the following integral: $$ I(a)=\int_{a}^{\infty} e^{t}\, t^{-a}\, \Gamma(a-1,t)\, dt, $$ where $a\in(0,1)$ is a fixed parameter and $\Gamma(\cdot,\cdot)$ denotes the upper incomplete Gamma function?
One idea is to use the integral representation $$ e^{t} t^{-a}\Gamma(a-1,t)=\dfrac{1}{t}\int_0^{\infty}\dfrac{e^{-ty}}{(1+y)^{2-a}}dy, $$ which gives $$ I(a)=\int_0^\infty\dfrac{E_{1}(ay)}{(1+y)^{2-a}}dy, $$ where $E_1(z)$ denotes the basic exponential integral function. I am not sure how to proceed from here though.
Another idea is to use the integral representation $$ e^{t} t^{-a}\Gamma(a-1,t)=\dfrac{1}{t\Gamma(2-a)}\int_0^{\infty}\dfrac{y^{1-a} e^{-y}}{(y+t)}dy, $$ which gives $$ I(a)=\dfrac{1}{\Gamma(2-a)}\int_0^\infty y^{-a} e^{-y}\log(1+y/a) dy; $$ note that $a\in(0,1)$ so the integral is convergent.