Skip to content
GitLab
Explore
Sign in
Primary navigation
Search or go to…
Project
T
tasks
Manage
Activity
Members
Labels
Plan
Issues
Issue boards
Milestones
Wiki
Code
Merge requests
Repository
Branches
Commits
Tags
Repository graph
Compare revisions
Snippets
Build
Pipelines
Jobs
Pipeline schedules
Artifacts
Deploy
Releases
Container Registry
Model registry
Operate
Environments
Monitor
Incidents
Analyze
Value stream analytics
Contributor analytics
CI/CD analytics
Repository analytics
Model experiments
Help
Help
Support
GitLab documentation
Compare GitLab plans
Community forum
Contribute to GitLab
Provide feedback
Keyboard shortcuts
?
Snippets
Groups
Projects
Show more breadcrumbs
biowdl
tasks
Commits
8b27e8aa
Commit
8b27e8aa
authored
6 years ago
by
Cats
Browse files
Options
Downloads
Patches
Plain Diff
styling and paths
parent
fb1a2c68
No related branches found
Branches containing commit
No related tags found
Tags containing commit
1 merge request
!17
Caching related changes
Changes
2
Hide whitespace changes
Inline
Side-by-side
Showing
2 changed files
gatk.wdl
+66
-66
66 additions, 66 deletions
gatk.wdl
picard.wdl
+11
-11
11 additions, 11 deletions
picard.wdl
with
77 additions
and
77 deletions
gatk.wdl
+
66
−
66
View file @
8b27e8aa
...
...
@@ -114,16 +114,16 @@ task GatherBqsrReports {
# Call variants on a single sample with HaplotypeCaller to produce a GVCF
task HaplotypeCallerGvcf {
String? preCommand
Array[File]+ input
_b
ams
Array[File]+ input
_b
ams
_i
ndex
Array[File]+ interval
_l
ist
String gvcf
_basename
File ref
_d
ict
File ref
_f
asta
File ref
_f
asta
_i
ndex
Array[File]+ input
B
ams
Array[File]+ input
B
ams
I
ndex
Array[File]+ interval
L
ist
String gvcf
Path
File ref
D
ict
File ref
F
asta
File ref
F
asta
I
ndex
Float? contamination
Int? compression
_l
evel
String gatk
_j
ar
Int? compression
L
evel
String gatk
J
ar
Float? memory
Float? memoryMultiplier
...
...
@@ -132,20 +132,20 @@ task HaplotypeCallerGvcf {
command {
set -e -o pipefail
${preCommand}
java ${"-Dsamjdk.compression_level=" + compression
_l
evel} \
-Xmx${mem}G -jar ${gatk
_j
ar} \
java ${"-Dsamjdk.compression_level=" + compression
L
evel} \
-Xmx${mem}G -jar ${gatk
J
ar} \
HaplotypeCaller \
-R ${ref
_f
asta} \
-O ${gvcf
_basename}.vcf.gz
\
-I ${sep=" -I " input
_b
ams} \
-L ${sep=' -L ' interval
_l
ist} \
-R ${ref
F
asta} \
-O ${gvcf
Path}
\
-I ${sep=" -I " input
B
ams} \
-L ${sep=' -L ' interval
L
ist} \
-contamination ${default=0 contamination} \
-ERC GVCF
}
output {
File output
_gvcf = "${gvcf_basename}.vcf.gz"
File output
_gvcf_
index =
"${
gvcf
_basename}.vcf.gz
.tbi"
File output
GVCF = gvcfPath
File output
GVCF
index = gvcf
Path + "
.tbi"
}
runtime {
...
...
@@ -155,22 +155,22 @@ task HaplotypeCallerGvcf {
task GenotypeGVCFs {
String? preCommand
File gvcf
_f
iles
File gvcf
_f
ile
_i
ndexes
File gvcf
F
iles
File gvcf
F
ile
I
ndexes
Array[File]+ intervals
String output
_basename
String output
Path
String gatk
_j
ar
String gatk
J
ar
File ref
_f
asta
File ref
_f
asta
_i
ndex
File ref
_d
ict
File ref
F
asta
File ref
F
asta
I
ndex
File ref
D
ict
File dbsnp
_vcf
File dbsnp
_vcf_
index
File dbsnp
VCF
File dbsnp
VCF
index
Int? compression
_l
evel
Int? compression
L
evel
Float? memory
Float? memoryMultiplier
...
...
@@ -179,22 +179,22 @@ task GenotypeGVCFs {
set -e -o pipefail
${preCommand}
java ${"-Dsamjdk.compression_level=" + compression
_l
evel} \
-Xmx${mem}G -jar ${gatk
_j
ar} \
java ${"-Dsamjdk.compression_level=" + compression
L
evel} \
-Xmx${mem}G -jar ${gatk
J
ar} \
GenotypeGVCFs \
-R ${ref
_f
asta} \
-O ${output
_basename + ".vcf.gz"
} \
-D ${dbsnp
_vcf
} \
-R ${ref
F
asta} \
-O ${output
Path
} \
-D ${dbsnp
VCF
} \
-G StandardAnnotation \
--only-output-calls-starting-in-intervals \
-new-qual \
-V ${gvcf
_f
iles} \
-V ${gvcf
F
iles} \
-L ${sep=' -L ' intervals}
}
output {
File output
_vcf
= output
_basename + ".vcf.gz"
File output
_vcf_
index = output
_basename + ".vcf.gz
.tbi"
File output
VCF
= output
Path
File output
VCF
index = output
Path + "
.tbi"
}
runtime{
...
...
@@ -204,19 +204,19 @@ task GenotypeGVCFs {
task CombineGVCFs {
String? preCommand
Array[File]+ gvcf
_f
iles
Array[File]+ gvcf
_f
ile
_i
ndexes
Array[File]+ gvcf
F
iles
Array[File]+ gvcf
F
ile
I
ndexes
Array[File]+ intervals
String output
_basename
String output
Path
String gatk
_j
ar
String gatk
J
ar
File ref
_f
asta
File ref
_f
asta
_i
ndex
File ref
_d
ict
File ref
F
asta
File ref
F
asta
I
ndex
File ref
D
ict
Int? compression
_l
evel
Int? compression
L
evel
Float? memory
Float? memoryMultiplier
...
...
@@ -225,23 +225,23 @@ task CombineGVCFs {
set -e -o pipefail
${preCommand}
if [ ${length(gvcf
_f
iles)} -gt 1 ]; then
java ${"-Dsamjdk.compression_level=" + compression
_l
evel} \
-Xmx${mem}G -jar ${gatk
_j
ar} \
if [ ${length(gvcf
F
iles)} -gt 1 ]; then
java ${"-Dsamjdk.compression_level=" + compression
L
evel} \
-Xmx${mem}G -jar ${gatk
J
ar} \
CombineGVCFs \
-R ${ref
_f
asta} \
-O ${output
_basename + ".vcf.gz"
} \
-V ${sep=' -V ' gvcf
_f
iles} \
-R ${ref
F
asta} \
-O ${output
Path
} \
-V ${sep=' -V ' gvcf
F
iles} \
-L ${sep=' -L ' intervals}
else
ln -sf ${select_first(gvcf
_f
iles)} ${output
_basename + ".vcf.gz"
}
ln -sf ${select_first(gvcf
_f
ile
s)}.tbi ${output_basename + ".vcf.gz
.tbi
"}
ln -sf ${select_first(gvcf
F
iles)} ${output
Path
}
ln -sf ${select_first(gvcf
F
ile
Indexes)} ${outputPath}
.tbi
fi
}
output {
File output
_gvcf
= output
_basename + ".vcf.gz"
File output
_gvcf_
index = output
_basename + ".vcf.gz
.tbi"
File output
GVCF
= output
Path
File output
GVCF
index = output
Path + "
.tbi"
}
runtime {
...
...
@@ -252,13 +252,13 @@ task CombineGVCFs {
task SplitNCigarReads {
String? preCommand
File input
_b
am
File input
_bam_i
ndex
File ref
_f
asta
File ref
_f
asta
_i
ndex
File ref
_d
ict
String output
_b
am
String gatk
_j
ar
File input
B
am
File input
BamI
ndex
File ref
F
asta
File ref
F
asta
I
ndex
File ref
D
ict
String output
B
am
String gatk
J
ar
Array[File]+ intervals
Float? memory
...
...
@@ -268,17 +268,17 @@ task SplitNCigarReads {
command {
set -e -o pipefail
${preCommand}
java -Xms${mem}G -jar ${gatk
_j
ar} \
java -Xms${mem}G -jar ${gatk
J
ar} \
SplitNCigarReads \
-I ${input
_b
am} \
-R ${ref
_f
asta} \
-O ${output
_b
am} \
-I ${input
B
am} \
-R ${ref
F
asta} \
-O ${output
B
am} \
-L ${sep=' -L ' intervals}
}
output {
File bam = output
_b
am
File bam_index = sub(output
_b
am, "\\.bam$", ".bai")
File bam = output
B
am
File bam_index = sub(output
B
am, "\\.bam$", ".bai")
}
runtime {
...
...
This diff is collapsed.
Click to expand it.
picard.wdl
+
11
−
11
View file @
8b27e8aa
...
...
@@ -120,11 +120,11 @@ task MarkDuplicates {
# Combine multiple VCFs or GVCFs from scattered HaplotypeCaller runs
task MergeVCFs {
String? preCommand
Array[File] input
_vcf
s
Array[File] input
_vcfs_i
ndexes
String output
_vcf_
path
Int? compression
_l
evel
String picard
_j
ar
Array[File] input
VCF
s
Array[File] input
VCFsI
ndexes
String output
VCF
path
Int? compression
L
evel
String picard
J
ar
Float? memory
Float? memoryMultiplier
...
...
@@ -135,16 +135,16 @@ task MergeVCFs {
command {
set -e -o pipefail
${preCommand}
java ${"-Dsamjdk.compression_level=" + compression
_l
evel} \
-Xmx${mem}G -jar ${picard
_j
ar} \
java ${"-Dsamjdk.compression_level=" + compression
L
evel} \
-Xmx${mem}G -jar ${picard
J
ar} \
MergeVcfs \
INPUT=${sep=' INPUT=' input
_vcf
s} \
OUTPUT=${output
_vcf_
path}
INPUT=${sep=' INPUT=' input
VCF
s} \
OUTPUT=${output
VCF
path}
}
output {
File output
_vcf
= output
_vcf_
path
File output
_vcf_
index = output
_vcf_
path + ".tbi"
File output
VCF
= output
VCF
path
File output
VCF
index = output
VCF
path + ".tbi"
}
runtime {
...
...
This diff is collapsed.
Click to expand it.
Preview
0%
Loading
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Save comment
Cancel
Please
register
or
sign in
to comment