Entering edit mode
10 months ago
dlera.lozano
▴
10
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
Looks like
canu
does not like the options that are being fed to it bycirclator
. Are you able to use a newer version ofcanu
(latest seems to be v.2.2) ?circlator
manual does not specify what version ofcanu
it works with. Alternate option is to useSPAdes
where compatible versions are defined bycirclator
.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.