Skip to content

Commit

Permalink
Revert "fix dipcall bug"
Browse files Browse the repository at this point in the history
This reverts commit 36b9768.
  • Loading branch information
fellen31 committed May 17, 2023
1 parent 36b9768 commit 713272e
Show file tree
Hide file tree
Showing 2 changed files with 7 additions and 20 deletions.
5 changes: 0 additions & 5 deletions modules.json
Original file line number Diff line number Diff line change
Expand Up @@ -21,11 +21,6 @@
"installed_by": ["modules"],
"patch": "modules/nf-core/gfastats/gfastats.diff"
},
"glnexus": {
"branch": "master",
"git_sha": "911696ea0b62df80e900ef244d7867d177971f73",
"installed_by": ["modules"]
},
"hifiasm": {
"branch": "master",
"git_sha": "911696ea0b62df80e900ef244d7867d177971f73",
Expand Down
22 changes: 7 additions & 15 deletions modules/local/dipcall.nf
Original file line number Diff line number Diff line change
Expand Up @@ -16,17 +16,9 @@ process DIPCALL {
tuple val(meta), path(haplotype_1), path(haplotype_2), path(reference), path(fai), path(mmi), path(par), val(sex)

output:
tuple val(meta), path("*.dip.vcf.gz") , emit: variant_calls
tuple val(meta), path("*.dip.bed") , emit: confident_regions
tuple val(meta), path("*.paf.gz") , emit: paf
tuple val(meta), path("*.sam.gz") , emit: sam
tuple val(meta), path("*.var.gz") , emit: var
tuple val(meta), path("*.bam") , emit: bam
tuple val(meta), path("*.hap1.bed") , emit: hap1_bed
tuple val(meta), path("*.hap2.bed") , emit: hap2_bed
tuple val(meta), path("*.pair.vcf.gz"), emit: pair
tuple val(meta), path("*.tmp") , emit: tmp, optional: true
path "versions.yml" , emit: versions
tuple val(meta), path("*.dip.vcf.gz"), emit: variant_calls
tuple val(meta), path("*.dip.bed") , emit: confident_regions
path "versions.yml" , emit: versions
//TODO: Add all outputs

when:
Expand All @@ -40,8 +32,8 @@ process DIPCALL {
minimap2 -c --paf-no-hit -xasm5 --cs -t ${task.cpus} ${mmi} ${haplotype_1} 2> ${prefix}.hap1.paf.gz.log | gzip > ${prefix}.hap1.paf.gz
minimap2 -c --paf-no-hit -xasm5 --cs -t ${task.cpus} ${mmi} ${haplotype_2} 2> ${prefix}.hap2.paf.gz.log | gzip > ${prefix}.hap2.paf.gz
minimap2 -a -xasm5 --cs -t ${task.cpus} ${mmi} ${haplotype_1} 2> ${prefix}.hap1.sam.gz.log | gzip > ${prefix}.hap1.sam.gz
minimap2 -a -xasm5 --cs -t ${task.cpus} ${mmi} ${haplotype_2} 2> ${prefix}.hap2.sam.gz.log | gzip > ${prefix}.hap2.sam.gz
minimap2 -a -xasm5 --cs -t ${task.cpus} ${mmi} ${prefix}.asm.bp.hap1.p_ctg.fa 2> ${prefix}.hap1.sam.gz.log | gzip > ${prefix}.hap1.sam.gz
minimap2 -a -xasm5 --cs -t ${task.cpus} ${mmi} ${prefix}.asm.bp.hap2.p_ctg.fa 2> ${prefix}.hap2.sam.gz.log | gzip > ${prefix}.hap2.sam.gz
gzip -dc ${prefix}.hap1.paf.gz | sort -k6,6 -k8,8n | k8 /usr/local/bin/paftools.js call - 2> ${prefix}.hap1.var.gz.vst | gzip > ${prefix}.hap1.var.gz
gzip -dc ${prefix}.hap2.paf.gz | sort -k6,6 -k8,8n | k8 /usr/local/bin/paftools.js call - 2> ${prefix}.hap2.var.gz.vst | gzip > ${prefix}.hap2.var.gz
Expand All @@ -68,8 +60,8 @@ process DIPCALL {
minimap2 -c --paf-no-hit -xasm5 --cs -t ${task.cpus} ${mmi} ${haplotype_1} 2> ${prefix}.hap1.paf.gz.log | gzip > ${prefix}.hap1.paf.gz
minimap2 -c --paf-no-hit -xasm5 --cs -t ${task.cpus} ${mmi} ${haplotype_2} 2> ${prefix}.hap2.paf.gz.log | gzip > ${prefix}.hap2.paf.gz
minimap2 -a -xasm5 --cs -t ${task.cpus} ${mmi} ${haplotype_1} 2> ${prefix}.hap1.sam.gz.log | gzip > ${prefix}.hap1.sam.gz
minimap2 -a -xasm5 --cs -t ${task.cpus} ${mmi} ${haplotype_2} 2> ${prefix}.hap2.sam.gz.log | gzip > ${prefix}.hap2.sam.gz
minimap2 -a -xasm5 --cs -t ${task.cpus} ${mmi} ${prefix}.asm.bp.hap1.p_ctg.fa 2> ${prefix}.hap1.sam.gz.log | gzip > ${prefix}.hap1.sam.gz
minimap2 -a -xasm5 --cs -t ${task.cpus} ${mmi} ${prefix}.asm.bp.hap2.p_ctg.fa 2> ${prefix}.hap2.sam.gz.log | gzip > ${prefix}.hap2.sam.gz
gzip -dc ${prefix}.hap1.paf.gz | sort -k6,6 -k8,8n | k8 /usr/local/bin/paftools.js call - 2> ${prefix}.hap1.var.gz.vst | gzip > ${prefix}.hap1.var.gz
gzip -dc ${prefix}.hap2.paf.gz | sort -k6,6 -k8,8n | k8 /usr/local/bin/paftools.js call - 2> ${prefix}.hap2.var.gz.vst | gzip > ${prefix}.hap2.var.gz
Expand Down

0 comments on commit 713272e

Please sign in to comment.