{"id":2508,"date":"2019-03-01T10:59:02","date_gmt":"2019-03-01T10:59:02","guid":{"rendered":"http:\/\/ri.itservices.manchester.ac.uk\/csf3\/?page_id=2508"},"modified":"2021-04-07T11:55:12","modified_gmt":"2021-04-07T10:55:12","slug":"canu","status":"publish","type":"page","link":"https:\/\/ri.itservices.manchester.ac.uk\/csf3\/software\/applications\/canu\/","title":{"rendered":"Canu"},"content":{"rendered":"

Overview<\/h2>\n

Canu<\/a> is a fork of the Celera Assembler, designed for high-noise single-molecule sequencing (such as the PacBio RS II\/Sequel or Oxford Nanopore MinION).<\/p>\n

The application is written in C++ and can be used in batch or interactively (via the qrsh command).<\/p>\n

Version 2.0 and 1.8 are installed on the CSF.<\/li>\n

Restrictions on use<\/h2>\n

There are no restrictions on using this software on the CSF. The application is released under the GPLv2 license.<\/p>\n

Set up procedure<\/h2>\n

To access the software you must load ONE of the modulefiles:<\/p>\n

\r\nmodule load apps\/binapps\/canu\/2.1.1\r\nmodule load apps\/binapps\/canu\/2.0\r\nmodule load apps\/canu\/1.8\r\n<\/pre>\n
We recommend that you load the modulefile in your jobscript – see examples below.<\/p>\n
Running the application<\/h2>\n
Please do not run canu on the login node. Jobs should be submitted to the compute nodes via batch or run interactively via the qrsh command. Note that canu has the ability to submit batch jobs on your behalf. However this does not work correctly on the CSF (because canu will try to submit jobs from other jobs which is not supported on the CSF). So you must run it from a jobscript and submit the job to the batch system in the usual way.<\/p>\n
Serial batch job submission<\/h3>\n
Example jobscript<\/p>\n
\r\n#!\/bin\/bash --login\r\n#$ -cwd             # Job will run from the current directory\r\n\r\n# Load the software\r\n\r\nmodule load apps\/binapps\/canu\/2.0\r\n\r\n# $NSLOTS is automatically set to 1 in a serial job\r\n# You MUST set the maxThreads, mhapThreads and maxMemory options\r\n# to tell canu how much resource it can use. We assume 4GB per core.\r\n# The useGrid=false option is needed to run from a batch job.\r\n\r\ncanu -d run1 -p godzilla genomeSize=1g maxThreads=$NSLOTS mhapThreads=$NSLOTS maxMemory=$((4*NSLOTS)) \\\r\n    useGrid=false -nanopore-raw reads\/*.fasta.gz\r\n<\/pre>\n
Submit the jobscript using:<\/p>\n
qsub scriptname<\/em><\/pre>\n
where scriptname<\/em> is the name of your jobscript.<\/p>\n
Parallel batch job submission<\/h3>\n
Example jobscript<\/p>\n
#!\/bin\/bash --login\r\n#$ -cwd             # Job will run from the current directory\r\n#$ -pe smp.pe 8     # Number of cores (2-32 permitted)\r\n\r\n# Load the software\r\n\r\nmodule load apps\/binapps\/canu\/2.0\r\n\r\n# $NSLOTS is automatically set to the number given on the -pe line above.\r\n# You MUST set the maxThreads, mhapThreads and maxMemory options\r\n# to tell canu how much resource it can use. We assume 4GB per core.\r\n# The useGrid=false option is needed to run from a batch job.\r\n\r\ncanu -d run1 -p godzilla genomeSize=1g maxThreads=$NSLOTS mhapThreads=$NSLOTS maxMemory=$((4*NSLOTS)) \\        \r\n   useGrid=false -nanopore-raw reads\/*.fasta.gz            #                    #\r\n                                                           #                    # Ensure 4GB\/core \r\n                                                           #                    # is used    \r\n                                                           #\r\n                                                           # If your pipeline uses the MHAP algorithm\r\n                                                           # this ensures Java will use the correct \r\n                                                           # number of cores\r\n<\/pre>\n
Submit the jobscript using:<\/p>\n
qsub scriptname<\/em><\/pre>\n
where scriptname<\/em> is the name of your jobscript.<\/p>\n
Further info<\/h2>\n