[85/89f50f] process > get_software_versions ...
# clipseq
l
[85/89f50f] process > get_software_versions [100%] 1 of 1 ✔️ [b0/f34aa6] process > generate_premap_index (Homo_sapiens.smallRNA.fa.gz) [100%] 1 of 1 ✔️ [d7/47f7cd] process > generate_fai (GRCh38.p13.genome.fa) [100%] 1 of 1 ✔️ [3a/dcce0b] process > generate_star_index (GRCh38.p13.genome.fa) [100%] 1 of 1 ✔️ [79/306e25] process > icount_segment (gencode.v34.annotation.gtf) [100%] 1 of 1 ✔️ [64/e3a3f6] process > fastqc (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [3f/e98335] process > move_umi (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [fa/0cde74] process > cutadapt (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [b2/a60dcb] process > premap (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [dd/843fe1] process > align (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [40/8114ad] process > preseq (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [e9/03c5ee] process > dedup (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [3d/41d84f] process > rseqc (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [48/db8344] process > get_crosslinks (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [6a/f7c583] process > icount_peak_call (PTBP1iCLIP_siCtrl_repl1_s) [ 25%] 1 of 4 [68/f63dbe] process > paraclu_peak_call (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [17/05d893] process > pureclip_peak_call (PTBP1iCLIP_siMATR3_rep2_s) [ 57%] 4 of 7, failed: 4, retries: 3 [d3/d99fb6] process > piranha_peak_call (PTBP1iCLIP_siCtrl_repl1_s) [100%] 4 of 4 ✔️ [- ] process > clipqc [ 0%] 0 of 1 [- ] process > multiqc - [ac/bbfcc0] process > output_documentation [100%] 1 of 1 ✔️ Error executing process > 'pureclip_peak_call (PTBP1iCLIP_siMATR3_rep2_s)' Caused by: Process
pureclip_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&lt;boostexception detailerror info injector&lt;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"