Genome-wide identification of overexpression and downregulation gene targets based on the sum of covariances of the outgoing reaction fluxes.