|
4
|
1 #peaks file
|
|
|
2 # This file is generated by MACS version 2.1.2
|
|
|
3 # Command line: callpeak -t /tmp/tmpKNHGvp/files/000/dataset_2.dat -c /tmp/tmpKNHGvp/files/000/dataset_1.dat --format=BED --name=test_MACS2.1.2 --bw=300 --gsize=775000000.0 --nomodel --extsize=243 --qvalue=0.05 -B --SPMR --mfold 5 50 --keep-dup 1
|
|
|
4 # ARGUMENTS LIST:
|
|
|
5 # name = test_MACS2.1.2
|
|
|
6 # format = BED
|
|
|
7 # ChIP-seq file = ['/tmp/tmpKNHGvp/files/000/dataset_2.dat']
|
|
|
8 # control file = ['/tmp/tmpKNHGvp/files/000/dataset_1.dat']
|
|
|
9 # effective genome size = 7.75e+08
|
|
|
10 # band width = 300
|
|
|
11 # model fold = [5, 50]
|
|
|
12 # qvalue cutoff = 5.00e-02
|
|
|
13 # The maximum gap between significant sites is assigned as the read length/tag size.
|
|
|
14 # The minimum length of peaks is assigned as the predicted fragment length "d".
|
|
|
15 # Larger dataset will be scaled towards smaller dataset.
|
|
|
16 # Range for calculating regional lambda is: 1000 bps and 10000 bps
|
|
|
17 # Broad region calling is off
|
|
|
18 # Paired-End mode is off
|
|
|
19 # MACS will save fragment pileup signal per million reads
|
|
|
20
|
|
|
21 # tag size is determined as 50 bps
|
|
|
22 # total tags in treatment: 50
|
|
|
23 # tags after filtering in treatment: 50
|
|
|
24 # maximum duplicate tags at the same position in treatment = 1
|
|
|
25 # Redundant rate in treatment: 0.00
|
|
|
26 # total tags in control: 50
|
|
|
27 # tags after filtering in control: 50
|
|
|
28 # maximum duplicate tags at the same position in control = 1
|
|
|
29 # Redundant rate in control: 0.00
|
|
|
30 # d = 243
|
|
|
31 #chr start end length abs_summit pileup -log10(pvalue) fold_enrichment -log10(qvalue) name
|
|
|
32 chr26 4118914 4119282 369 4119130 9.00 9.13132 6.31632 2.51561 test_MACS2.1.2_peak_1
|
