We calculate the QCD cross-over temperature, the equation of state and fluctuations of conserved charges at finite density by analytical continuation from imaginary to real chemical potentials. Our calculations are based on new continuum extrapolated lattice simulations using the 4stout staggered actions with a lattice resolution up to $N_t=16$. The simulation parameters are tuned such that the strangeness neutrality is maintained, as it is in heavy ion collisions.