abyss mpirun non zero code
1
0
Entering edit mode
10.2 years ago
Siva ▴ 20

Hi guys,

I've a problem with mpirun and abyss :(

I've tried to assemble test reads in abyss directory:

mpirun -np 3 abyss-pe contigs clean name=test mp="" j=3 k=20 in='test-data/reads1.fastq test-data/reads2.fastq'

And I got output file ".fa", ".hist" ... But also an error :

-------------------------------------------------------
**Primary job terminated normally, but 1 process returned
a non-zero exit code**.. Per user-direction, the job has been aborted.
-------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:
  Process name: [[25897,1],1]
**  Exit code:    2**
-------------------------------------------------------*

In that example I got outputs but with my own data I haven't got any even if I can see in the log make: *** [mappedR.fa] Completed.

Advance thanks,
Siva

Assembly cluster abyss mpirun • 4.9k views
ADD COMMENT
0
Entering edit mode

Even without mpirun i got that: "make: *** [mappedR_XRQrefSF099_SEPE-1.fa] Processus arrete"

I've got 105G data (fastq file)

ADD REPLY
0
Entering edit mode
10.2 years ago
benv ▴ 730

Hi Siva,

Just a note: The second way you tried to run abyss-pe (without the mpirun prefix) is the correct way, i.e.:

$ abyss-pe contigs clean name=test mp="" j=3 k=20 in='test-data/reads1.fastq test-data/reads2.fastq'

The mp="" part is probably not needed.

Your abyss-pe command looks fine. Unfortunately, the output you have provided does not provide any indication about what is causing the error. Can you please post the full log of your job to a github gist?

ADD COMMENT
0
Entering edit mode

Also: If possible, please generate the log with the verbose option turned on, by adding the argument "v=-v" to the abyss-pe command line.

ADD REPLY
0
Entering edit mode

Hi Benv,

Actually I haven't got more informations, I've tried with verbose mode and with my real data and have this file.

ABYSS -k60 -q3 -v --coverage-hist=coverage.hist -s mappedR_XRQrefSF099_SEPE-bubbles.fa -o mappedR_XRQrefSF099_SEPE-1.fa final_bwa_refXRQ_SF099_primALN_mapped_PE.fq final_bwa_refXRQ_SF099_primALN_mapped_SE.fq
ABySS 1.5.2
Reading `final_bwa_refXRQ_SF099_primALN_mapped_PE.fq'...
Read 100000 reads. Hash load: 3812904 / 1073741824 = 0.00355 using 482 MB
Read 5100000 reads. Hash load: 156329408 / 1073741824 = 0.146 using 4.72 GB
...
Read 4800000 reads. Hash load: 148139019 / 1073741824 = 0.138 using 4.55 GB
Read 4800000 reads. Hash load: 148139019 / 1073741824 = 0.138 using 4.55 GB
Read 4500000 reads. Hash load: 139824327 / 1073741824 = 0.13 using 4.36 GB
Read 4800000 reads. Hash load: 148139019 / 1073741824 = 0.138 using 4.55 GB
Read 4300000 reads. Hash load: 134225708 / 1073741824 = 0.125 using 4.23 GB
-------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
Read 4800000 reads. Hash load: 148139019 / 1073741824 = 0.138 using 4.55 GB
Read 4900000 reads. Hash load: 150877230 / 1073741824 = 0.141 using 4.6 GB
Read 4900000 reads. Hash load: 150877230 / 1073741824 = 0.141 using 4.6 GB
Read 5000000 reads. Hash load: 153597524 / 1073741824 = 0.143 using 4.66 GB
Read 5200000 reads. Hash load: 159060550 / 1073741824 = 0.148 using 4.78 GB
Read 5300000 reads. Hash load: 161777571 / 1073741824 = 0.151 using 4.84 GB
Read 4900000 reads. Hash load: 150877230 / 1073741824 = 0.141 using 4.6 GB
Read 4400000 reads. Hash load: 137029057 / 1073741824 = 0.128 using 4.3 GB
Read 4900000 reads. Hash load: 150877230 / 1073741824 = 0.141 using 4.6 GB
Read 4900000 reads. Hash load: 150877230 / 1073741824 = 0.141 using 4.6 GB
Read 4600000 reads. Hash load: 142605340 / 1073741824 = 0.133 using 4.42 GB
Read 3900000 reads. Hash load: 122979009 / 1073741824 = 0.115 using 3.95 GB
Read 5100000 reads. Hash load: 156329408 / 1073741824 = 0.146 using 4.72 GB
Read 5200000 reads. Hash load: 159060550 / 1073741824 = 0.148 using 4.78 GB
Read 4700000 reads. Hash load: 145375696 / 1073741824 = 0.135 using 4.48 GB
Read 4800000 reads. Hash load: 148139019 / 1073741824 = 0.138 using 4.55 GB
Read 4500000 reads. Hash load: 139824327 / 1073741824 = 0.13 using 4.36 GB
Read 4600000 reads. Hash load: 142605340 / 1073741824 = 0.133 using 4.42 GB
Read 4700000 reads. Hash load: 145375696 / 1073741824 = 0.135 using 4.48 GB
Read 4800000 reads. Hash load: 148139019 / 1073741824 = 0.138 using 4.55 GB
Read 5000000 reads. Hash load: 153597524 / 1073741824 = 0.143 using 4.66 GB
Read 5100000 reads. Hash load: 156329408 / 1073741824 = 0.146 using 4.72 GB
Read 5200000 reads. Hash load: 159060550 / 1073741824 = 0.148 using 4.78 GB
Read 5000000 reads. Hash load: 153597524 / 1073741824 = 0.143 using 4.66 GB
Read 5100000 reads. Hash load: 156329408 / 1073741824 = 0.146 using 4.72 GB
Read 5000000 reads. Hash load: 153597524 / 1073741824 = 0.143 using 4.66 GB
Read 5100000 reads. Hash load: 156329408 / 1073741824 = 0.146 using 4.72 GB
Read 5200000 reads. Hash load: 159060550 / 1073741824 = 0.148 using 4.78 GB
Read 5000000 reads. Hash load: 153597524 / 1073741824 = 0.143 using 4.66 GB
Read 5100000 reads. Hash load: 156329408 / 1073741824 = 0.146 using 4.72 GB
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:
Process name: [[25065,1],5]
Exit code: 2
--------------------------------------------------------------------------
view raw log.txt hosted with ❤ by GitHub

ps: I use mp="" to avoid other assembly like scaff or unitigs, I just want to have contigs.

ADD REPLY
0
Entering edit mode

I think I've maybe found the error:

error: sequence and quality must be the same length near
ANCGTTTTCTCGGTCGACATGGAGGTAAACAAAGGAACACTAGATCCGAAGATCGGATGGTCTTGCATATTATTAAAGGGTGTTAACAAGTACGAATACAA
@#1==DDAFFFD=EFFB1CGEEFECFFGIIIIIEGII?DFECEEFDEAFEEII;?A1([M::main_bam2fq] discarded 8840174 singletons
make: *** [mappedR_XRQrefSF099_SEPE-1.fa] Erreur 1
ADD REPLY

Login before adding your answer.

Traffic: 2628 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6