Discrepancy when calculating LTI system output using inverse z-Transform
As Matt L pointed out, your initial condition is wrong. y[0] is NOT 0, it's 4. This is obvious from the difference equation. Once you use y[0] = 4, both results come out to be the same.