AnantharamanLab / METABOLIC

A scalable high-throughput metabolic and biogeochemical functional trait profiler
173 stars 43 forks source link

Hi, when running METABOLIC-C (either test or my data), it fails to create the All_gene_collections.gene file, and all the downstream mapping fails. After "Drawing element cycling diagrams..." it states "No output file specified!". Please help - many thanks in advance! #65

Open MaxRubinBlum opened 2 years ago

MaxRubinBlum commented 2 years ago

Describe the bug A clear and concise description of what the bug is.

To Reproduce Steps to reproduce the behavior:

  1. Go to '...'
  2. Click on '....'
  3. Scroll down to '....'
  4. See error

Expected behavior A clear and concise description of what you expected to happen.

Screenshots If applicable, add screenshots to help explain your problem.

Desktop (please complete the following information):

Smartphone (please complete the following information):

Additional context Add any other context about the problem here.

ChaoLab commented 2 years ago

Hi, It can be due to some dependencies are not well-installed. Were the gene files generated in the genome folder? The "All_gene_collections.gene" is the concatenated file from gene files in the genome folder.

MaxRubinBlum commented 2 years ago

Yes, the gene files were generated.

MaxRubinBlum commented 2 years ago

Hi, thanks for the response! Still cannot make this work. I tried to override this by concatenating the files manually and re-running the script, but the file apparently gets removed, and the issue still remains. Bowtie 2 thus fails to build an index. I am running a manual installation, as my conda is quite broken and fails to build environments due to conflict issues. I think all the dependencies were installed. What am I missing? Thanks for the help.

ChaoLab commented 2 years ago

It will be much easier to run within conda since there are many dependencies

MaxRubinBlum commented 2 years ago

I understand, but this is not feasible at the moment. Which dependency could be missing to cause this issue?

MaxRubinBlum commented 2 years ago

Dear Zhichao,

I am still having a problem resolving this issue, your help is highly appreciated.

Many thanks, Maxim

From: Zhichao Zhou @.> Sent: Saturday, 12 March 2022 21:37 To: AnantharamanLab/METABOLIC @.> Cc: Maxim Rubin Blum @.>; Author @.> Subject: Re: [AnantharamanLab/METABOLIC] Hi, when running METABOLIC-C (either test or my data), it fails to create the All_gene_collections.gene file, and all the downstream mapping fails. After "Drawing element cycling diagrams..." it states "No output file ...

It will be much easier to run within conda since there are many dependencies

- Reply to this email directly, view it on GitHubhttps://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2FAnantharamanLab%2FMETABOLIC%2Fissues%2F65%23issuecomment-1065949814&data=04%7C01%7Cmrubin%40ocean.org.il%7C49e746f52cbf49c6b51f08da045f9c93%7Cb28ad7a9bfc042c2ae1506fa7cf73a5c%7C0%7C0%7C637827105949376363%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=Eghu%2Fw2Y4pBj6F7gkbRNWJxUu7H5uerM2Bs9vMMZzM8%3D&reserved=0, or unsubscribehttps://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fnotifications%2Funsubscribe-auth%2FAQGZGJ2VQOBL6XFILPWXX7DU7TW35ANCNFSM5QQILXCQ&data=04%7C01%7Cmrubin%40ocean.org.il%7C49e746f52cbf49c6b51f08da045f9c93%7Cb28ad7a9bfc042c2ae1506fa7cf73a5c%7C0%7C0%7C637827105949376363%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=L0fagIIrXqQ6P9oWDJxhR7OYUG%2Bfn2pgnIQM1APSTiw%3D&reserved=0. Triage notifications on the go with GitHub Mobile for iOShttps://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fapps.apple.com%2Fapp%2Fapple-store%2Fid1477376905%3Fct%3Dnotification-email%26mt%3D8%26pt%3D524675&data=04%7C01%7Cmrubin%40ocean.org.il%7C49e746f52cbf49c6b51f08da045f9c93%7Cb28ad7a9bfc042c2ae1506fa7cf73a5c%7C0%7C0%7C637827105949376363%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=RodXn63njhxR9DJafwB4LH732Uq%2BR1UxvKNj%2BwXXjic%3D&reserved=0 or Androidhttps://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fplay.google.com%2Fstore%2Fapps%2Fdetails%3Fid%3Dcom.github.android%26referrer%3Dutm_campaign%253Dnotification-email%2526utm_medium%253Demail%2526utm_source%253Dgithub&data=04%7C01%7Cmrubin%40ocean.org.il%7C49e746f52cbf49c6b51f08da045f9c93%7Cb28ad7a9bfc042c2ae1506fa7cf73a5c%7C0%7C0%7C637827105949376363%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=RQUbqoh%2BeI0uGuk6LUbAy2d1nVhXA1Y%2B7s6YwBSVqA4%3D&reserved=0. You are receiving this because you authored the thread.Message ID: @.**@.>>

ChaoLab commented 2 years ago

Can you try to use Bowtie2 to make the index independent of the METABOLIC? Maybe you can also try to run the following Bowtie2 mapping steps for troubleshooting. Sometimes, it could be just a dependency issue

H1889 commented 2 years ago

Hi, I have had the same problem and I think it is an error of the perl script "METABOLIC-C.pl". The line 1904 is:

system ("bowtie2-build $output/All_gene_collections.gene $output/All_gene_collections.gene.scaffold --quiet");

but must be written:

system ("bowtie2-build --quiet $output/All_gene_collections.gene $output/All_gene_collections.gene.scaffold");

The bowtie2 manual says: Usage: bowtie2-build [options]*

so the options must be written BEFORE the reference and base name not AFTER.

Best regards

ChaoLab commented 2 years ago

@H1889 Hi, I changed the script accordingly

patriciatran commented 2 years ago

Hi @MaxRubinBlum , have you been able to test with our updates METABOLI-c.pl scripts (see comments above) ? Thanks for letting us know.

MaxRubinBlum commented 2 years ago

Hi @patriciatran, not yet, will certainly give it a go. Many thanks!

H1889 commented 2 years ago

I tested it and worked, however this error does not: ...: [2022-04-10 01:00:51] Drawing element cycling diagrams finished [2022-04-10 01:00:51] Drawing metabolic handoff diagrams... mv: no se puede efectuar stat' sobre 'COMMUNITY_BACTERIA_OUT/newdir/Bar_plot/bar_plot_input_1.pdf'mv: no se puede efectuarstat' sobre 'COMMUNITY_BACTERIA_OUT/newdir/Bar_plot/bar_plot_input_1.pdf': No existe el archivo o el directorio : No existe el archivo o el directorio mv: no se puede efectuar stat' sobre 'COMMUNITY_BACTERIA_OUT/newdir/Bar_plot/bar_plot_input_2.pdf'mv: no se puede efectuarstat' sobre 'COMMUNITY_BACTERIA_OUT/newdir/Bar_plot/bar_plot_input_2.pdf': No existe el archivo o el directorio : No existe el archivo o el directorio ....