Let n ⩾ 2 and s ∈ (n − 2, n). Assume that Ω ⊂ ℝn is a one-sided bounded non-tangentially accessible domain with s-Ahlfors regular boundary and σ is the surface measure on the boundary of Ω, denoted by ∂Ω. Let β be a non-negative measurable function on ∂Ω satisfying \(\beta \in {{L}^{{q_0}}}(\partial\Omega, \, \sigma)\) with \({q_0} \in ({s \over {s+2-n}},\, \infty]\) and β ⩾ a0 on E0 ⊂ ∂Ω, where a0 is a given positive constant and E0 ⊂ ∂Ω is a σ-measurable set with σ(E0) > 0. In this paper, for any f ∈ Lp(∂Ω, σ) with p ∈ (s/(s + 2 − n), ∞], we obtain the existence and uniqueness, the global Hölder regularity, and the boundary Harnack inequality of the weak solution to the Robin problem \(\begin{cases}{-{\rm div}(A \nabla u) = 0} & {\rm in} \; \Omega,\\A \nabla u \cdot \nu + \beta u = f & {\rm on} \; \partial\Omega,\end{cases}\) where the coefficient matrix A is real-valued, bounded, and measurable, which satisfies the uniform ellipticity condition, and where ν denotes the outward unit normal to ∂Ω. Furthermore, we establish the existence, upper bound pointwise estimates, and the Hölder regularity of Green’s functions associated with this Robin problem. As applications, we further prove that the harmonic measure associated with this Robin problem is mutually absolutely continuous with respect to the surface measure σ and also provide a quantitative characterization of the mutual absolute continuity at small scales. These results extend the corresponding results established by David et al. (2024) via weakening their assumption that β is a given positive constant.