I used Mercury, Venus, Earth, Moon, Jupiter barycenter, Saturn barycenter, Uranus barycenter and Pluto barycenter. I'm using the same value of mu (G * sun mass) as in the Apophis simulation where our results were identical.
Can you try the first 4000 years at a slower time step? That probably won't take too long using Picard.
Also, very important in this simulation are the masses of Earth and Mars. I used (kilograms):
5.97369125232006E+24 Earth
6.41853544762209E+23 Mars