add stub to modules used in raredisease pipeline (#1206)

* add stub section

* add stub for bcftools norm

* add stub to more modules

* fix expansionhunter tests

* revert changes -picard

* Update stub to write version no.s to a file

* add picard

* revert picard again

* add stubs to more modules

* fix bwamem2

* add bcftools view

* add stubs

* fix svdb query

* review suggestions
krakentools/kreport2krona
Ramprasad Neethiraj 2 years ago committed by GitHub
parent 240ee4328c
commit 49b18b1639
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23

@ -34,4 +34,15 @@ process BCFTOOLS_NORM {
bcftools: \$(bcftools --version 2>&1 | head -n1 | sed 's/^.*bcftools //; s/ .*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.vcf.gz
cat <<-END_VERSIONS > versions.yml
"${task.process}":
bcftools: \$(bcftools --version 2>&1 | head -n1 | sed 's/^.*bcftools //; s/ .*\$//')
END_VERSIONS
"""
}

@ -41,4 +41,15 @@ process BCFTOOLS_VIEW {
bcftools: \$(bcftools --version 2>&1 | head -n1 | sed 's/^.*bcftools //; s/ .*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.vcf.gz
cat <<-END_VERSIONS > versions.yml
"${task.process}":
bcftools: \$(bcftools --version 2>&1 | head -n1 | sed 's/^.*bcftools //; s/ .*\$//')
END_VERSIONS
"""
}

@ -31,4 +31,19 @@ process BWAMEM2_INDEX {
bwamem2: \$(echo \$(bwa-mem2 version 2>&1) | sed 's/.* //')
END_VERSIONS
"""
stub:
"""
mkdir bwamem2
touch bwamem2/${fasta}.0123
touch bwamem2/${fasta}.ann
touch bwamem2/${fasta}.pac
touch bwamem2/${fasta}.amb
touch bwamem2/${fasta}.bwt.2bit.64
cat <<-END_VERSIONS > versions.yml
"${task.process}":
bwamem2: \$(echo \$(bwa-mem2 version 2>&1) | sed 's/.* //')
END_VERSIONS
"""
}

@ -43,4 +43,15 @@ process BWAMEM2_MEM {
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.bam
cat <<-END_VERSIONS > versions.yml
"${task.process}":
bwamem2: \$(echo \$(bwa-mem2 version 2>&1) | sed 's/.* //')
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
}

@ -47,4 +47,14 @@ process CAT_CAT {
pigz: \$( pigz --version 2>&1 | sed 's/pigz //g' )
END_VERSIONS
"""
stub:
"""
touch $prefix
cat <<-END_VERSIONS > versions.yml
"${task.process}":
pigz: \$( pigz --version 2>&1 | sed 's/pigz //g' )
END_VERSIONS
"""
}

@ -44,4 +44,16 @@ process DEEPVARIANT {
deepvariant: \$(echo \$(/opt/deepvariant/bin/run_deepvariant --version) | sed 's/^.*version //; s/ .*\$//' )
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.vcf.gz
touch ${prefix}.g.vcf.gz
cat <<-END_VERSIONS > versions.yml
"${task.process}":
deepvariant: \$(echo \$(/opt/deepvariant/bin/run_deepvariant --version) | sed 's/^.*version //; s/ .*\$//' )
END_VERSIONS
"""
}

@ -37,4 +37,15 @@ process EXPANSIONHUNTER {
expansionhunter: \$( echo \$(ExpansionHunter --version 2>&1) | sed 's/^.*ExpansionHunter v//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.vcf
cat <<-END_VERSIONS > versions.yml
"${task.process}":
expansionhunter: \$( echo \$(ExpansionHunter --version 2>&1) | sed 's/^.*ExpansionHunter v//')
END_VERSIONS
"""
}

@ -44,4 +44,16 @@ process FASTQC {
END_VERSIONS
"""
}
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.html
touch ${prefix}.zip
cat <<-END_VERSIONS > versions.yml
"${task.process}":
fastqc: \$( fastqc --version | sed -e "s/FastQC v//g" )
END_VERSIONS
"""
}

@ -39,4 +39,15 @@ process GATK4_BEDTOINTERVALLIST {
gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.interval_list
cat <<-END_VERSIONS > versions.yml
"${task.process}":
gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
END_VERSIONS
"""
}

@ -37,4 +37,15 @@ process GATK4_CREATESEQUENCEDICTIONARY {
gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.dict
cat <<-END_VERSIONS > versions.yml
"${task.process}":
gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
END_VERSIONS
"""
}

@ -51,4 +51,22 @@ process GATK4_INTERVALLISTTOOLS {
gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
mkdir -p ${prefix}_split/temp_0001_of_6
mkdir -p ${prefix}_split/temp_0002_of_6
mkdir -p ${prefix}_split/temp_0003_of_6
mkdir -p ${prefix}_split/temp_0004_of_6
touch ${prefix}_split/temp_0001_of_6/1scattered.interval_list
touch ${prefix}_split/temp_0002_of_6/2scattered.interval_list
touch ${prefix}_split/temp_0003_of_6/3scattered.interval_list
touch ${prefix}_split/temp_0004_of_6/4scattered.interval_list
cat <<-END_VERSIONS > versions.yml
"${task.process}":
gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
END_VERSIONS
"""
}

@ -42,4 +42,15 @@ process GLNEXUS {
glnexus: \$( echo \$(glnexus_cli 2>&1) | head -n 1 | sed 's/^.*release v//; s/ .*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.bcf
cat <<-END_VERSIONS > versions.yml
"${task.process}":
glnexus: \$( echo \$(glnexus_cli 2>&1) | head -n 1 | sed 's/^.*release v//; s/ .*\$//')
END_VERSIONS
"""
}

@ -28,4 +28,16 @@ process MULTIQC {
multiqc: \$( multiqc --version | sed -e "s/multiqc, version //g" )
END_VERSIONS
"""
stub:
"""
touch multiqc_data
touch multiqc_plots
touch multiqc_report.html
cat <<-END_VERSIONS > versions.yml
"${task.process}":
multiqc: \$( multiqc --version | sed -e "s/multiqc, version //g" )
END_VERSIONS
"""
}

@ -48,4 +48,15 @@ process PICARD_COLLECTHSMETRICS {
picard: \$(echo \$(picard CollectHsMetrics --version 2>&1) | grep -o 'Version:.*' | cut -f2- -d:)
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}_collecthsmetrics.txt
cat <<-END_VERSIONS > versions.yml
"${task.process}":
picard: \$(echo \$(picard CollectHsMetrics --version 2>&1) | grep -o 'Version:.*' | cut -f2- -d:)
END_VERSIONS
"""
}

@ -42,4 +42,24 @@ process PICARD_COLLECTMULTIPLEMETRICS {
picard: \$(picard CollectMultipleMetrics --version 2>&1 | grep -o 'Version.*' | cut -f2- -d:)
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.CollectMultipleMetrics.alignment_summary_metrics
touch ${prefix}.CollectMultipleMetrics.insert_size_metrics
touch ${prefix}.CollectMultipleMetrics.quality_distribution.pdf
touch ${prefix}.CollectMultipleMetrics.base_distribution_by_cycle_metrics
touch ${prefix}.CollectMultipleMetrics.quality_by_cycle_metrics
touch ${prefix}.CollectMultipleMetrics.read_length_histogram.pdf
touch ${prefix}.CollectMultipleMetrics.base_distribution_by_cycle.pdf
touch ${prefix}.CollectMultipleMetrics.quality_by_cycle.pdf
touch ${prefix}.CollectMultipleMetrics.insert_size_histogram.pdf
touch ${prefix}.CollectMultipleMetrics.quality_distribution_metrics
cat <<-END_VERSIONS > versions.yml
"${task.process}":
picard: \$(echo \$(picard CollectMultipleMetrics --version 2>&1) | grep -o 'Version:.*' | cut -f2- -d:)
END_VERSIONS
"""
}

@ -42,4 +42,17 @@ process PICARD_MARKDUPLICATES {
picard: \$(echo \$(picard MarkDuplicates --version 2>&1) | grep -o 'Version:.*' | cut -f2- -d:)
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.bam
touch ${prefix}.bam.bai
touch ${prefix}.MarkDuplicates.metrics.txt
cat <<-END_VERSIONS > versions.yml
"${task.process}":
picard: \$(echo \$(picard MarkDuplicates --version 2>&1) | grep -o 'Version:.*' | cut -f2- -d:)
END_VERSIONS
"""
}

@ -46,4 +46,17 @@ process PICARD_SORTVCF {
picard: \$(picard SortVcf --version 2>&1 | grep -o 'Version:.*' | cut -f2- -d:)
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}_sorted.vcf.gz
touch ${prefix}.bam.bai
touch ${prefix}.MarkDuplicates.metrics.txt
cat <<-END_VERSIONS > versions.yml
"${task.process}":
picard: \$(picard SortVcf --version 2>&1 | grep -o 'Version:.*' | cut -f2- -d:)
END_VERSIONS
"""
}

@ -52,4 +52,72 @@ process QUALIMAP_BAMQC {
qualimap: \$(echo \$(qualimap 2>&1) | sed 's/^.*QualiMap v.//; s/Built.*\$//')
END_VERSIONS
"""
stub:
prefix = task.ext.suffix ? "${meta.id}${task.ext.suffix}" : "${meta.id}"
"""
mkdir -p $prefix/css
mkdir $prefix/images_qualimapReport
mkdir $prefix/raw_data_qualimapReport
cd $prefix/css
touch agogo.css
touch basic.css
touch bgtop.png
touch comment-close.png
touch doctools.js
touch down-pressed.png
touch jquery.js
touch plus.png
touch qualimap_logo_small.png
touch searchtools.js
touch up.png
touch websupport.js
touch ajax-loader.gif
touch bgfooter.png
touch comment-bright.png
touch comment.png
touch down.png
touch file.png
touch minus.png
touch pygments.css
touch report.css
touch underscore.js
touch up-pressed.png
cd ../images_qualimapReport/
touch genome_coverage_0to50_histogram.png
touch genome_coverage_quotes.png
touch genome_insert_size_across_reference.png
touch genome_mapping_quality_histogram.png
touch genome_uniq_read_starts_histogram.png
touch genome_coverage_across_reference.png
touch genome_gc_content_per_window.png
touch genome_insert_size_histogram.png
touch genome_reads_clipping_profile.png
touch genome_coverage_histogram.png
touch genome_homopolymer_indels.png
touch genome_mapping_quality_across_reference.png
touch genome_reads_content_per_read_position.png
cd ../raw_data_qualimapReport
touch coverage_across_reference.txt
touch genome_fraction_coverage.txt
touch insert_size_histogram.txt
touch mapped_reads_nucleotide_content.txt
touch coverage_histogram.txt
touch homopolymer_indels.txt
touch mapped_reads_clipping_profile.txt
touch mapping_quality_across_reference.txt
touch duplication_rate_histogram.txt
touch insert_size_across_reference.txt
touch mapped_reads_gc-content_distribution.txt
touch mapping_quality_histogram.txt
cd ../
touch genome_results.txt
touch qualimapReport.html
cd ../
cat <<-END_VERSIONS > versions.yml
"${task.process}":
qualimap: \$(echo \$(qualimap 2>&1) | sed 's/^.*QualiMap v.//; s/Built.*\$//')
END_VERSIONS
"""
}

@ -29,4 +29,14 @@ process SAMTOOLS_FAIDX {
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
stub:
"""
touch ${fasta}.fai
cat <<-END_VERSIONS > versions.yml
"${task.process}":
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
}

@ -33,4 +33,16 @@ process SAMTOOLS_INDEX {
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
stub:
"""
touch ${input}.bai
touch ${input}.crai
touch ${input}.csi
cat <<-END_VERSIONS > versions.yml
"${task.process}":
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
}

@ -38,4 +38,16 @@ process SAMTOOLS_MERGE {
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
stub:
prefix = task.ext.suffix ? "${meta.id}${task.ext.suffix}" : "${meta.id}"
def file_type = input_files[0].getExtension()
"""
touch ${prefix}.${file_type}
cat <<-END_VERSIONS > versions.yml
"${task.process}":
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
}

@ -28,4 +28,15 @@ process SAMTOOLS_SORT {
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.bam
cat <<-END_VERSIONS > versions.yml
"${task.process}":
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
}

@ -34,4 +34,15 @@ process SAMTOOLS_STATS {
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${input}.stats
cat <<-END_VERSIONS > versions.yml
"${task.process}":
samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
END_VERSIONS
"""
}

@ -70,5 +70,4 @@ process SVDB_QUERY {
svdb: \$( echo \$(svdb) | head -1 | sed 's/usage: SVDB-\\([0-9]\\.[0-9]\\.[0-9]\\).*/\\1/' )
END_VERSIONS
"""
}

@ -30,4 +30,16 @@ process TABIX_BGZIPTABIX {
tabix: \$(echo \$(tabix -h 2>&1) | sed 's/^.*Version: //; s/ .*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.gz
touch ${prefix}.gz.tbi
cat <<-END_VERSIONS > versions.yml
"${task.process}":
tabix: \$(echo \$(tabix -h 2>&1) | sed 's/^.*Version: //; s/ .*\$//')
END_VERSIONS
"""
}

@ -27,4 +27,15 @@ process TABIX_TABIX {
tabix: \$(echo \$(tabix -h 2>&1) | sed 's/^.*Version: //; s/ .*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${tab}.tbi
cat <<-END_VERSIONS > versions.yml
"${task.process}":
tabix: \$(echo \$(tabix -h 2>&1) | sed 's/^.*Version: //; s/ .*\$//')
END_VERSIONS
"""
}

@ -36,4 +36,16 @@ process TIDDIT_COV {
tiddit: \$(echo \$(tiddit 2>&1) | sed 's/^.*TIDDIT-//; s/ .*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch $prefix.wig
touch $prefix.tab
cat <<-END_VERSIONS > versions.yml
"${task.process}":
tiddit: \$(echo \$(tiddit 2>&1) | sed 's/^.*TIDDIT-//; s/ .*\$//')
END_VERSIONS
"""
}

@ -38,4 +38,17 @@ process TIDDIT_SV {
tiddit: \$(echo \$(tiddit 2>&1) | sed 's/^.*TIDDIT-//; s/ .*\$//')
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch $prefix.vcf
touch $prefix.ploidy.tab
touch $prefix.signals.tab
cat <<-END_VERSIONS > versions.yml
"${task.process}":
tiddit: \$(echo \$(tiddit 2>&1) | sed 's/^.*TIDDIT-//; s/ .*\$//')
END_VERSIONS
"""
}

@ -35,4 +35,15 @@ process UCSC_WIGTOBIGWIG {
ucsc: $VERSION
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}.bw
cat <<-END_VERSIONS > versions.yml
"${task.process}":
ucsc: $VERSION
END_VERSIONS
"""
}

@ -33,4 +33,15 @@ process UNTAR {
untar: \$(echo \$(tar --version 2>&1) | sed 's/^.*(GNU tar) //; s/ Copyright.*\$//')
END_VERSIONS
"""
stub:
untar = archive.toString() - '.tar.gz'
"""
touch $untar
cat <<-END_VERSIONS > versions.yml
"${task.process}":
untar: \$(echo \$(tar --version 2>&1) | sed 's/^.*(GNU tar) //; s/ Copyright.*\$//')
END_VERSIONS
"""
}

@ -39,4 +39,15 @@ process VCFANNO {
vcfanno: \$(echo \$(vcfanno 2>&1 | grep version | cut -f3 -d' ' ))
END_VERSIONS
"""
stub:
def prefix = task.ext.prefix ?: "${meta.id}"
"""
touch ${prefix}_annotated.vcf
cat <<-END_VERSIONS > versions.yml
"${task.process}":
vcfanno: \$(echo \$(vcfanno 2>&1 | grep version | cut -f3 -d' ' ))
END_VERSIONS
"""
}

@ -7,8 +7,8 @@ include { EXPANSIONHUNTER } from '../../../modules/expansionhunter/main.nf'
workflow test_expansionhunter {
input = [ [ id:'test', gender:'male' ], // meta map
file(params.test_data['homo_sapiens']['illumina']['test_paired_end_recalibrated_sorted_bam'], checkIfExists: true),
file(params.test_data['homo_sapiens']['illumina']['test_paired_end_recalibrated_sorted_bam_bai'], checkIfExists: true),
file(params.test_data['homo_sapiens']['illumina']['test_paired_end_sorted_bam'], checkIfExists: true),
file(params.test_data['homo_sapiens']['illumina']['test_paired_end_sorted_bam_bai'], checkIfExists: true),
]
fasta = file(params.test_data['homo_sapiens']['genome']['genome_fasta'], checkIfExists: true)
variant_catalog = file(params.test_data['homo_sapiens']['genome']['repeat_expansions'], checkIfExists: true)

@ -4,4 +4,4 @@
- expansionhunter
files:
- path: output/expansionhunter/test.vcf
md5sum: ef6c2101d7bd67211bb5a5a132690e02
md5sum: cfd4a1d35c0e469b99eb6aaa6d22de76

Loading…
Cancel
Save