carveme icon indicating copy to clipboard operation
carveme copied to clipboard

GPR information is missing in output file

Open Shakra-A opened this issue 3 years ago • 6 comments

Hi. I am new to all this work. I used CarveMe for making a metabolic model from protein fasta file. I used the following code. $ carve genome.faa --output model.xml, from the documentation. The fasta file was consist of locus tags and their protein sequences. The result I got should have locus tags in the column GPR but the column is empty. (I converted the sbml file to excel file for data visualization). Kindly help me to solve this problem.

Shakra-A avatar Jun 09 '21 04:06 Shakra-A

Sorry for the late reply. Can you send me the SBML file?

cdanielmachado avatar Jun 30 '21 12:06 cdanielmachado

yes I can but here in reply the file format is not supported. Kindly send me your email address so I could send it there.

Shakra-A avatar Jul 01 '21 01:07 Shakra-A

cdanielmachado at gmail

cdanielmachado avatar Jul 01 '21 07:07 cdanielmachado

If you run with -v you will see the error message: The input genome did not match sufficient genes/reactions in the database.

There is something strange with your fasta file, most of the proteins have a really short AA sequence. Was this fasta file built from a fully assembled genome? Also, do you know which organism is this?

cdanielmachado avatar Jul 08 '21 07:07 cdanielmachado

This organism does not have a GBK file. I used automatic pipe lines for annotation of the whole genome. It is a bacteria. The fasta file is from the automatic annotation proteine seq.

On Thu, Jul 8, 2021 at 4:22 PM Daniel Machado @.***> wrote:

If you run with -v you will see the error message: The input genome did not match sufficient genes/reactions in the database.

There is something strange with your fasta file, most of the proteins have a really short AA sequence. Was this fasta file built from a fully assembled genome? Also, do you know which organism is this?

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/cdanielmachado/carveme/issues/131#issuecomment-876197961, or unsubscribe https://github.com/notifications/unsubscribe-auth/ASUSI6JXMPZMVQU4VZLYUSLTWVG4HANCNFSM46LIBWSA .

-- Regards, ....SHAKRA....

Shakra-A avatar Jul 08 '21 08:07 Shakra-A

Maybe you can run it in eggnog-mapper and see how many metabolic enzymes it predicts: http://eggnog-mapper.embl.de/

Or try to run a reconstruction with modelSEED: https://modelseed.org/

cdanielmachado avatar Jul 08 '21 08:07 cdanielmachado