Question: denovo pacbio assembly
0
gravatar for 1234anjalianjali1234
9 weeks ago by
India
1234anjalianjali12340 wrote:

Hi everyone,

I have tried to run the tutorial on supplied on the CANU readthedocs http://canu.readthedocs.io/en/latest/quick-start.html. I downloaded the supplied pacbio file and run the script, but it is giving the error:

ERROR: *Don't panic, but a mostly harmless error occurred and Canu stopped.
Canu release v1.5 failed with:
ERROR:  no reads supplied, and can't find any library information in gkpStore*

my command:

kunal@kunal-CELSIUS-M740Bpower:~/canu$ canu-1.5/Linux-amd64/bin/canu -p ecoli -d ecoli-auto genomeSize=4.8 pacbio.fastq

please help.

assembly • 230 views
ADD COMMENTlink modified 9 weeks ago by colindaven410 • written 9 weeks ago by 1234anjalianjali12340

What's the output when you run:

pwd
readlink -f pacbio.fastq
head -n 4 pacbio.fastq
ADD REPLYlink modified 9 weeks ago • written 9 weeks ago by Ram12k

kunal@kunal-CELSIUS-M740Bpower:~/canu$ pwd

/home/kunal/canu

kunal@kunal-CELSIUS-M740Bpower:~/canu$ readlink -f pacbio.fastq

/home/kunal/canu/pacbio.fastq

kunal@kunal-CELSIUS-M740Bpower:~/canu$ head -n 4 pacbio.fastq

@m141013_011508_sherri_c100709962550000001823135904221533_s1_p0/54519/26233_32397
TGTAAATTGCGGGGAGTTGAGTAAAAGCGTTTTGTAACGTATGGACAGCGTTATGGCGGCATTTG
GTTCTGGGGTTAAGTATTGTTACCAGTTTGCATACACGGGCAAGTTTTCAGGAAGCGACGAGCTG
ATCGTTAAAGGATGCTAACGGCAATCTGCTGGCGACGGCGACAGCGTTATCCATCATTAAAGATCT
GAAGGTGAAAGGTAGCTCT
+
,/-,*..,..,**)..,,./,.++*./.).,+./-%)/////.+$..--.---/.+.-*"/-*-+-,.(/,+,.-/+/...+..////+-//.//,'./-*)./)/#.,),*'-$.)*+).)././*,)-,/)-.(////-.+//,/./,*.,,')/.,-///).,/,)*//((/../**.+.*./.+.,*/,--//,/,-$..-*)///../..,+..#/(//./////-)(/.../-.//-/''/-*./$./,&-.),',+*,)..#%..-*/./..&-.'..,%&.+.(---.+-/)(./)#,/..--.....--/./+-.//)/($*,)*)-$../&)..)'...((/'///#-(//..*,-..+*,.//,,../,/..%-../.*....//*."-./')'-///.-..)//)/-+.)'-&-.--.//-*/,,*.,./,*+.#,.(*//-..//-//-$#%&/$-//)/.-%./*+,/-//-/,//////*-/.*&.../-/*//+/*..)..//.*/*/*//../.'.-.$*+'-)/-././/-..".,*,()//-+//./",,-.(.-/-.-.*//)..,,(./()...+/.+,%.,.)/,.-,//'////)-///..../*/,./'/--,//,.//(-))
ADD REPLYlink modified 9 weeks ago • written 9 weeks ago by 1234anjalianjali12340

With my new command i am getting new error.

command:

canu-1.5/Linux-amd64/bin/canu -correct  -s pacbio.fastq -p ecoli -d ecoli-auto genomeSize=4.8\ pacbio-raw*fastq

output:

ERROR:  File not found or unknown specFile option line '@m141013_011508_sherri_c100709962550000001823135904221533_s1_p0/54519/26233_32397'.
ERROR:  File not found in spec file option '@m141013_011508_sherri_c100709962550000001823135904221533_s1_p0/54519/26233_32397'
ERROR:  File not found or unknown specFile option line CGCCGGGCCTGGGTACGGCTGAAGAAAGATCAGTACGTGGCAAACCAGATATAAAGAAAGCTGCTGACCACGCAAGCACGCTGAACACGGTTAGTTAACTTT
ERROR:  File not found in spec file option 'TGTAAATTGCGGGGAGTTGAGTAAAAGCGTTTTGTAACGTATGGACAGCGTTATGGCGGCATTTGGTTCTGGGGTTAAGTATTGTCAGCCGATTACCAGTTTGC
ERROR:  File not found or unknown specFile option line '+'.
ERROR:  File not found in spec file option '+'
ERROR:  File not found or unknown specFile option line ',/-,*..,..,**)..,,./,.++*./.).,+./-%)/////.+$..--.---/.+.-*"/-*-+-,.(/,+,.-/+/...+..////+-//.//,'./-*)./)/#.,),*'-$.)*+).)././*,)-,/)-.(////-.+//,/./,*.,,')/.,-///).,/,)*//((/../**.+.*./.+.,*/,--//,/,-$..-*)///../..,+..#/(//./////-)(/.../-.//-/''/-*./$./,&-.),',+*,)..#%..-*/./..&-.'..,%&.+.(---.+-/)(./)#,/..--.....--/./+-.//)/($*,)*)-$../&)..)'...((
ERROR:  File not found in spec file option ',/-,*..,..,**)..,,./,.++*./.).,+./-%)/////.+$..--.---/.+.-*"/-*-+-,.(/,+,.-/+/...+..////+-//.//,'./-*)./)/#.,),*'-$.)*+).)././*,)-,/)-.(////-.+//,/./,*.,,')/.,-///).,/,)*//((/../**.+.*./.+.,*/,--//,/,-$..-*)///../..,+..#/(//./////-)(/.../-.//-/''/-*./$./,&-.),',+*,)..#%..-*/./..&-.'..,%&.+.(---.+-/)(./)#,/..--.....--/./+-.//)/($ERROR:  File not found or unknown specFile option line '@m141013_011508_sherri_c100709962550000001823135904221533_s1_p0/54532/2817_4395'.
ERROR:  File not found in spec file option '@m141013_011508_sherri_c100709962550000001823135904221533_s1_p0/54532/2817_4395'
ERROR:  File not found or unknown specFile option line 'CGCTATCACCAGGCTTTTTCCCTTGGCGTGGGCAGCACGAATAAGGCAGGGGTTTTATGACGACGAAAATGGTTTGAGCGTGGCACTCTCCGGCTTCTTGCTGATATGG............

......

ADD REPLYlink modified 9 weeks ago • written 9 weeks ago by 1234anjalianjali12340

The spec file (with -s) is not a fastq file. HAve a look at example .spec files in the docs. Cheers

ADD REPLYlink written 8 weeks ago by colindaven410
0
gravatar for colindaven
9 weeks ago by
colindaven410
colindaven410 wrote:

Canu did not find any reads. Check whether the file ~/canu/pacbio.fastq is present and contains reads.

Also check using "less pacbio.fastq" if the read file is ok and contains enough reads to do the autocorrection step (should be fairly large, i.e. more than 50 -100 MB I would guess).

Good luck.

ADD COMMENTlink written 9 weeks ago by colindaven410
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 607 users visited in the last hour