File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -721,14 +721,15 @@ process samtools_filter {
721721 prefix= " $bam " - ~/ (\. bam)? /
722722 rm_type = " ${ params.bam_unmapped_keep_type} " == ' bam' ? ' bam' : ' fastq'
723723 rm_unmapped = " ${ params.bam_discard_unmapped} " ? " rm *.unmapped.${ rm_type} " : ' '
724- fq_convert = " ${ params.bam_unmapped_to_fastq} " ? " samtools fastq -tn " ${prefix}. unmapped. bam" | gzip > " ${prefix}. unmapped. fq. gz" " :
724+ fq_convert = " ${ params.bam_unmapped_to_fastq} " ? " samtools fastq -tn ${ prefix} .unmapped.bam | gzip > ${ prefix} .unmapped.fq.gz" : ' '
725725
726726
727727 if (" ${ params.bam_separate_unmapped} " ){
728728 """
729729 samtools view -h $bam | tee >(samtools view - -@ ${ task.cpus} -f4 -q ${ params.bam_mapping_quality_threshold} -o ${ prefix} .unmapped.bam) >(samtools view - -@ ${ task.cpus} -F4 -q ${ params.bam_mapping_quality_threshold} -o ${ prefix} .filtered.bam)
730730 samtools index ${ prefix} .filtered.bam
731731 $fq_convert
732+ $rm_unmapped
732733 """
733734 } else {
734735 """
You can’t perform that action at this time.
0 commit comments