diff --git a/src/nf_modules/gatk4/main.nf b/src/nf_modules/gatk4/main.nf index bf7c0f248c78c82d72e3c14bbcab695d995f2808..ca8999ff8d09a679836c862b1697a6d1c09dd39a 100644 --- a/src/nf_modules/gatk4/main.nf +++ b/src/nf_modules/gatk4/main.nf @@ -98,7 +98,7 @@ gatk --java-options "-Xmx${xmx_memory}G" SelectVariants \ } params.high_confidence_snp_filter = "(QD < 2.0) || (FS > 60.0) || (MQ < 40.0) || (MQRankSum < -12.5) || (ReadPosRankSum < -8.0) || (SOR > 4.0)" -params.high_confidence_snp = "" +params.high_confidence_snp = "--filter-expression \"${params.high_confidence_snp_filter}\" --filter-name \"basic_snp_filter\"" params.high_confidence_snp_out = "" process high_confidence_snp { container = "${container_url}" @@ -124,16 +124,13 @@ process high_confidence_snp { gatk --java-options "-Xmx${xmx_memory}G" VariantFiltration \ -R ${fasta} \ -V ${vcf} \ - ${params.high_confidence_snp_filter} \ - --filter-expression "${high_confidence_snp_filter}" \ - --filter-name "basic_snp_filter" \ + ${params.high_confidence_snp} \ -O ${vcf.simpleName}_filtered_snp.vcf """ } -high_confidence_indel_filter = "QD < 3.0 || FS > 200.0 || ReadPosRankSum < -20.0 || SOR > 10.0" - -params.high_confidence_indels = "" +params.high_confidence_indel_filter = "QD < 3.0 || FS > 200.0 || ReadPosRankSum < -20.0 || SOR > 10.0" +params.high_confidence_indels = "--filter-expression \"${params.high_confidence_indel_filter}\" --filter-name \"basic_indel_filter\"" params.high_confidence_indels_out = "" process high_confidence_indels { container = "${container_url}" @@ -160,8 +157,6 @@ gatk --java-options "-Xmx${xmx_memory}G" VariantFiltration \ -R ${fasta} \ -V ${vcf} \ ${params.high_confidence_indels} \ - --filter-expression "${high_confidence_indel_filter}" \ - --filter-name "basic_indel_filter" \ -O ${vcf.simpleName}_filtered_indel.vcf """ } @@ -394,3 +389,5 @@ gatk --java-options "-Xmx${xmx_memory}G" FastaAlternateReferenceMaker\ """ } + +