nanopolish icon indicating copy to clipboard operation
nanopolish copied to clipboard

Can you restart nano polish?

Open yamkela-mg opened this issue 2 years ago • 1 comments

Hi there,

I am running nano polish with this command:

python3 nanopolish_makerange.py draft.fa | parallel --results nanopolish.results -P 20
nanopolish variants --consensus -o polished.{1}.vcf -w {1} -r reads.fa -b reads.sorted.bam -g draft.fa -t 56 --min-candidate-frequency 0.1

  1. I am polishing a plant genome (1.2 GB) and it's taking more than a week, is this normal? Is there a way to speed it?

  2. If the polishing job gets killed because it has exceeded the set wall time, is there a way of restarting? If I re-run it in the same directory as the previous run, will it continue or start from the beginning?

yamkela-mg avatar May 12 '22 15:05 yamkela-mg

HI @yamkela-mg,

I don't recommend nanopolish for polishing genome assemblies anymore, medaka gives consistently better results. If you wish to continue your nanopolish run I would suggest checking which .vcf files are complete (perhaps using the job's exit status), then remove any incomplete .vcf files and wrapping the commands in a shell script that only runs nanopolish if the target .vcf file does not exist.

I hope this helps. Jared

jts avatar May 12 '22 20:05 jts