Canu ERROR: File supplied on command line; use -s, -pacbio-raw, -pacbio-corrected, -nanopore-raw, or -nanopore-corrected.
0
1
Entering edit mode
4.8 years ago
lbrutscher ▴ 10

Hi everyone,

I am trying to assemble a bacterial genome using PacBio reads and short illumina reads from a metagenomic pool of both eurkaryotic host and bacterial sequences.

In an ideal world, I would have been able to culture the bacteria and obtain pure DNA to sequence and assemble from, but its cultivation has proved challenging. So, instead we have a ton of metagenomic data sequenced obtained from whole host DNA. We have roughly ~70 million Illumina MiSeq paired end reads (2x150) and ~120,000 PacBio reads (5-20 KB). I suspect that I have roughly 3x coverage of the bacterial genome with the PacBio reads and ~20-30x coverage with the illumina reads (with 80-90% of the reads likely being host). This project is further complicated by the fact that the host genome is not fully sequenced, so I cannot remove known host sequence to simplify analysis and de novo assembly.

I have been trying to use Canu 1.4 to correct my PacBio reads using fastq as input and the following code:

./canu -correct -pacbio-raw m170212_104155_42146_c101160022550000001823267305221711_s1_p0.3.subreads.fastq -d subreads3-auto -p subreads3 genomeSize=1.5m gnuplotTested=true  corMinCoverage=0 errorRate=0.035 corMaxEvidenceErate=0.15


But I keep getting this error: ERROR: File supplied on command line; use -s, -pacbio-raw, -pacbio-corrected, -nanopore-raw, or -nanopore-corrected.

It appears that Canu is not recognizing my input fastq file "m170212_104155_42146_c101160022550000001823267305221711_s1_p0.3.subreads.fastq". Has anyone else received this error and know what to do about it?

Assembly canu pacbio • 2.3k views
0
Entering edit mode

Have you tried moving -pacbio-raw m170212_104155_42146_c101160022550000001823267305221711_s1_p0.3.subreads.fastq to the end of the command-line? Some programs are finicky about the order of input options. Also all those - are hyphens (i.e. not smart hyphens etc).

0
Entering edit mode

Thanks for the response. Moving -pacbio-raw m170212_104155_42146_c101160022550000001823267305221711_s1_p0.3.subreads.fastq to the end does not help.

0
Entering edit mode

Check your java version. Use full path to canu /Full/Path/to/bin/Canu

It will be better to show all your error log

0
Entering edit mode

Hi, I am experiencing exactly the same problem, when I try to run the tutorial on ecoli supplied on the CANU readthedocs http://canu.readthedocs.io/en/latest/quick-start.html. I download the supplied pacbio file and run the script, and then: ERROR: File supplied on command line; use -s, -pacbio-raw, -pacbio-corrected, -nanopore-raw, or -nanopore-corrected. Have you cracked this problem?

0
Entering edit mode

Any particular reason for using 1.4 instead of the latest release (1.6)?

edit. just realized OP posted quite a long time ago

0
Entering edit mode

When I'm experiencing the same problem now, it is with v1.6.