Skip to content
Snippets Groups Projects
Commit cdd00894 authored by Peter van 't Hof's avatar Peter van 't Hof
Browse files

Added dbsnp

parent dd862229
No related branches found
No related tags found
1 merge request!24Changes after testing
...@@ -193,8 +193,8 @@ task GenotypeGVCFs { ...@@ -193,8 +193,8 @@ task GenotypeGVCFs {
File refFastaIndex File refFastaIndex
File refDict File refDict
File dbsnpVCF File? dbsnpVCF
File dbsnpVCFindex File? dbsnpVCFindex
Int? compressionLevel Int? compressionLevel
Float? memory Float? memory
...@@ -215,7 +215,7 @@ task GenotypeGVCFs { ...@@ -215,7 +215,7 @@ task GenotypeGVCFs {
GenotypeGVCFs \ GenotypeGVCFs \
-R ${refFasta} \ -R ${refFasta} \
-O ${outputPath} \ -O ${outputPath} \
-D ${dbsnpVCF} \ ${"-D " + dbsnpVCF} \
-G StandardAnnotation \ -G StandardAnnotation \
--only-output-calls-starting-in-intervals \ --only-output-calls-starting-in-intervals \
-new-qual \ -new-qual \
...@@ -247,6 +247,9 @@ task HaplotypeCallerGvcf { ...@@ -247,6 +247,9 @@ task HaplotypeCallerGvcf {
Int? compressionLevel Int? compressionLevel
String? gatkJar String? gatkJar
File? dbsnpVCF
File? dbsnpVCFindex
Float? memory Float? memory
Float? memoryMultiplier Float? memoryMultiplier
Int mem = ceil(select_first([memory, 4.0])) Int mem = ceil(select_first([memory, 4.0]))
...@@ -265,6 +268,7 @@ task HaplotypeCallerGvcf { ...@@ -265,6 +268,7 @@ task HaplotypeCallerGvcf {
-O ${gvcfPath} \ -O ${gvcfPath} \
-I ${sep=" -I " inputBams} \ -I ${sep=" -I " inputBams} \
-L ${sep=' -L ' intervalList} \ -L ${sep=' -L ' intervalList} \
${"-D " + dbsnpVCF} \
-contamination ${default=0 contamination} \ -contamination ${default=0 contamination} \
-ERC GVCF -ERC GVCF
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment