The nucleon(N)-Omega(Ω) system in the S-wave and spin-2 channel (S25) is studied from the (2+1)-flavor lattice QCD with nearly physical quark masses (mπ≃146MeV and mK≃525MeV). The time-dependent HAL QCD method is employed to convert the lattice QCD data of the two-baryon correlation function to the baryon-baryon potential and eventually to the scattering observables. The NΩ(S25) potential, obtained under the assumption that its couplings to the D-wave octet-baryon pairs are small, is found to be attractive in all distances and to produce a quasi-bound state near unitarity: In this channel, the scattering length, the effective range and the binding energy from QCD alone read a0=5.30(0.44)(−0.01+0.16)fm, reff=1.26(0.01)(−0.01+0.02)fm, B=1.54(0.30)(−0.10+0.04)MeV, respectively. Including the extra Coulomb attraction, the binding energy of pΩ−(S25) becomes BpΩ−=2.46(0.34)(−0.11+0.04)MeV. Such a spin-2 pΩ− state could be searched through two-particle correlations in p-p, p-nucleus and nucleus-nucleus collisions.