@@ -117,11 +117,13 @@ function Run(){
117117 function process (){
118118 $ header = $ this ->GetReadFile ()->Read (200000 );
119119 $ header_arr = explode ("\t" , $ header );
120-
121- if (count ($ header_arr ) != 40 )
120+ $ n = 41 ;
121+ $ c = count ($ header_arr );
122+ if ($ c != $ n )
122123 {
123124 echo PHP_EOL ;
124- trigger_error ("Header format is different than expected, please update the script " ,E_USER_ERROR );
125+ print_r ($ header_arr );
126+ trigger_error ("Expected $ n columns, found $ c . please update the script " ,E_USER_ERROR );
125127 exit ;
126128 }
127129
@@ -164,9 +166,10 @@ function process(){
164166 $ refseq_mappeddatasuppliedbyNCBI = $ fields [34 ];
165167 $ uniprot_id_mappeddatasuppliedbyUniProt = $ fields [35 ];
166168 $ ensembl_id_mappeddatasuppliedbyEnsembl = $ fields [36 ];
167- $ ucsc_id_mappeddatasuppliedbyUCSC = $ fields [37 ];
168- $ mouse_genome_database_id_mappeddatasuppliedbyMGI = $ fields [38 ];
169- $ rat_genome_database_id_mappeddatasuppliedbyRGD = $ fields [39 ];
169+ $ vega_id_mappeddatasuppliedbyVega = $ fields [37 ];
170+ $ ucsc_id_mappeddatasuppliedbyUCSC = $ fields [38 ];
171+ $ mouse_genome_database_id_mappeddatasuppliedbyMGI = $ fields [39 ];
172+ $ rat_genome_database_id_mappeddatasuppliedbyRGD = $ fields [40 ];
170173
171174 $ id_res = $ id ;
172175 $ id_label = "Gene Symbol for " .$ approved_symbol ;
@@ -464,6 +467,15 @@ function process(){
464467 }
465468 }
466469
470+ if (!empty ($ ucsc_id_mappeddatasuppliedbyVega )){
471+ $ ucsc_id_mappeddatasuppliedbyVega = explode (", " , $ ucsc_id_mappeddatasuppliedbyVega );
472+ foreach ($ ucsc_id_mappeddatasuppliedbyVega as $ vega_id ) {
473+ parent ::AddRDF (
474+ parent ::triplify ($ id_res , $ this ->getVoc ()."x-vega " , "vega: " .trim ($ vega_id )).
475+ parent ::describeProperty ($ this ->getVoc ()."x-vega " , "Vega entry " )
476+ );
477+ }
478+ }
467479 if (!empty ($ ucsc_id_mappeddatasuppliedbyUCSC )){
468480 $ ucsc_id_mappeddatasuppliedbyUCSC = explode (", " , $ ucsc_id_mappeddatasuppliedbyUCSC );
469481 foreach ($ ucsc_id_mappeddatasuppliedbyUCSC as $ ucsc_id ) {
0 commit comments