The bulk photovoltaic effect is an experimentally verified phenomenon by which a direct charge current is induced within a non-centrosymmetric material by light illumination. Calculations of its intrinsic contribution, the shift current, are nowadays amenable from first-principles employing plane-wave bases. In this work, we present a general method for evaluating the shift conductivity in the framework of localized Gaussian basis sets that can be employed in both the length and velocity gauges, carrying the idiosyncrasies of the quantum-chemistry approach. The (possibly magnetic) symmetry of the system is exploited in order to fold the reciprocal space summations to the representation domain, allowing us to reduce computation time and unveiling the complete symmetry properties of the conductivity tensor under general light polarization.