gene deletion computed in separate function

The logic used to determine which reactions are inactivated by a gene
deletion has been moved from delete_model_genes into a standalone
function find_gene_knockout_reactions