Skip to content
Snippets Groups Projects
Verified Commit e0b5c50b authored by Laurent Modolo's avatar Laurent Modolo
Browse files

gatk4: change output to avoid using file_id var as it can be a dict

parent 324f5505
Branches
Tags
No related merge requests found
...@@ -18,7 +18,7 @@ process variant_calling { ...@@ -18,7 +18,7 @@ process variant_calling {
gatk --java-options "-Xmx${xmx_memory}G" HaplotypeCaller \ gatk --java-options "-Xmx${xmx_memory}G" HaplotypeCaller \
-R ${fasta} \ -R ${fasta} \
-I ${bam} \ -I ${bam} \
-O ${file_id}.vcf -O ${bam.simpleName}.vcf
""" """
} }
...@@ -39,7 +39,7 @@ gatk --java-options "-Xmx${xmx_memory}G" SelectVariants \ ...@@ -39,7 +39,7 @@ gatk --java-options "-Xmx${xmx_memory}G" SelectVariants \
-R ${fasta} \ -R ${fasta} \
-V ${vcf} \ -V ${vcf} \
-select-type SNP \ -select-type SNP \
-O ${file_id}_snp.vcf -O ${vcf.simpleName}_snp.vcf
""" """
} }
...@@ -60,7 +60,7 @@ gatk --java-options "-Xmx${xmx_memory}G" SelectVariants \ ...@@ -60,7 +60,7 @@ gatk --java-options "-Xmx${xmx_memory}G" SelectVariants \
-R ${fasta} \ -R ${fasta} \
-V ${vcf} \ -V ${vcf} \
-select-type INDEL \ -select-type INDEL \
-O ${file_id}_indel.vcf -O ${vcf.simpleName}_indel.vcf
""" """
} }
...@@ -84,7 +84,7 @@ gatk --java-options "-Xmx${xmx_memory}G" VariantFiltration \ ...@@ -84,7 +84,7 @@ gatk --java-options "-Xmx${xmx_memory}G" VariantFiltration \
-V ${vcf} \ -V ${vcf} \
--filter-expression "${high_confidence_snp_filter}" \ --filter-expression "${high_confidence_snp_filter}" \
--filter-name "basic_snp_filter" \ --filter-name "basic_snp_filter" \
-O ${file_id}_filtered_snp.vcf -O ${vcf.simpleName}_filtered_snp.vcf
""" """
} }
...@@ -108,7 +108,7 @@ gatk --java-options "-Xmx${xmx_memory}G" VariantFiltration \ ...@@ -108,7 +108,7 @@ gatk --java-options "-Xmx${xmx_memory}G" VariantFiltration \
-V ${vcf} \ -V ${vcf} \
--filter-expression "${high_confidence_indel_filter}" \ --filter-expression "${high_confidence_indel_filter}" \
--filter-name "basic_indel_filter" \ --filter-name "basic_indel_filter" \
-O ${file_id}_filtered_indel.vcf -O ${vcf.simpleName}_filtered_indel.vcf
""" """
} }
...@@ -155,7 +155,7 @@ gatk --java-options "-Xmx${xmx_memory}G" ApplyBQSR \ ...@@ -155,7 +155,7 @@ gatk --java-options "-Xmx${xmx_memory}G" ApplyBQSR \
-R ${fasta} \ -R ${fasta} \
-I ${bam} \ -I ${bam} \
--bqsr-recal-file recal_data_table \ --bqsr-recal-file recal_data_table \
-O ${file_id}_recal.bam -O ${bam.simpleName}_recal.bam
""" """
} }
...@@ -176,7 +176,7 @@ gatk --java-options "-Xmx${xmx_memory}G" HaplotypeCaller \ ...@@ -176,7 +176,7 @@ gatk --java-options "-Xmx${xmx_memory}G" HaplotypeCaller \
-R ${fasta} \ -R ${fasta} \
-I ${bam} \ -I ${bam} \
-ERC GVCF \ -ERC GVCF \
-O ${file_id}.gvcf -O ${bam.simpleName}.gvcf
""" """
} }
...@@ -196,7 +196,7 @@ process gvcf_genotyping { ...@@ -196,7 +196,7 @@ process gvcf_genotyping {
gatk --java-options "-Xmx${xmx_memory}G" GenotypeGVCFs \ gatk --java-options "-Xmx${xmx_memory}G" GenotypeGVCFs \
-R ${fasta} \ -R ${fasta} \
-V ${gvcf} \ -V ${gvcf} \
-O ${file_id}_joint.vcf.gz -O ${bam.simpleName}_joint.vcf.gz
""" """
} }
...@@ -217,7 +217,7 @@ gatk --java-options "-Xmx${xmx_memory}GG" SelectVariants \ ...@@ -217,7 +217,7 @@ gatk --java-options "-Xmx${xmx_memory}GG" SelectVariants \
-R ${fasta} \ -R ${fasta} \
-V ${vcf} \ -V ${vcf} \
-select-type SNP \ -select-type SNP \
-O ${file_id}_joint_snp.vcf -O ${vcf.simpleName}_joint_snp.vcf
""" """
} }
...@@ -259,7 +259,7 @@ process personalized_genome { ...@@ -259,7 +259,7 @@ process personalized_genome {
gatk --java-options "-Xmx${xmx_memory}G" FastaAlternateReferenceMaker\ gatk --java-options "-Xmx${xmx_memory}G" FastaAlternateReferenceMaker\
-R ${reference} \ -R ${reference} \
-V ${vcf} \ -V ${vcf} \
-O ${file_id[0]}_genome.fasta -O ${vcf.simpleName}_genome.fasta
""" """
} }
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment