Skip to content

Commit

Permalink
interim commit
Browse files Browse the repository at this point in the history
  • Loading branch information
abhi18av committed Oct 22, 2020
1 parent 778acd7 commit a4bc8bb
Show file tree
Hide file tree
Showing 2 changed files with 8 additions and 8 deletions.
12 changes: 6 additions & 6 deletions modules/gatk/haplotype_caller/haplotype_caller.nf
Original file line number Diff line number Diff line change
Expand Up @@ -4,8 +4,8 @@ import java.nio.file.Paths


params.gatk_haplotype_caller_contamination = 0
params.gatk_haplotype_caller_make_gvcf = true
params.gatk_haplotype_caller_make_bamout = true
params.make_gvcf = true
params.make_bamout = true

params.gatk_haplotype_caller_memory = '4'
params.gatk_haplotype_caller_java_opts = ''
Expand All @@ -28,10 +28,10 @@ process GATK_HAPLOTYPE_CALLER {
path("${bam_basename}.bamout.bam") optional true

script:
output_suffix = params.gatk_haplotype_caller_make_gvcf ? ".g.vcf.gz" : ".vcf.gz"
output_suffix = params.make_gvcf ? ".g.vcf.gz" : ".vcf.gz"
bam_basename = input_bam.getBaseName()
output_filename = input_bam.getBaseName() + output_suffix
bamout_arg = params.gatk_haplotype_caller_make_bamout ? "-bamout ${bam_basename}.bamout.bam" : ""
bamout_arg = params.make_bamout ? "-bamout ${bam_basename}.bamout.bam" : ""

"""
set -e
Expand All @@ -45,9 +45,9 @@ process GATK_HAPLOTYPE_CALLER {
-contamination ${params.gatk_haplotype_caller_contamination} \
-G StandardAnnotation \
-G StandardHCAnnotation \
${params.gatk_haplotype_caller_make_gvcf ? "-G AS_StandardAnnotation" : ""} \
${params.make_gvcf ? "-G AS_StandardAnnotation" : ""} \
-GQB 10 -GQB 20 -GQB 30 -GQB 40 -GQB 50 -GQB 60 -GQB 70 -GQB 80 -GQB 90 \
${params.gatk_haplotype_caller_make_gvcf ? "-ERC GVCF" : ""} \
${params.make_gvcf ? "-ERC GVCF" : ""} \
${bamout_arg}
"""
}
Expand Down
4 changes: 2 additions & 2 deletions modules/gatk/merge_gvcfs/merge_gvcfs.nf
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
nextflow.enable.dsl = 2

params.gatk_merge_gvcfs_make_gvcf = true
params.make_gvcf = false
params.gatk_merge_gvcfs_memory = '10'


Expand All @@ -21,7 +21,7 @@ process GATK_MERGE_GVCFS {
input_str = input_vcf
.reduce("") { list_of_vcfs, new_vcf -> list_of_vcfs + " --INPUT ${new_vcf}" }
.join(' ')
output_suffix = params.gatk_merge_gvcfs_make_gvcf ? ".g.vcf.gz" : ".vcf.gz"
output_suffix = params.make_gvcf ? ".g.vcf.gz" : ".vcf.gz"
output_filename = input_vcf.getBaseName() + output_suffix

"""
Expand Down

0 comments on commit a4bc8bb

Please sign in to comment.