SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Pacific Biosciences



Similar Threads
Thread Thread Starter Forum Replies Last Post
Wig Bed Overlap Nick Bioinformatics 1 10-10-2014 08:54 AM
Need some suggestion for overlap assembler ljhwahaha Bioinformatics 5 03-25-2014 08:27 AM
Overlap Graph navin_elango Bioinformatics 2 05-20-2013 11:56 AM
overlap -> contig software? gtoal Bioinformatics 0 05-20-2013 10:46 AM
Fragment an assembly for later overlap k-gun12 Bioinformatics 2 05-27-2011 12:16 PM

Reply
 
Thread Tools
Old 07-07-2014, 03:09 AM   #1
Medhat
Member
 
Location: Poland

Join Date: Jun 2013
Posts: 37
Exclamation Overlap prep job failed

I am trying to do the example here

PacBio

specially this type "Correction using complementary data";

I used this commands

Code:
java -jar convertFastaAndQualToFastq.jar pacbio.filtered_subreads.fasta > pacbio.filtered_subreads.fastq

/bin/fastqToCA -libraryname illumina -technology illumina -type sanger -innie -reads illumina.fastq > illumina.frg
/bin/PBcR -length 500 -partitions 200 -l lambdaIll -s pacbio.spec -fastq pacbio.filtered_subreads.fastq genomeSize=50000 illumina.frg
but it gives me this Error

Dumping 7550 fragments from library IID 2
/home/medhat/Downloads/sampleData//templambdaIll/1-overlapper/ovlprep.sh 2
Scanning store to find libraries used and reads to dump.
Added 0 reads to maintain mate relationships.
Dumping 0 fragments from unknown library (version 1 has these)
Dumping 25000 fragments from library IID 1
Dumping 0 fragments from library IID 2
----------------------------------------END CONCURRENT Mon Jul 7 11:48:04 2014 (0 seconds)
ERROR: Overlap prep job /home/medhat/Downloads/sampleData//templambdaIll/1-overlapper/long_reads_part 1 FAILED.

1 overlap partitioning jobs failed.


@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@2

I used totally different machine and it gives me this error

Fragment correction job 0001 failed.
================================================================================

runCA failed.

----------------------------------------
Stack trace:

at /home/medhat/source/celera/wgs-8.2alpha/Linux-amd64/bin/runCA line 1568
main::caFailure('1 overlap jobs failed; remove /home/medhat/lambda-example/sam...', undef) called at /home/medhat/source/celera/wgs-8.2alpha/Linux-amd64/bin/runCA line 4327
main:verlapCorrection() called at /home/medhat/source/celera/wgs-8.2alpha/Linux-amd64/bin/runCA line 6476

----------------------------------------
Failure message:

1 overlap jobs failed; remove /home/medhat/lambda-example/sampleData/lambdaIll/3-overlapcorrection/frgcorr.sh to try again

----------------------------------------END Mon Jul 14 12:36:54 2014 (6 seconds)
Failed to execute /home/medhat/source/celera/wgs-8.2alpha/Linux-amd64/bin/runCA -s /home/medhat/lambda-example/sampleData//templambdaIll/lambdaIll.spec -p asm -d lambdaIll ovlRefBlockLength=100000000000 ovlRefBlockSize=0 useGrid=0 scriptOnGrid=0 unitigger=bogart ovlErrorRate=0.03 utgErrorRate=0.013 cgwErrorRate=0.1 cnsErrorRate=0.1 utgGraphErrorLimit=0 utgGraphErrorRate=0.013 utgMergeErrorLimit=0 utgMergeErrorRate=0.02 frgCorrBatchSize=100000 doOverlapBasedTrimming=1 obtErrorRate=0.03 obtErrorLimit=4.5 frgMinLen=605.993313521545 ovlMinLen=80 consensus=pbutgcns merSize=22 cnsMaxCoverage=1 cnsReuseUnitigs=1 sgePropagateHold="pBcR_asm" lambdaIll.longest25.frg


any idea?

Last edited by Medhat; 07-14-2014 at 03:41 AM.
Medhat is offline   Reply With Quote
Old 07-18-2014, 02:05 AM   #2
wikander
Junior Member
 
Location: Uppsala, Sweden

Join Date: Jul 2014
Posts: 1
Default

Can you paste the contents of your pacbio.spec file?
wikander is offline   Reply With Quote
Old 07-24-2014, 03:44 AM   #3
Medhat
Member
 
Location: Poland

Join Date: Jun 2013
Posts: 37
Default

It is just simple

merSize=14
Medhat is offline   Reply With Quote
Old 12-02-2014, 12:26 PM   #4
cvargasc
Junior Member
 
Location: Valencia, Spain

Join Date: Aug 2011
Posts: 5
Default

Hello,

I'm having the same issue.

Partitioning for hash: 210616671-210877666 ref: 1,210616670
HASH 210616671- 210785993 REFR 1- 210616670 STRINGS 169323 BASES 1000004712 JOB 1
HASH 210785994- 210877666 REFR 1- 210616670 STRINGS 91673 BASES 533324435 JOB 2
----------------------------------------END Tue Dec 2 21:16:07 2014 (55 seconds)
ERROR: Overlap prep job /gnmP04/Soryzae/RawandFilteredReads/Illumina/PE/CA//temppacbio/1-overlapper/long_reads_part 1 FAILED.
ERROR: Overlap prep job /gnmP04/Soryzae/RawandFilteredReads/Illumina/PE/CA//temppacbio/1-overlapper/long_reads_part 2 FAILED.

2 overlap partitioning jobs failed.

I have tried modifying the pacbio.spec file, the only line I'm keeping is merylMemory = 32000.

And the command I'm using is PBcR -length 500 -partitions 200 -libraryname pacbio -threads 20 -s pacbio.spec -fastq allreads.fastq illuminaPE.frg

Does anyone have any idea of what could be the issue?

Thanks in advance!
cvargasc is offline   Reply With Quote
Old 02-04-2015, 02:37 PM   #5
cwisch88
Member
 
Location: St. Louis

Join Date: Jan 2012
Posts: 15
Default

cvargasc,

I wanted to respond to his still unanswered question because I ran into this issue this week.

I discovered that there is an error file in the 1-overlapper directory ending with *.err. Inside it said that the binary sawriter was missing.

I suggest reading through any error files produced in that directory, if there is any mention of sawriter, you'll need to install blasr and perhaps the library it depends on hdf5.

After installing the PBcR stopped giving me that particular error.

Here's the thread that helped me:

http://sourceforge.net/p/wgs-assembl...sage/33127830/
cwisch88 is offline   Reply With Quote
Reply

Tags
assembly, pacbio, sequence

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 11:59 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO