Stress on seismogenic faults provides critical information about how much elastic energy is stored in the crust and released by earthquakes, which is crucial in understanding earthquake energetics and recurrence. However, determining post-earthquake stress states on faults remains challenging because current borehole methods are rarely applicable to damaged fault zone rocks. We applied neutron texture analysis to gouge samples of the 1999 Chi-Chi earthquake in Taiwan to infer the stress state after the earthquake. Results indicate that the clay fabric within the principal slip zone is orthogonal to the fault plane, whereas outside the principal slip zone the fabric is predominantly parallel to the bedding-parallel fault plane. We suggest that the clay fabric in the slip zone was first neutralized by the coseismic fluidization caused by thermal pressurization and later re-oriented to the new direction of post-earthquake principal stress. Such stress orientation is consistent with the orientations inferred from core-scale fault slip data and dislocation models constrained from global navigation satellite system displacements. If thermal pressurization is a ubiquitous process during earthquakes, gouge fabrics can be used to help probe the post-earthquake stress state of faults.