gbouras13/hybracter

[BUG] No plasmid in final output

Closed this issue · 4 comments

Hi,

I have a problem with the assembly since yesterday. I assemble bacterial genomes which contain plasmids with Hybracter and since yesterday, I haven't any plasmids in the final output. It worked very well last week because I assembled this same genome with hybracter and I had the plasmid in final output.
Do you know why it does that ? And do you know how to fix it ?
Thanks

Corentin

Hi @CorentinEscobar ,

What versions of hybracter and plassembler are you using (before and after)?

George

Hi George,

Looking in the intermediate log files I noticed that the problem came from Plassembler, and more precisely from SPADES when unicycler is running. SPADES goes into error (OS return -11) and the plasmids are therefore not present in the output file. So in the final output of hybracter doesn't find the plasmids.
According to the issues on the SPADES github, this could be due to lack of memory. By reducing the number of threads I no longer have this problem. However, it takes a little longer ;)

You may have another way to resolve this problem while keeping a high number of threads for all steps except unicycler ?

Thank you !

Corentin

Hi @CorentinEscobar ,

Glad to hear it worked in the end - the best way is to modify this line

threads: config.resources.big.cpu

to e.g. (for 1 thread)

threads: 1

George

Ok perfect ! Thank you @gbouras13
Corentin