Error with Chip-Seq Singularity Pipeline at Macs2 Peak Calling
See original GitHub issueDuring 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:
- Created 4 years ago
- Comments:5 (3 by maintainers)
Top 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 >Top Related Medium Post
No results found
Top Related StackOverflow Question
No results found
Troubleshoot Live Code
Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start FreeTop Related Reddit Thread
No results found
Top Related Hackernoon Post
No results found
Top Related Tweet
No results found
Top Related Dev.to Post
No results found
Top Related Hashnode Post
No results found
Top GitHub Comments
@tiagobrc : TF ChIP-seq cannot run without inputs but histone ChIP-seq can do .
Use : instead of =
“chip.fraglen”: [200, 150],