computing all minimal Markov bases of a configuration matrix