The pole-zero plot shows several poles (marked by x) that are stable, but relatively close to the imaginary axis. In your application, you might want to ensure that balred does not discard these near-unstable poles.
To do so, calculate a reduced-order system that preserves the two stable poles within 0.0005 $s^{-1}$ of the imaginary axis.