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

alntools: double escape awk script

parent f2f617c5
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);
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
......
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