annotate lib/strelka_config_isaac_default.ini @ 19:bddc6981a6a7

change version
author mini
date Wed, 15 Oct 2014 15:30:28 +0200
parents 3c10d88b55ad
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
18
3c10d88b55ad improved user interface
mini
parents:
diff changeset
1
3c10d88b55ad improved user interface
mini
parents:
diff changeset
2 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
3 ; User configuration options for Strelka somatic small-variant caller
3c10d88b55ad improved user interface
mini
parents:
diff changeset
4 ; workflow:
3c10d88b55ad improved user interface
mini
parents:
diff changeset
5 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
6
3c10d88b55ad improved user interface
mini
parents:
diff changeset
7 [user]
3c10d88b55ad improved user interface
mini
parents:
diff changeset
8
3c10d88b55ad improved user interface
mini
parents:
diff changeset
9 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
10 ; isSkipDepthFilters should be set to 1 to skip depth filtration for
3c10d88b55ad improved user interface
mini
parents:
diff changeset
11 ; whole exome or other targeted sequencing data
3c10d88b55ad improved user interface
mini
parents:
diff changeset
12 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
13 isSkipDepthFilters = 0
3c10d88b55ad improved user interface
mini
parents:
diff changeset
14
3c10d88b55ad improved user interface
mini
parents:
diff changeset
15 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
16 ; strelka will not accept input reads above this depth (they will be skipped
3c10d88b55ad improved user interface
mini
parents:
diff changeset
17 ; until the depth drops below this value). Set this value <= 0 to disable
3c10d88b55ad improved user interface
mini
parents:
diff changeset
18 ; this feature. Using this filter will bound memory usage given extremely high
3c10d88b55ad improved user interface
mini
parents:
diff changeset
19 ; depth input, but may be problematic in high-depth targeted sequencing
3c10d88b55ad improved user interface
mini
parents:
diff changeset
20 ; applications.
3c10d88b55ad improved user interface
mini
parents:
diff changeset
21 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
22 maxInputDepth = 10000
3c10d88b55ad improved user interface
mini
parents:
diff changeset
23
3c10d88b55ad improved user interface
mini
parents:
diff changeset
24 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
25 ; If the depth filter is not skipped, all variants which occur at a
3c10d88b55ad improved user interface
mini
parents:
diff changeset
26 ; depth greater than depthFilterMultiple*chromosome mean depth will be
3c10d88b55ad improved user interface
mini
parents:
diff changeset
27 ; filtered out.
3c10d88b55ad improved user interface
mini
parents:
diff changeset
28 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
29 depthFilterMultiple = 3.0
3c10d88b55ad improved user interface
mini
parents:
diff changeset
30
3c10d88b55ad improved user interface
mini
parents:
diff changeset
31 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
32 ; Somatic SNV calls are filtered at sites where greater than this
3c10d88b55ad improved user interface
mini
parents:
diff changeset
33 ; fraction of basecalls have been removed by the mismatch density
3c10d88b55ad improved user interface
mini
parents:
diff changeset
34 ; filter in either sample.
3c10d88b55ad improved user interface
mini
parents:
diff changeset
35 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
36 snvMaxFilteredBasecallFrac = 0.4
3c10d88b55ad improved user interface
mini
parents:
diff changeset
37
3c10d88b55ad improved user interface
mini
parents:
diff changeset
38 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
39 ; Somatic SNV calls are filtered at sites where greater than this
3c10d88b55ad improved user interface
mini
parents:
diff changeset
40 ; fraction of overlapping reads contain deletions which span the SNV
3c10d88b55ad improved user interface
mini
parents:
diff changeset
41 ; call site.
3c10d88b55ad improved user interface
mini
parents:
diff changeset
42 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
43 snvMaxSpanningDeletionFrac = 0.75
3c10d88b55ad improved user interface
mini
parents:
diff changeset
44
3c10d88b55ad improved user interface
mini
parents:
diff changeset
45 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
46 ; Somatic indel calls are filtered if they represent an expansion or
3c10d88b55ad improved user interface
mini
parents:
diff changeset
47 ; contraction of a repeated pattern with a repeat count greater than
3c10d88b55ad improved user interface
mini
parents:
diff changeset
48 ; indelMaxRefRepeat in the reference (ie. if indelMaxRefRepeat is 8,
3c10d88b55ad improved user interface
mini
parents:
diff changeset
49 ; then the indel is filtered when it is an expansion/contraction of a
3c10d88b55ad improved user interface
mini
parents:
diff changeset
50 ; homopolymer longer than 8 bases, a dinucleotide repeat longer than
3c10d88b55ad improved user interface
mini
parents:
diff changeset
51 ; 16 bases, etc.)
3c10d88b55ad improved user interface
mini
parents:
diff changeset
52 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
53 indelMaxRefRepeat = 8
3c10d88b55ad improved user interface
mini
parents:
diff changeset
54
3c10d88b55ad improved user interface
mini
parents:
diff changeset
55 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
56 ; Somatic indel calls are filtered if greater than this fraction of
3c10d88b55ad improved user interface
mini
parents:
diff changeset
57 ; basecalls in a window extending 50 bases to each side of an indel's
3c10d88b55ad improved user interface
mini
parents:
diff changeset
58 ; call position have been removed by the mismatch density filter.
3c10d88b55ad improved user interface
mini
parents:
diff changeset
59 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
60 indelMaxWindowFilteredBasecallFrac = 0.3
3c10d88b55ad improved user interface
mini
parents:
diff changeset
61
3c10d88b55ad improved user interface
mini
parents:
diff changeset
62 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
63 ; Somatic indels are filtered if they overlap ’interrupted
3c10d88b55ad improved user interface
mini
parents:
diff changeset
64 ; homopolymers’ greater than this length. The term 'interrupted
3c10d88b55ad improved user interface
mini
parents:
diff changeset
65 ; homopolymer' is used to indicate the longest homopolymer which can
3c10d88b55ad improved user interface
mini
parents:
diff changeset
66 ; be found intersecting or adjacent to the called indel when a single
3c10d88b55ad improved user interface
mini
parents:
diff changeset
67 ; non-homopolymer base is allowed.
3c10d88b55ad improved user interface
mini
parents:
diff changeset
68 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
69 indelMaxIntHpolLength = 14
3c10d88b55ad improved user interface
mini
parents:
diff changeset
70
3c10d88b55ad improved user interface
mini
parents:
diff changeset
71 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
72 ; prior probability of a somatic snv or indel
3c10d88b55ad improved user interface
mini
parents:
diff changeset
73 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
74 ssnvPrior = 0.000001
3c10d88b55ad improved user interface
mini
parents:
diff changeset
75 sindelPrior = 0.000001
3c10d88b55ad improved user interface
mini
parents:
diff changeset
76
3c10d88b55ad improved user interface
mini
parents:
diff changeset
77 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
78 ; probability of an snv or indel noise allele
3c10d88b55ad improved user interface
mini
parents:
diff changeset
79 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
80 ; NB: in the calling model a noise allele is shared in tumor and
3c10d88b55ad improved user interface
mini
parents:
diff changeset
81 ; normal samples, but occurs at any frequency.
3c10d88b55ad improved user interface
mini
parents:
diff changeset
82 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
83 ssnvNoise = 0.0000005
3c10d88b55ad improved user interface
mini
parents:
diff changeset
84 sindelNoise = 0.000001
3c10d88b55ad improved user interface
mini
parents:
diff changeset
85
3c10d88b55ad improved user interface
mini
parents:
diff changeset
86 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
87 ; Fraction of snv noise attributed to strand-bias.
3c10d88b55ad improved user interface
mini
parents:
diff changeset
88 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
89 ; It is not recommended to change this setting. However, if it is
3c10d88b55ad improved user interface
mini
parents:
diff changeset
90 ; essential to turn the strand bias penalization off, the following is
3c10d88b55ad improved user interface
mini
parents:
diff changeset
91 ; recommended:
3c10d88b55ad improved user interface
mini
parents:
diff changeset
92 ; Assuming the current value of ssnvNoiseStrandBiasFrac is 0.5,
3c10d88b55ad improved user interface
mini
parents:
diff changeset
93 ; (1) set ssnvNoiseStrandBiasFrac = 0
3c10d88b55ad improved user interface
mini
parents:
diff changeset
94 ; (2) divide the current ssnvNoise value by 2
3c10d88b55ad improved user interface
mini
parents:
diff changeset
95 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
96 ssnvNoiseStrandBiasFrac = 0.5
3c10d88b55ad improved user interface
mini
parents:
diff changeset
97
3c10d88b55ad improved user interface
mini
parents:
diff changeset
98 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
99 ; minimum MAPQ score for PE reads at tier1:
3c10d88b55ad improved user interface
mini
parents:
diff changeset
100 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
101 minTier1Mapq = 20
3c10d88b55ad improved user interface
mini
parents:
diff changeset
102
3c10d88b55ad improved user interface
mini
parents:
diff changeset
103 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
104 ; minimum MAPQ score for PE and SE reads at tier2:
3c10d88b55ad improved user interface
mini
parents:
diff changeset
105 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
106 minTier2Mapq = 0
3c10d88b55ad improved user interface
mini
parents:
diff changeset
107
3c10d88b55ad improved user interface
mini
parents:
diff changeset
108 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
109 ; Somatic quality score (QSS_NT, NT=ref) below which somatic SNVs are
3c10d88b55ad improved user interface
mini
parents:
diff changeset
110 ; marked as filtered:
3c10d88b55ad improved user interface
mini
parents:
diff changeset
111 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
112 ssnvQuality_LowerBound = 15
3c10d88b55ad improved user interface
mini
parents:
diff changeset
113
3c10d88b55ad improved user interface
mini
parents:
diff changeset
114 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
115 ; Somatic quality score (QSI_NT, NT=ref) below which somatic indels
3c10d88b55ad improved user interface
mini
parents:
diff changeset
116 ; are marked as filtered:
3c10d88b55ad improved user interface
mini
parents:
diff changeset
117 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
118 sindelQuality_LowerBound = 30
3c10d88b55ad improved user interface
mini
parents:
diff changeset
119
3c10d88b55ad improved user interface
mini
parents:
diff changeset
120 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
121 ; Optionally write out read alignments which were altered during the
3c10d88b55ad improved user interface
mini
parents:
diff changeset
122 ; realignment step. At the completion of the workflow run, the
3c10d88b55ad improved user interface
mini
parents:
diff changeset
123 ; realigned reads can be found in:
3c10d88b55ad improved user interface
mini
parents:
diff changeset
124 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
125 ; ${ANALYSIS_DIR}/realigned/{normal,tumor}.realigned.bam
3c10d88b55ad improved user interface
mini
parents:
diff changeset
126 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
127 isWriteRealignedBam = 0
3c10d88b55ad improved user interface
mini
parents:
diff changeset
128
3c10d88b55ad improved user interface
mini
parents:
diff changeset
129 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
130 ; Jobs are parallelized over segments of the reference genome no larger
3c10d88b55ad improved user interface
mini
parents:
diff changeset
131 ; than this size:
3c10d88b55ad improved user interface
mini
parents:
diff changeset
132 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
133 binSize = 25000000
3c10d88b55ad improved user interface
mini
parents:
diff changeset
134
3c10d88b55ad improved user interface
mini
parents:
diff changeset
135 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
136 ; Additional arguments passed to strelka.
3c10d88b55ad improved user interface
mini
parents:
diff changeset
137 ;
3c10d88b55ad improved user interface
mini
parents:
diff changeset
138 extraStrelkaArguments = --remap-input-softclip
3c10d88b55ad improved user interface
mini
parents:
diff changeset
139