Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Canu skips where it ended when restart #2104

Closed
Yutang-ETH opened this issue Mar 18, 2022 · 4 comments
Closed

Canu skips where it ended when restart #2104

Yutang-ETH opened this issue Mar 18, 2022 · 4 comments

Comments

@Yutang-ETH
Copy link

Yutang-ETH commented Mar 18, 2022

Hi,

I am using Canu 2.2 and here's my command:

canu \
    -trim-assemble \
    -p mygenome -d mygenome_canu \
    useGrid=false \
    -maxMemory=490g \
    -maxThreads=48 \
    'maxInputCoverage=100' \
    'batOptions=-eg 0.12 -sb 0.01 -dg 3 -db 3 -dr 1 -ca 500 -cp 50' \
    'correctedErrorRate=0.12' \
    genomeSize=2.8g \
    -pacbio-hifi myreads.fa 

Due to an outage of our machine, my Canu program was terminated. When I restarted it, I found it just skips the trimming stage to next unitigging stage and doesn't resume the unfinished overlapper step in trimming. How could I fix this?

Thank you very much.
Best wishes,
Yutang

@skoren
Copy link
Member

skoren commented Mar 18, 2022

Are you sure trimming wasn't already completed? What was the latest in the logs before it terminated? What is the contents of my genome.seqStore folder?

@Yutang-ETH
Copy link
Author

Thank you very much for your quick reply.
Yes, trimming wasn't completed because the total overlaps should be 1345, but only 900 overlaps were completed.
The last overlap file I saw is 1050.out, however there are files missing between 900.out and 1050.out. I guess this is associated with our storage outage, which somehow made canu skip some jobs.
I also checked seqStore folder, files or the folder called version.001 in there are all created by the first day I started Canu, no recent files.
Is it possible to tell Canu to start from overlap 901?

@skoren
Copy link
Member

skoren commented Mar 18, 2022

What's the full log from the restart say? It's possible you're hitting #1781, try using -untrimmed instead of -trim-assemble.

@Yutang-ETH
Copy link
Author

Thank you very much for your suggestion, I will try.

Best wishes,
Yutang

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants