In this study, we propose a THz computed tomography (CT) method based on phase contrast, which retrieves the phase shift information at each data point through a phase modulation technique using a Mach-Zehnder interferometer with a continuous wave (CW) source. The THz CT is based on first-generation CT, which acquires a set of projections by translational and rotational scans using a thin beam. From the phase-shift projections, we reconstruct a spatial distribution of refractive indices in a cross section of interest. We constructed a preliminary system using a highly coherent CW THz source with a frequency of 0.54 THz to prove the concept and performed an imaging experiment using phantoms to investigate its imaging features such as artifact-immune imaging, quantitative measurement, and selective detection.