We report measurements of the NMR lineshape of 4.685% abundant 29 Si in single-crystal silicon for four different crystallographic orientations relative to the applied magnetic field. To avoid significant inhomogeneous broadening, the sample crystals are immersed in a susceptibility matched liquid and the proton NMR lineshape in the liquid is used to measure the residual susceptibility mismatch and to shim the applied field. The observed lineshapes are in good agreement with disorder-averaged spin dynamics simulations performed on a 4 × 4 × 4 unit cell lattice with periodic boundary conditions. The splitting of resolved doublet features is used to measure the asymmetry ∆J of the nearest-neighbor indirect dipole or J-coupling tensor, with the result ∆J = 90 ± 20 Hz. All resolved doublets can be identified with specific nearest and next-nearest neighbor isolated spin pairs that occur in the disorder ensemble.