% for the pathway info and how it contradict GE data -- these can be post-processing



topCla(rs-{ReactionID,LimitingType,State,Time},[reaction_state0(ReactionID,LimitingType,State,Time)]).

topCla(met-{SID,Change,Time},[metConcentration(SID,Change,Time)]). % use the same predicate as concentration(MID,Change,Time)? Will it be used for recursive call? no,