Li Las
04/18/2025, 9:17 AMpureclip_peak_call (PTBP1iCLIP_siMATR3_rep2_s) terminated with an error exit status (134)
Command executed:
pureclip \
-i PTBP1iCLIP_siMATR3_rep2_s.dedup.bam \
-bai PTBP1iCLIP_siMATR3_rep2_s.dedup.bam.bai \
-g GRCh38.p13.genome.fa \
-nt 20 \
-bc 0 -dm 8 -iv 'all' \
-o "PTBP1iCLIP_siMATR3_rep2_s.sigxl.bed" \
-or "PTBP1iCLIP_siMATR3_rep2_s.8nt.peaks.bed"
pigz PTBP1iCLIP_siMATR3_rep2_s.sigxl.bed PTBP1iCLIP_siMATR3_rep2_s.8nt.peaks.bed
Command exit status:
134
Command output:
b0:0
mean:1
k:1
tp:0.00807885
initial transition frequencies:
0: 1 1 1 1
1: 1 1 1 1
2: 1 1 1 1
3: 1 1 1 1
Build HMM ...
* Transition probabilitites *
0: 0.25 0.25 0.25 0.25
1: 0.25 0.25 0.25 0.25
2: 0.25 0.25 0.25 0.25
3: 0.25 0.25 0.25 0.25
Baum-Welch ...
learn binomial parameter
.. 0th iteration
computeEmissionProbs()
computeStatePosteriorsFB()
updateDensityParams()
** Convergence ! **
learn gamma parameter
.. 0th iteration
computeEmissionProbs()
computeStatePosteriorsFB()
updateDensityParams()
** Convergence ! **
learn binomial parameter
.. 0th iteration
computeEmissionProbs()
computeStatePosteriorsFB()
updateDensityParams()
** Convergence ! **
learn gamma parameter
.. 0th iteration
computeEmissionProbs()
computeStatePosteriorsFB()
updateDensityParams()
** Convergence ! **
Apply learned parameters to whole dataset ...
chr1
Get preprocessed intervals:
Compute KDEs ...
Estiamte Ns ... 1
Apply learned HMM:
applyParameters
Command error:
INFO: Environment variable SINGULARITYENV_TMPDIR is set, but APPTAINERENV_TMPDIR is preferred
INFO: Environment variable SINGULARITYENV_NXF_DEBUG is set, but APPTAINERENV_NXF_DEBUG is preferred
ERROR: Can't find contigName in store.contigNameStore!all
terminate called after throwing an instance of 'boost:exception detailclone impl<boostexception detailerror info injector<std:domain_error> >'
what(): Error in function boost:math:binomial_distribution<long double>:binomial distribution Success fraction argument is -nan, but must be >= 0 and <= 1 !
.command.sh: line 9: 226464 Aborted (core dumped) pureclip -i PTBP1iCLIP_siMATR3_rep2_s.dedup.bam -bai PTBP1iCLIP_siMATR3_rep2_s.dedup.bam.bai -g GRCh38.p13.genome.fa -nt 20 -bc 0 -dm 8 -iv 'all' -o "PTBP1iCLIP_siMATR3_rep2_s.sigxl.bed" -or "PTBP1iCLIP_siMATR3_rep2_s.8nt.peaks.bed"