We present results of experimental and theoretical studies of polarization-resolved light transmission through optical fiber with disorder generated in its germanium-doped core via UV radiation transmitted through a diffuser. In samples longer than a certain characteristic length, the power transmitted with preserved polarization is observed to be distributed over all forward-propagating modes, as evidenced by the Rayleigh negative exponential distribution of the near-field intensity at the output surface of the fiber. Furthermore, the transmitted power becomes also equally distributed over both polarizations. To describe the optical properties of the fibers with the experimentally induced disorder, a theoretical model based on coupled-mode theory is developed. The obtained analytical expression for the correlation function describing spatial properties of the disorder shows that it is highly anisotropic. Our calculations demonstrate that this experimentally controllable anisotropy can lead to suppression of the radiative leakage of the propagating modes, so that intermode coupling becomes the dominant scattering process. The obtained theoretical expressions for the polarization-resolved transmission fit very well with the experimental data, and the information extracted from the fit shows that radiative leakage is indeed small. The reported technique provides an easy way to fabricate different configurations of controlled disorder in optical fibers suitable for such applications as random fiber lasers.