Mercurial > repos > mini > strelka
changeset 19:bddc6981a6a7
change version
author | mini |
---|---|
date | Wed, 15 Oct 2014 15:30:28 +0200 |
parents | 3c10d88b55ad |
children | d3ca3d3e5dc7 |
files | strelka.xml strelka_config.sample |
diffstat | 2 files changed, 1 insertions(+), 140 deletions(-) [+] |
line wrap: on
line diff
--- a/strelka.xml Wed Oct 15 14:43:12 2014 +0200 +++ b/strelka.xml Wed Oct 15 15:30:28 2014 +0200 @@ -1,4 +1,4 @@ -<tool id="strelka" name="Strelka" version="1.0.0"> +<tool id="strelka" name="Strelka" version="1.0.1"> <!-- Made by Gregoire Seguin-Henry for geviteam in 2014 --> <description>Strelka</description> <requirements>
--- a/strelka_config.sample Wed Oct 15 14:43:12 2014 +0200 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,139 +0,0 @@ - -; -; User configuration options for Strelka somatic small-variant caller -; workflow: -; - -[user] - -; -; isSkipDepthFilters should be set to 1 to skip depth filtration for -; whole exome or other targeted sequencing data -; -isSkipDepthFilters = 1 - -; -; strelka will not accept input reads above this depth (they will be skipped -; until the depth drops below this value). Set this value <= 0 to disable -; this feature. Using this filter will bound memory usage given extremely high -; depth input, but may be problematic in high-depth targeted sequencing -; applications. -; -maxInputDepth = 10000 - -; -; If the depth filter is not skipped, all variants which occur at a -; depth greater than depthFilterMultiple*chromosome mean depth will be -; filtered out. -; -depthFilterMultiple = 3.0 - -; -; Somatic SNV calls are filtered at sites where greater than this -; fraction of basecalls have been removed by the mismatch density -; filter in either sample. -; -snvMaxFilteredBasecallFrac = 0.4 - -; -; Somatic SNV calls are filtered at sites where greater than this -; fraction of overlapping reads contain deletions which span the SNV -; call site. -; -snvMaxSpanningDeletionFrac = 0.75 - -; -; Somatic indel calls are filtered if they represent an expansion or -; contraction of a repeated pattern with a repeat count greater than -; indelMaxRefRepeat in the reference (ie. if indelMaxRefRepeat is 8, -; then the indel is filtered when it is an expansion/contraction of a -; homopolymer longer than 8 bases, a dinucleotide repeat longer than -; 16 bases, etc.) -; -indelMaxRefRepeat = 8 - -; -; Somatic indel calls are filtered if greater than this fraction of -; basecalls in a window extending 50 bases to each side of an indel's -; call position have been removed by the mismatch density filter. -; -indelMaxWindowFilteredBasecallFrac = 0.3 - -; -; Somatic indels are filtered if they overlap ’interrupted -; homopolymers’ greater than this length. The term 'interrupted -; homopolymer' is used to indicate the longest homopolymer which can -; be found intersecting or adjacent to the called indel when a single -; non-homopolymer base is allowed. -; -indelMaxIntHpolLength = 14 - -; -; prior probability of a somatic snv or indel -; -ssnvPrior = 0.000001 -sindelPrior = 0.000001 - -; -; probability of an snv or indel noise allele -; -; NB: in the calling model a noise allele is shared in tumor and -; normal samples, but occurs at any frequency. -; -ssnvNoise = 0.0000005 -sindelNoise = 0.000001 - -; -; Fraction of snv noise attributed to strand-bias. -; -; It is not recommended to change this setting. However, if it is -; essential to turn the strand bias penalization off, the following is -; recommended: -; Assuming the current value of ssnvNoiseStrandBiasFrac is 0.5, -; (1) set ssnvNoiseStrandBiasFrac = 0 -; (2) divide the current ssnvNoise value by 2 -; -ssnvNoiseStrandBiasFrac = 0.5 - -; -; minimum MAPQ score for PE reads at tier1: -; -minTier1Mapq = 20 - -; -; minimum MAPQ score for PE and SE reads at tier2: -; -minTier2Mapq = 5 - -; -; Somatic quality score (QSS_NT, NT=ref) below which somatic SNVs are -; marked as filtered: -; -ssnvQuality_LowerBound = 15 - -; -; Somatic quality score (QSI_NT, NT=ref) below which somatic indels -; are marked as filtered: -; -sindelQuality_LowerBound = 30 - -; -; Optionally write out read alignments which were altered during the -; realignment step. At the completion of the workflow run, the -; realigned reads can be found in: -; -; ${ANALYSIS_DIR}/realigned/{normal,tumor}.realigned.bam -; -isWriteRealignedBam = 0 - -; -; Jobs are parallelized over segments of the reference genome no larger -; than this size: -; -binSize = 25000000 - -; -; Additional arguments passed to strelka. -; -extraStrelkaArguments = -