Skip to content
New issue

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

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

Already on GitHub? Sign in 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.

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

No branches or pull requests

2 participants
@niehu2018 and others