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

emase-zero: fix gene_to_transcript file to add strain info to gene name

parent 231d2aaf
No related branches found
No related tags found
No related merge requests found
......@@ -45,7 +45,9 @@ process emase {
script:
"""
grep ">" ${fasta} | sed 's/>//' > tr_list.txt
grep -Fw -f tr_list.txt ${gene_to_transcript} > gene_to_transcript.txt
grep -Fw -f tr_list.txt ${gene_to_transcript} | \
sed -E 's/(.+gn[0-9]+)(.+)(_.+)/\\1\\3\\2\\3/g' > \
gene_to_transcript.txt
emase-zero ${params.count} \
-o ${bin.simpleName}.quantified \
-l ${transcript_length} \
......
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