-
Notifications
You must be signed in to change notification settings - Fork 20
/
configfile_BiasCorrection_ICEBias
116 lines (88 loc) · 3.55 KB
/
configfile_BiasCorrection_ICEBias
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
#====================================
# Sample configuration file for running FitHiChIP
#====================================
##********
##***** parameters to provide the input HiChIP alignment files
##********
##============
## option 1: provide the valid pairs from HiC-Pro pipeline - can be gzipped as well
##============
ValidPairs=./TestData/Sample_ValidPairs.txt.gz
##============
## option 2: provide the bin interval and contact matrix files from HiC-Pro pipeline
##============
## bin interval file (of the format *_abs.bed from HiC-pro output)
Interval=
## matrix file (of the format *.matrix from HiC-pro output)
Matrix=
##============
## option 3: If HiChIP was processed by aligners other than HiC-Pro
## a) provide the locus pairs as a .bed formatted file with the following format (7 fields):
## chr1 start1 end1 chr2 start2 end2 contactcounts
##============
Bed=
##============
## option 4: If HiChIP data is provided in .hic format
## Make sure that the .hic file contains the target resolution which is provided in the BINSIZE parameter (below)
##============
HIC=
##============
## option 5: If HiChIP data is provided in .cool / .mcool format
## Make sure that the .cool or .mcool file contains the target resolution which is provided in the BINSIZE parameter (below)
##============
COOL=
##********
## File containing chromomosome size information corresponding to the reference genome.
##********
ChrSizeFile=./TestData/chrom_hg19.sizes
##********
## Mandatory parameter - Reference ChIP-seq / HiChIP peaks (in .bed format) - can be gzipped as well
## We recommend using reference ChIP-seq peaks (if available)
## Otherwise, peaks can be computed from HiChIP data.
## See the documentation: https://ay-lab.github.io/FitHiChIP/usage/Utilities.html#inferring-peaks-from-hichip-data-for-use-in-the-hichip-pipeline
##********
PeakFile=./TestData/Sample.Peaks.gz
##********
## Mandatory parameter - Output directory to contain all the results
##********
OutDir=./TestData/results/
##********
## Mandatory parameter - Boolean variable indicating if the reference genome is circular
## 0, by default. If 1 (circular genome), calculation of genomic distance is slightly different
##********
CircularGenome=0
##********
##***** Various FitHiChIP loop calling related parameters
##********
##Interaction type
## 1: peak to peak
## 2: peak to non peak
## 3: peak to all (default - both peak-to-peak and peak-to-nonpeak)
## 4: all to all (similar to Hi-C)
## 5: All of the modes 1 to 4 are computed.
IntType=3
## Bin size, in bases, for the interactions. Default = 5000 (5 Kb).
BINSIZE=5000
## Lower distance threshold of loops - default = 20000 (20 Kb)
LowDistThr=20000
## Upper distance threshold of loops - default = 2000000 (2 Mb)
UppDistThr=2000000
## Values 0/1 - Applicable if IntType = 3 (peak to all output interactions)
## 1 indicates FitHiChIP(S) model - uses only peak to peak loops for background modeling
## 0 corresponds to FitHiChIP(L) - uses both peak to peak and peak to nonpeak loops for background modeling
UseP2PBackgrnd=0
## type of bias - values: 1 / 2
## 1: coverage bias regression
## 2: ICE bias regression
BiasType=2
## if 1 (default), merge filtering (corresponding to either FitHiChIP(L+M) or FitHiChIP(S+M)
## depending on the parameter UseP2PBackgrnd) is enabled
MergeInt=1
## FDR (q-value) threshold for loop significance
QVALUE=0.01
## prefix string of all the output files (Default = 'FitHiChIP').
PREFIX=FitHiChIP
## Binary variable 1/0:
## if 1, overwrites any existing output file.
## otherwise (0), does not overwrite any output file.
OverWrite=0