diff --git a/src/nf_modules/alntools/main.nf b/src/nf_modules/alntools/main.nf index 148e976b8c8c07d41b081d7513eb1e51b916efbe..a302ec2dedede714c0dd36a3f39d70590603416c 100644 --- a/src/nf_modules/alntools/main.nf +++ b/src/nf_modules/alntools/main.nf @@ -21,14 +21,14 @@ process bam2ec { script: """ -awk -F"[\\t;]" ' +awk -F"[\t;]" ' \$3=="exon" { ID=gensub(/transcript_id \"(.*)\"/, "\\1", "g", \$11); LEN[ID]+=\$5-\$4+1; } END{ for(i in LEN) - {print i"\\t"LEN[i]} + {print i"\t"LEN[i]} } ' ${gtf} > ${gtf.simpleName}_transcripts_lengths.tsv alntools bam2ec ${params.bam2ec} -t ${gtf.simpleName}_transcripts_lengths.tsv -c 8 ${bam} ${bam.simpleName}.bin diff --git a/src/nf_modules/kb/main.nf b/src/nf_modules/kb/main.nf index d174b19c09304156ef1268568f17f5ed25da5548..e5895ba68f4222fcb15365819bcb68b278fe38d1 100644 --- a/src/nf_modules/kb/main.nf +++ b/src/nf_modules/kb/main.nf @@ -102,7 +102,7 @@ workflow count { switch(params.kb_protocol) { case "marsseq": split(fastq, config) - kb_marseq(index, split.out.fastq.view(), transcript_to_gene, whitelist_optional) + kb_marseq(index.collect(), split.out.fastq.view(), transcript_to_gene, whitelist_optional) kb_marseq.out.counts.set{res_counts} kb_marseq.out.report.set{res_report} break;