Skip to content
Snippets Groups Projects
Verified Commit f2f617c5 authored by Laurent Modolo's avatar Laurent Modolo
Browse files

nf_modules: update alntools awk

parent d2c5476d
No related branches found
No related tags found
No related merge requests found
......@@ -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
......
......@@ -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;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment