  Using a dynamic cluster quantum Monte Carlo approximation, we study a
two-dimensional Hubbard model with a small orthorhombic distortion in the
nearest neighbor hopping integrals. We find a large nematic response in the
low-frequency single-particle scattering rate which develops with decreasing
temperature and doping as the pseudogap region is entered. At the same time,
the d-wave superconducting gap function develops an s-wave component and its
amplitude becomes anisotropic. The strength of the pairing correlations,
however, is found to be unaffected by the strong anisotropy, indicating that
d-wave superconductivity can coexist with strong nematicity in the system.
