This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Attempting to Use circlator running into issues with CANU

Hi!

I am attempting to circularize my assembly and have been running into the following issue, I also posted the slurm script I am trying to run it with.

I am using Canu 1.7.1 if that is of any help as the latest version would be recognized by Circlator.

canu -useGrid=false gnuplotTested=true -assemble genomeSize=0.1m -d /storage/work/dfl5440/Kim_Ritchie_Bacteria/circulization/BC1_INITIAL/03.assemble -p canu -pacbio-corrected /storage/work/dfl5440/Kim_Ritchie_Bacteria/circulization/BC1_INITIAL/02.bam2reads.fasta

The output was:

-- WARNING:
-- WARNING:  Option '-pacbio-corrected' is deprecated.
-- WARNING:  Use options '-corrected -pacbio <files>' in the future.
-- WARNING:

usage:   canu [-version] [-citation] \
              [-haplotype | -correct | -trim | -assemble | -trim-assemble] \
              [-s <assembly-specifications-file>] \
               -p <assembly-prefix> \
               -d <assembly-directory> \
               genomeSize=<number>[g|m|k] \
              [other-options] \
              [-haplotype{NAME} illumina.fastq.gz] \
              [-corrected] \
              [-trimmed] \
              [-pacbio |
               -nanopore |
               -pacbio-hifi] file1 file2 ...

example: canu -d run1 -p godzilla genomeSize=1g -nanopore-raw reads/*.fasta.gz


  To restrict canu to only a specific stage, use:
    -haplotype     - generate haplotype-specific reads
    -correct       - generate corrected reads
    -trim          - generate trimmed reads
    -assemble      - generate an assembly
#!/bin/bash
#SBATCH --job-name=circle_01
#SBATCH --nodes=1
#SBATCH --ntasks=12
#SBATCH --time=20:00:00
#SBATCH --mem=60G

#Variables
export ASSEMBLY=path/to/assembly.fasta
export READS=path/to/barcode01_trimmed100.fastq
export OUTPUT=path/to/directory/BC1_INITIAL

#Command
circlator all --assembler canu --threads $SLURM_NTASKS $ASSEMBLY $READS $OUTPUT
canu circlator

Looks like canu does not like the options that are being fed to it by circlator. Are you able to use a newer version of canu (latest seems to be v.2.2) ? circlator manual does not specify what version of canu it works with. Alternate option is to use SPAdes where compatible versions are defined by circlator.

I attempted to run it with canu v2.2 and 2.1 but circlator could not find the version number when I installed those versions. I would prefer to use canu but at this point I have tried a fair bit of solutions only to continue to run into errors.

0 answers

No answers yet.

Log in to answer this question.