INFO  @ Tue, 05 Sep 2017 14:18:24: 
# Command line: callpeak -t /scratch/AG_Ohler/rebecca/poznan_ChIPseq_analysis/atacseq_chr21/rep1_hg19_Aligned.out.coordSrt.rmdup.uniq.chr21.nameSrt.filtered.1bp.bedpe.bam /scratch/AG_Ohler/rebecca/poznan_ChIPseq_analysis/atacseq_chr21/rep2_hg19_Aligned.out.coordSrt.rmdup.uniq.chr21.nameSrt.filtered.1bp.bedpe.bam -f BAM -g hs -q 0.01 -n /scratch/AG_Ohler/rebecca/poznan_ChIPseq_analysis/atacseq_chr21/testmacs2/macs2_atac --nomodel --shift -75 --extsize 150 --trackline
# ARGUMENTS LIST:
# name = /scratch/AG_Ohler/rebecca/poznan_ChIPseq_analysis/atacseq_chr21/testmacs2/macs2_atac
# format = BAM
# ChIP-seq file = ['/scratch/AG_Ohler/rebecca/poznan_ChIPseq_analysis/atacseq_chr21/rep1_hg19_Aligned.out.coordSrt.rmdup.uniq.chr21.nameSrt.filtered.1bp.bedpe.bam', '/scratch/AG_Ohler/rebecca/poznan_ChIPseq_analysis/atacseq_chr21/rep2_hg19_Aligned.out.coordSrt.rmdup.uniq.chr21.nameSrt.filtered.1bp.bedpe.bam']
# control file = None
# effective genome size = 2.70e+09
# band width = 300
# model fold = [5, 50]
# qvalue cutoff = 1.00e-02
# Larger dataset will be scaled towards smaller dataset.
# Range for calculating regional lambda is: 10000 bps
# Broad region calling is off
 
INFO  @ Tue, 05 Sep 2017 14:18:24: #1 read tag files... 
INFO  @ Tue, 05 Sep 2017 14:18:24: #1 read treatment tags... 
INFO  @ Tue, 05 Sep 2017 14:18:27: #1 tag size is determined as 0 bps 
INFO  @ Tue, 05 Sep 2017 14:18:27: #1 tag size = 0 
INFO  @ Tue, 05 Sep 2017 14:18:27: #1  total tags in treatment: 495785 
INFO  @ Tue, 05 Sep 2017 14:18:27: #1 user defined the maximum tags... 
INFO  @ Tue, 05 Sep 2017 14:18:27: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) 
INFO  @ Tue, 05 Sep 2017 14:18:27: #1  tags after filtering in treatment: 452737 
INFO  @ Tue, 05 Sep 2017 14:18:27: #1  Redundant rate of treatment: 0.09 
INFO  @ Tue, 05 Sep 2017 14:18:27: #1 finished! 
INFO  @ Tue, 05 Sep 2017 14:18:27: #2 Build Peak Model... 
INFO  @ Tue, 05 Sep 2017 14:18:27: #2 Skipped... 
INFO  @ Tue, 05 Sep 2017 14:18:27: #2 Sequencing ends will be shifted towards 5' by 75 bp(s) 
INFO  @ Tue, 05 Sep 2017 14:18:27: #2 Use 150 as fragment length 
INFO  @ Tue, 05 Sep 2017 14:18:27: #3 Call peaks... 
INFO  @ Tue, 05 Sep 2017 14:18:27: #3 Pre-compute pvalue-qvalue table... 
INFO  @ Tue, 05 Sep 2017 14:18:28: #3 Call peaks for each chromosome... 
INFO  @ Tue, 05 Sep 2017 14:18:28: #4 Write output xls file... /scratch/AG_Ohler/rebecca/poznan_ChIPseq_analysis/atacseq_chr21/testmacs2/macs2_atac_peaks.xls 
INFO  @ Tue, 05 Sep 2017 14:18:28: #4 Write peak in narrowPeak format file... /scratch/AG_Ohler/rebecca/poznan_ChIPseq_analysis/atacseq_chr21/testmacs2/macs2_atac_peaks.narrowPeak 
INFO  @ Tue, 05 Sep 2017 14:18:28: #4 Write summits bed file... /scratch/AG_Ohler/rebecca/poznan_ChIPseq_analysis/atacseq_chr21/testmacs2/macs2_atac_summits.bed 
INFO  @ Tue, 05 Sep 2017 14:18:28: Done! 
