We present the numerical code TAURUS_pav that computes projected expectation values of various operators between two real general Bogoliubov quasi-particle states represented in a Spherical Harmonic Oscillator (SHO) basis. The present version of the code allows one to project onto a good: number of protons Z, number of neutrons N, number of nucleons A, total angular momentum J and its third component M_J, parity P. The model space considered is invariant under spatial and isospin rotations butno specific SHO basis is assumed such that the code can carry out either valence-space or no-core calculations (depending on the hamiltonian file used as input). In addition, no number parity is assumed for the Bogoliubov quasi-particle states such that the code can be used to handle even-even, odd-even and odd-odd nuclei.