add uniprot id instead of replacement

README.md

@@ -33,7 +33,7 @@ This repository contains the current consensus genome-scale metabolic model of _
 
 | Taxonomy | Latest update | Version | Reactions | Metabolites | Genes |
 |:-------|:--------------|:------|:------|:----------|:-----|
-| _Saccharomyces cerevisiae_ | 04-Aug-2023 | 8.6.3 | 4099 | 2768 | 1161 |
+| _Saccharomyces cerevisiae_ | 06-Aug-2023 | 8.6.3 | 4099 | 2768 | 1161 |
 
 # Installation & usage
 

add_uniprot_id/addUniprotID.m

@@ -4,10 +4,23 @@
 fid = fopen('SGD_with_Uniprot.csv');
 uniprot = textscan(fid,'%q %q %q','Delimiter',',','HeaderLines',1);
 for i = 1:length(model.genes)
-    ids.name{1, 1} = char("Uniprot ID");
-    ind = find(strcmp(model.genes{i, 1}, uniprot{1, 1}));
-    ids.value{1, 1} = uniprot{1, 3}{ind, 1};
-    model.geneMiriams{i, 1} = ids;
+    gM = model.geneMiriams(i);
+    gM = gM{1,1};
+    if ~isa(gM, 'struct')
+        ids.name{1, 1} = char('kegg.genes');
+        ids.name{2, 1} = char('ncbigene');
+        ids.name{3, 1} = char('refseq');
+        ids.name{4, 1} = char('uniprot');
+        ids.name{5, 1} = char('ncbiprotein');
+        ind = find(strcmp(model.genes{i, 1}, uniprot{1, 1}));
+        ids.value{4, 1} = uniprot{1, 3}{ind, 1};
+        ids.value{1, 1} = char('NaN');
+        ids.value{2, 1} = char('NaN');
+        ids.value{3, 1} = char('NaN');
+        ids.value{5, 1} = char('NaN');
+        model.geneMiriams{i, 1} = ids;
+        clear ids
+    end
 end
 cd ../code/
 saveYeastModel(model)

model/dependencies.txt

@@ -1,4 +1,4 @@
 MATLAB	9.11.0.1769968 (R2021b)
 libSBML	5.19.0
 RAVEN_toolbox	2.8.3
-COBRA_toolbox	commit fatal: 
+COBRA_toolbox	commit f602a23