@@ -1357,18 +1357,15 @@ process samtools_flagstat_after_filter {
13571357* Step 4c: Keep unmapped/remove unmapped reads flagstat
13581358*/
13591359
1360- // merge tuples, for when filtered has been run
1361-
1362- // ch_bam_filtered_flagstat_for_endorspy.view {it -> "DoubleBefore: $it"}
1363-
13641360if (params. run_bam_filtering) {
13651361 ch_flagstat_for_endorspy
13661362 .join(ch_bam_filtered_flagstat_for_endorspy)
13671363 .set{ ch_allflagstats_for_endorspy }
13681364
13691365} else {
1366+ // Add a file entry to match expected no. tuple elements for endorS.py even if not giving second file
13701367 ch_flagstat_for_endorspy
1371- .groupTuple( by : 0 )
1368+ .map { it -> [it[ 0 ], it[ 1 ], file( ' dummy_postfilterflagstat.stats ' )] }
13721369 .set{ ch_allflagstats_for_endorspy }
13731370}
13741371
@@ -1388,9 +1385,16 @@ process endorSpy {
13881385
13891386 script:
13901387 prefix = " ${ name} "
1391- """
1392- endorS.py -o json -n ${ name} ${ stats} ${ poststats}
1393- """
1388+
1389+ if (params. run_bam_filtering) {
1390+ """
1391+ endorS.py -o json -n ${ name} ${ stats} ${ poststats}
1392+ """
1393+ } else {
1394+ """
1395+ endorS.py -o json -n ${ name} ${ stats}
1396+ """
1397+ }
13941398}
13951399
13961400
@@ -1463,7 +1467,7 @@ process markDup{
14631467 prefix = " ${ bam.baseName} "
14641468 size = " ${ params.large_ref} " ? ' -c' : ' '
14651469 """
1466- picard -Xmx${ task.memory.toMega()} M -Xms${ task.memory.toMega()} M MarkDuplicates INPUT=$bam OUTPUT=${ prefix} . _rmdup.bam REMOVE_DUPLICATES=TRUE AS=TRUE METRICS_FILE="${ prefix} .rmdup.metrics" VALIDATION_STRINGENCY=SILENT
1470+ picard -Xmx${ task.memory.toMega()} M -Xms${ task.memory.toMega()} M MarkDuplicates INPUT=$bam OUTPUT=${ prefix} _rmdup.bam REMOVE_DUPLICATES=TRUE AS=TRUE METRICS_FILE="${ prefix} .rmdup.metrics" VALIDATION_STRINGENCY=SILENT
14671471 samtools index "${ size} " ${ prefix} _rmdup.bam
14681472 """
14691473}
@@ -1830,8 +1834,8 @@ if ( params.gatk_ug_jar != '' ) {
18301834 prefix= " ${ bam.baseName} "
18311835 skip_coverage = " ${ params.freebayes_g} " == 0 ? " " : " -g ${ params.freebayes_g} "
18321836 """
1833- freebayes -f ${ fasta} -p ${ params.freebayes_p} -C ${ params.freebayes_C} ${ skip_coverage} ${ bam} > ${ bam.baseName} .vcf
1834- pigz -p ${ task.cpus} ${ bam.baseName} .vcf
1837+ freebayes -f ${ fasta} -p ${ params.freebayes_p} -C ${ params.freebayes_C} ${ skip_coverage} ${ bam} > ${ bam.baseName} .freebayes. vcf
1838+ pigz -p ${ task.cpus} ${ bam.baseName} .freebayes. vcf
18351839 """
18361840 }
18371841
0 commit comments