Fix for printing list of inactive reactions for use in pathway gapfilling

The original code used the biochemistry reaction ids instead of the modelreaction ids.