Reproducible example for https://github.com/opencobra/cobrapy/issues/245

bug245.py

#/usr/bin/env python

from cobra.test import create_test_model
from cobra.flux_analysis import single_gene_deletion

cobra_model = create_test_model("textbook")
dels = {"b0008": 0.87, "b0114": 0.71, "b0116": 0.56, "b2276": 0.11, "b1779": 0.00}

rates, statuses = single_gene_deletion(cobra_model, gene_list=dels.keys(),
    method="moma", solver="mosek")
for gene in dels:
    if statuses[gene] != 'optimal':
        raise ValueError("Non-optimal status: {}".format(statuses[gene]))