Gapseq returns models based on SEED reactions even when database is specified as vmh

[anubhavdas0907]

Hello,

I was generating a model using gapseq with database as vmh and pathway database as kegg.

Although gapseq returned a valid model, the reactions in the model have ids of seed and not vmh.

Can you please let me know, how can I get models with reaction ids of vmh and not seed?

code used.

gapseq find -p all -f ./test_out/ -d vmh -l KEGG -v 0 -K 5 -O -g ./genomic.fna

Thanks.