question-mark
Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

Error with Chip-Seq Singularity Pipeline at Macs2 Peak Calling

See original GitHub issue

During MACS2 processing with Singularity Chip-Seq pipeline, the pipeline throws up an error. I Tried but couldn’t identify the origin of the error.

Can you help me? What am I missing?

I appreciate a lot thanks!

singularity --version
2.6.1-dist

Input.json:

{
    "chip.genome_tsv" : "mm10.tsv",
    "chip.paired_end" : false,

    "chip.fastqs_rep1_R1" : [ "ChIP_R1.fastq.gz"],
    "chip.ctl_fastqs_rep1_R1" : [ "ChIP_Input_R1.fastq.gz" ],

    "chip.title" : "ChIP",
    "chip.description" : "ChIP Project",

    "chip.pipeline_type" : "tf",
    "chip.peak_caller" : "spp",

    "chip.align_only" : false,
    "chip.true_rep_only" : false,

    "chip.disable_fingerprint" : false,
    "chip.enable_count_signal_track" : false,

    "chip.xcor_pe_trim_bp" : 50,

    "chip.dup_marker" : "picard",
    "chip.mapq_thresh" : 30,
    "chip.no_dup_removal" : false,

    "chip.mito_chr_name" : "chrM",
    "chip.regex_filter_reads" : "chrM",
    "chip.subsample_reads" : 0,
    "chip.ctl_subsample_reads" : 0,
    "chip.xcor_subsample_reads" : 15000000,

    "chip.keep_irregular_chr_in_bfilt_peak" : false,
    
    "chip.always_use_pooled_ctl" : true,
    "chip.ctl_depth_ratio" : 1.2,

    "chip.macs2_cap_num_peak" : 500000,
    "chip.pval_thresh" : 0.01,
    "chip.idr_thresh" : 0.05,
    "chip.spp_cap_num_peak" : 300000,

    "chip.bwa_cpu" : 12,
    "chip.bwa_mem_mb" : 200000,
    "chip.bwa_time_hr" : 480,
    "chip.bwa_disks" : "local-disk 100 HDD",

    "chip.filter_cpu" : 12,
    "chip.filter_mem_mb" : 20000,
    "chip.filter_time_hr" : 24,
    "chip.filter_disks" : "local-disk 100 HDD",

    "chip.bam2ta_cpu" : 12,
    "chip.bam2ta_mem_mb" : 100000,
    "chip.bam2ta_time_hr" : 6,
    "chip.bam2ta_disks" : "local-disk 100 HDD",

    "chip.spr_mem_mb" : 16000,

    "chip.fingerprint_cpu" : 12,
    "chip.fingerprint_mem_mb" : 12000,
    "chip.fingerprint_time_hr" : 6,
    "chip.fingerprint_disks" : "local-disk 100 HDD",

    "chip.xcor_cpu" : 12,
    "chip.xcor_mem_mb" : 16000,
    "chip.xcor_time_hr" : 24,
    "chip.xcor_disks" : "local-disk 100 HDD",

    "chip.macs2_mem_mb" : 16000,
    "chip.macs2_time_hr" : 24,
    "chip.macs2_disks" : "local-disk 100 HDD",

    "chip.spp_cpu" : 12,
    "chip.spp_mem_mb" : 16000,
    "chip.spp_time_hr" : 72,
    "chip.spp_disks" : "local-disk 100 HDD"
}
find -name stderr -not -empty | xargs tail -n +1
Traceback (most recent call last):
  File "/software/chip-seq-pipeline/src/encode_macs2_chip.py", line 225, in <module>
    main()
  File "/software/chip-seq-pipeline/src/encode_macs2_chip.py", line 200, in main
    args.out_dir)
  File "/software/chip-seq-pipeline/src/encode_macs2_chip.py", line 88, in macs2
    run_shell_cmd(cmd0)
  File "/software/chip-seq-pipeline/src/encode_common.py", line 252, in run_shell_cmd
    raise Exception(err_str)
Exception: PID=28849, PGID=28849, RC=1
STDERR=ERROR:root:--extsize must >= 1!
STDOUT=

Issue Analytics

  • State:closed
  • Created 4 years ago
  • Comments:5 (3 by maintainers)

github_iconTop GitHub Comments

1reaction
leepc12commented, May 30, 2019

@tiagobrc : TF ChIP-seq cannot run without inputs but histone ChIP-seq can do .

1reaction
leepc12commented, May 24, 2019

Use : instead of =

“chip.fraglen”: [200, 150],

Read more comments on GitHub >

github_iconTop Results From Across the Web

nf-core/ chipseq
For other genomes, if this parameter is not specified then the MACS2 peak-calling and differential analysis will be skipped.
Read more >
Peak calling with MACS2 - GitHub Pages
A commonly used tool for identifying transcription factor binding sites is named Model-based Analysis of ChIP-seq (MACS). The MACS algorithm captures the ...
Read more >
ChIP Sequencing Pipeline — GenPipes 4.3.2 documentation
Peaks are called using Model based Analysis for Chip Sequencing (MACS2) and annotated using Homer. ... The ChIP-Seq pipeline can also be used...
Read more >
ChIP-AP: an integrated analysis pipeline for unbiased ... - NCBI
Of all the programs required for a ChIP-seq analysis, the choice of peak caller is the most critical [7]. In 2016, Steinhauser et...
Read more >
https://dockstore.org/api/ga4gh/trs/v2/tools/%23wo...
TF ChIP-seq requires control for peak calling but histone ChIP-seq does not. ... String docker = 'encodedcc/chip-seq-pipeline:v2.1.0' String singularity ...
Read more >

github_iconTop Related Medium Post

No results found

github_iconTop Related StackOverflow Question

No results found

github_iconTroubleshoot Live Code

Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start Free

github_iconTop Related Reddit Thread

No results found

github_iconTop Related Hackernoon Post

No results found

github_iconTop Related Tweet

No results found

github_iconTop Related Dev.to Post

No results found

github_iconTop Related Hashnode Post

No results found