Skip to content
New issue

Have a question about this project? # for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “#”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? # to your account

Segmentation fault (core dumped) #202

Open
niehu2018 opened this issue Jul 25, 2017 · 1 comment
Open

Segmentation fault (core dumped) #202

niehu2018 opened this issue Jul 25, 2017 · 1 comment

Comments

@niehu2018
Copy link

$ macs2 callpeak -t /public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/align/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.tagAlign.gz -f BED -n "/public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/peak/macs2/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.pf" -g "mm" -p 0.01 --nomodel --shift -75 --extsize 150 -B --SPMR --keep-dup all --call-summits
INFO @ Tue, 25 Jul 2017 12:33:24:

Command line: callpeak -t /public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/align/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.tagAlign.gz -f BED -n /public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/peak/macs2/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.pf -g mm -p 0.01 --nomodel --shift -75 --extsize 150 -B --SPMR --keep-dup all --call-summits

ARGUMENTS LIST:

name = /public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/peak/macs2/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.pf

format = BED

ChIP-seq file = ['/public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/align/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.tagAlign.gz']

control file = None

effective genome size = 1.87e+09

band width = 300

model fold = [5, 50]

pvalue cutoff = 1.00e-02

qvalue will not be calculated and reported as -1 in the final output.

Larger dataset will be scaled towards smaller dataset.

Range for calculating regional lambda is: 10000 bps

Broad region calling is off

Paired-End mode is off

Searching for subpeak summits is on

MACS will save fragment pileup signal per million reads

INFO @ Tue, 25 Jul 2017 12:33:24: #1 read tag files...
INFO @ Tue, 25 Jul 2017 12:33:24: #1 read treatment tags...
INFO @ Tue, 25 Jul 2017 12:33:28: 1000000
INFO @ Tue, 25 Jul 2017 12:33:31: 2000000
INFO @ Tue, 25 Jul 2017 12:33:35: 3000000
INFO @ Tue, 25 Jul 2017 12:33:38: 4000000
INFO @ Tue, 25 Jul 2017 12:33:42: 5000000
INFO @ Tue, 25 Jul 2017 12:33:45: 6000000
INFO @ Tue, 25 Jul 2017 12:33:49: 7000000
INFO @ Tue, 25 Jul 2017 12:33:52: 8000000
INFO @ Tue, 25 Jul 2017 12:33:55: 9000000
INFO @ Tue, 25 Jul 2017 12:33:58: 10000000
INFO @ Tue, 25 Jul 2017 12:34:01: 11000000
INFO @ Tue, 25 Jul 2017 12:34:04: 12000000
INFO @ Tue, 25 Jul 2017 12:34:07: 13000000
INFO @ Tue, 25 Jul 2017 12:34:09: 14000000
INFO @ Tue, 25 Jul 2017 12:34:12: 15000000
INFO @ Tue, 25 Jul 2017 12:34:15: 16000000
INFO @ Tue, 25 Jul 2017 12:34:18: 17000000
INFO @ Tue, 25 Jul 2017 12:34:21: 18000000
INFO @ Tue, 25 Jul 2017 12:34:24: 19000000
INFO @ Tue, 25 Jul 2017 12:34:27: 20000000
INFO @ Tue, 25 Jul 2017 12:34:29: 21000000
INFO @ Tue, 25 Jul 2017 12:34:32: 22000000
INFO @ Tue, 25 Jul 2017 12:34:36: #1 tag size is determined as 40 bps
INFO @ Tue, 25 Jul 2017 12:34:36: #1 tag size = 40
INFO @ Tue, 25 Jul 2017 12:34:36: #1 total tags in treatment: 22535426
INFO @ Tue, 25 Jul 2017 12:34:36: #1 finished!
INFO @ Tue, 25 Jul 2017 12:34:36: #2 Build Peak Model...
INFO @ Tue, 25 Jul 2017 12:34:36: #2 Skipped...
INFO @ Tue, 25 Jul 2017 12:34:36: #2 Sequencing ends will be shifted towards 5' by 75 bp(s)
INFO @ Tue, 25 Jul 2017 12:34:36: #2 Use 150 as fragment length
INFO @ Tue, 25 Jul 2017 12:34:36: #3 Call peaks...
INFO @ Tue, 25 Jul 2017 12:34:36: #3 Going to call summits inside each peak ...
INFO @ Tue, 25 Jul 2017 12:34:36: #3 Call peaks with given -log10pvalue cutoff: 2.00000 ...
INFO @ Tue, 25 Jul 2017 12:34:36: #3 Pre-compute pvalue-qvalue table...
INFO @ Tue, 25 Jul 2017 12:35:46: #3 In the peak calling step, the following will be performed simultaneously:
INFO @ Tue, 25 Jul 2017 12:35:46: #3 Write bedGraph files for treatment pileup (after scaling if necessary)... /public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/peak/macs2/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.pf_treat_pileup.bdg
INFO @ Tue, 25 Jul 2017 12:35:46: #3 Write bedGraph files for control lambda (after scaling if necessary)... /public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/peak/macs2/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.pf_control_lambda.bdg
INFO @ Tue, 25 Jul 2017 12:35:46: #3 --SPMR is requested, so pileup will be normalized by sequencing depth in million reads.
INFO @ Tue, 25 Jul 2017 12:35:46: #3 Call peaks for each chromosome...
INFO @ Tue, 25 Jul 2017 12:40:39: #4 Write output xls file... /public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/peak/macs2/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.pf_peaks.xls
INFO @ Tue, 25 Jul 2017 12:40:41: #4 Write peak in narrowPeak format file... /public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/peak/macs2/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.pf_peaks.narrowPeak
INFO @ Tue, 25 Jul 2017 12:40:43: #4 Write summits bed file... /public/home/niehu/projects/Lulab20170712/pipeline/E9_5/out/peak/macs2/rep2/E9_5_2_R1.trim.PE2SE.nodup.tn5.pf_summits.bed
INFO @ Tue, 25 Jul 2017 12:40:44: Done!
Segmentation fault (core dumped)

How to fix this problem.

# for free to join this conversation on GitHub. Already have an account? # to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants
@niehu2018 and others