Skip to content

Commit

Permalink
Fix line length
Browse files Browse the repository at this point in the history
  • Loading branch information
rernst committed May 24, 2022
1 parent d9c3db1 commit 3a945bc
Showing 1 changed file with 27 additions and 5 deletions.
32 changes: 27 additions & 5 deletions GATK/4.2.1.0/VariantFiltration.nf
Original file line number Diff line number Diff line change
Expand Up @@ -15,12 +15,34 @@ process VariantFiltrationSnpIndel {
ext_vcf = params.compress || vcf_file.getExtension() == ".gz" ? ".vcf.gz" : ".vcf"
ext_vcf_index = params.compress || vcf_file.getExtension() == ".gz" ? ".tbi" : ".idx"
"""
gatk --java-options "-Xmx${task.memory.toGiga()-4}G" SelectVariants --reference ${params.genome} --variant $vcf_file --output ${vcf_file.simpleName}.snp${ext_vcf} --select-type-to-exclude INDEL
gatk --java-options "-Xmx${task.memory.toGiga()-4}G" SelectVariants --reference ${params.genome} --variant $vcf_file --output ${vcf_file.simpleName}.indel${ext_vcf} --select-type-to-include INDEL
gatk --java-options "-Xmx${task.memory.toGiga()-4}G" SelectVariants \
--reference ${params.genome} \
--variant $vcf_file \
--output ${vcf_file.simpleName}.snp${ext_vcf} \
--select-type-to-exclude INDEL
gatk --java-options "-Xmx${task.memory.toGiga()-4}G" VariantFiltration --reference ${params.genome} --variant ${vcf_file.simpleName}.snp${ext_vcf} --output ${vcf_file.simpleName}.snp_filter${ext_vcf} ${params.snp_filter} ${params.snp_cluster}
gatk --java-options "-Xmx${task.memory.toGiga()-4}G" VariantFiltration --reference ${params.genome} --variant ${vcf_file.simpleName}.indel${ext_vcf} --output ${vcf_file.simpleName}.indel_filter${ext_vcf} ${params.indel_filter}
gatk --java-options "-Xmx${task.memory.toGiga()-4}G" SelectVariants \
--reference ${params.genome} \
--variant $vcf_file \
--output ${vcf_file.simpleName}.indel${ext_vcf} \
--select-type-to-include INDEL
gatk --java-options "-Xmx${task.memory.toGiga()-4}G" MergeVcfs --INPUT ${vcf_file.simpleName}.snp_filter${ext_vcf} --INPUT ${vcf_file.simpleName}.indel_filter${ext_vcf} --OUTPUT ${vcf_file.simpleName}.filter${ext_vcf}
gatk --java-options "-Xmx${task.memory.toGiga()-4}G" VariantFiltration \
--reference ${params.genome} \
--variant ${vcf_file.simpleName}.snp${ext_vcf} \
--output ${vcf_file.simpleName}.snp_filter${ext_vcf} \
${params.snp_filter} \
${params.snp_cluster}
gatk --java-options "-Xmx${task.memory.toGiga()-4}G" VariantFiltration \
--reference ${params.genome} \
--variant ${vcf_file.simpleName}.indel${ext_vcf} \
--output ${vcf_file.simpleName}.indel_filter${ext_vcf} \
${params.indel_filter}
gatk --java-options "-Xmx${task.memory.toGiga()-4}G" MergeVcfs \
--INPUT ${vcf_file.simpleName}.snp_filter${ext_vcf} \
--INPUT ${vcf_file.simpleName}.indel_filter${ext_vcf} \
--OUTPUT ${vcf_file.simpleName}.filter${ext_vcf}
"""
}

0 comments on commit 3a945bc

Please sign in to comment.