diff --git a/workflow/rna-seq.nf b/workflow/rna-seq.nf index e7ae99807e08cae3f3555406885b3987ee683e12..e4427f61a1d8074dcf71bf8d120b00a18754b8b9 100755 --- a/workflow/rna-seq.nf +++ b/workflow/rna-seq.nf @@ -401,9 +401,9 @@ process makeFeatureCounts { script: """ if [ "${endsManual_featureCounts }" == "se" ]; then - featureCounts -R SAM -p -G ${genome} -T `nproc` -s 1 -a ${gtf} -o ${repRID}.featureCounts ${repRID}.sorted.deduped.bam; + featureCounts -R SAM -p -G ${genome} -T `nproc` -s 1 -a ${gtf} -o ${repRID}.featureCounts -g 'gene_name' ${repRID}.sorted.deduped.bam; elif [ "${endsManual_featureCounts }" == "pe" ]; then - featureCounts -R SAM -p -G ${genome} -T `nproc` -s 1 -p -a ${gtf} -o ${repRID}.featureCounts ${repRID}.sorted.deduped.bam; + featureCounts -R SAM -p -G ${genome} -T `nproc` -s 1 -p -a ${gtf} -o ${repRID}.featureCounts -g 'gene_name' ${repRID}.sorted.deduped.bam; fi; Rscript calculateTPM.R --count "${repRID}.featureCounts"; """