Yes, it is the Bohr radius (which I'll call 'a'). The ground state wavefunction is:
ψ1s(r)=(1/π1/2a3/2)e-r/a.
The probability density is |ψ1s(r)|2, which is:
ρ1s(r)=|ψ1s(r)|2=(1/πa3)e-2r/a.
But that function is not going to give you the most probable radius. You have to take into account the fact that ρ1s is in spherical coordinates, whose volume element is:
dV=r2sin(φ)dr dθ dφ.
So, when you integrate ρ1s over all space, it gets multiplied by r2. Furthermore, since ψ1s is spherically symmetric, you can integrate over θ and φ to get what is called the radial probability density P1s(r):
P1s(r)=(4/a3)r2e-2r/a.
If you optimize this function, you will find that it has a relative maximum at r=a, the Bohr radius.