-
Notifications
You must be signed in to change notification settings - Fork 4
/
gatk-haplotype_caller.cwl
102 lines (96 loc) · 1.98 KB
/
gatk-haplotype_caller.cwl
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
cwlVersion: v1.0
class: CommandLineTool
id: gatk_haplotypecaller
requirements:
- class: InlineJavascriptRequirement
- class: DockerRequirement
dockerPull: broadinstitute/gatk3:3.6-0
- class: InitialWorkDirRequirement
listing:
- entry: $(inputs.reference_genome)
- entry: $(inputs.dict)
- class: ResourceRequirement
outdirMin: 7500
tmpdirMin: 7700
hints:
- class: ResourceRequirement
coresMin: 14
ramMin: 54000
baseCommand:
- java
- -jar
- /usr/GenomeAnalysisTK.jar
- '-T'
- HaplotypeCaller
- --never_trim_vcf_format_field
inputs:
- id: reference_genome
type: File
inputBinding:
position: 1
prefix: '-R'
secondaryFiles:
- .fai
- id: dict
type: File
# from gatk-base_recalibration_print_reads.cwl
- id: input
type: File
inputBinding:
position: 2
prefix: '-I'
secondaryFiles:
- ^.bai
# default is to analyse the complete genome
- id: chromosome
type: string?
inputBinding:
position: 3
prefix: '-L'
- id: ploidy
type: int?
inputBinding:
position: 5
prefix: '-ploidy'
- id: gqb
type:
- "null"
- type: array
items: int
inputBinding: { prefix: "--GVCFGQBands" }
inputBinding:
position: 12
- id: threads
type: string?
default: "2"
arguments:
- position: 0
prefix: '--num_cpu_threads_per_data_thread'
valueFrom: $(inputs.threads)
- position: 0
prefix: '-dt'
valueFrom: 'NONE'
- position: 0
prefix: '-rf'
valueFrom: 'BadCigar'
- position: 0
prefix: '-ERC'
valueFrom: 'GVCF'
- position: 0
prefix: '-variant_index_type'
valueFrom: 'LINEAR'
- position: 0
prefix: '-variant_index_parameter'
valueFrom: '128000'
- position: 0
prefix: '-o'
valueFrom: $(inputs.input.nameroot).vcf.gz
outputs:
- id: gvcf
type: File
outputBinding:
glob: "*.gz"
secondaryFiles:
#- .idx
- .tbi
label: gatk3-haplotypecaller