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

alntools: try fix for bam bai

parent 972e601e
No related branches found
No related tags found
No related merge requests found
...@@ -21,6 +21,8 @@ process bam2ec { ...@@ -21,6 +21,8 @@ process bam2ec {
script: script:
""" """
cp ${bam} file_bam.bam
cp ${bam_idx} file_bam.bam.bai
awk -F"[\\t;]" ' awk -F"[\\t;]" '
\$3=="exon" { \$3=="exon" {
ID=gensub(/transcript_id \\"(.*)\\"/, "\\\\1", "g", \$11); ID=gensub(/transcript_id \\"(.*)\\"/, "\\\\1", "g", \$11);
...@@ -31,6 +33,6 @@ END{ ...@@ -31,6 +33,6 @@ END{
{print i"\\t"LEN[i]} {print i"\\t"LEN[i]}
} }
' ${gtf} > ${gtf.simpleName}_transcripts_lengths.tsv ' ${gtf} > ${gtf.simpleName}_transcripts_lengths.tsv
alntools bam2ec ${params.bam2ec} -t ${gtf.simpleName}_transcripts_lengths.tsv ${bam} ${bam.simpleName}.bin alntools bam2ec ${params.bam2ec} -t ${gtf.simpleName}_transcripts_lengths.tsv file_bam.bam ${bam.simpleName}.bin
""" """
} }
\ No newline at end of file
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