Skip to content

Commit 49956c1

Browse files
committed
Remove -b flag that @mahesh-panchal correctly pointed out was the problem
1 parent 8fafdbc commit 49956c1

1 file changed

Lines changed: 2 additions & 2 deletions

File tree

main.nf

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -2037,12 +2037,12 @@ process pmdtools {
20372037
"""
20382038
## Note all `samtools view` steps ARE required for some reason
20392039
#Run Filtering step
2040-
samtools calmd -b ${bam} ${fasta} | samtools view -h - | pmdtools --threshold ${params.pmdtools_threshold} ${treatment} ${snpcap} --header | samtools view -Sb - > "${libraryid}".pmd.bam
2040+
samtools calmd ${bam} ${fasta} | pmdtools --threshold ${params.pmdtools_threshold} ${treatment} ${snpcap} --header | samtools view -Sb - > "${libraryid}".pmd.bam
20412041
20422042
#Run Calc Range step
20432043
## To allow early shut off of pipe: https://github.com/nextflow-io/nextflow/issues/1564
20442044
trap 'if [[ \$? == 141 ]]; then echo "Shutting samtools early due to -n parameter" && samtools index ${libraryid}.pmd.bam ${size}; exit 0; fi' EXIT
2045-
samtools calmd -b ${bam} ${fasta} | samtools view -h - | pmdtools --deamination ${platypus} --range ${params.pmdtools_range} ${treatment} ${snpcap} -n ${params.pmdtools_max_reads} > "${libraryid}".cpg.range."${params.pmdtools_range}".txt
2045+
samtools calmd ${bam} ${fasta} | pmdtools --deamination ${platypus} --range ${params.pmdtools_range} ${treatment} ${snpcap} -n ${params.pmdtools_max_reads} > "${libraryid}".cpg.range."${params.pmdtools_range}".txt
20462046
20472047
samtools index ${libraryid}.pmd.bam ${size}
20482048
"""

0 commit comments

Comments
 (0)