Skip to content
This repository was archived by the owner on May 22, 2023. It is now read-only.

Commit a36763d

Browse files
committed
generalise human-specific (simplistic) transgene ENSEMBL/ERCC test
1 parent 186b013 commit a36763d

File tree

1 file changed

+1
-1
lines changed

1 file changed

+1
-1
lines changed

buildindex.nf

+1-1
Original file line numberDiff line numberDiff line change
@@ -151,7 +151,7 @@ process makeTransGeneMatrix {
151151
'''
152152
perl -ne 'if (/^>(\\w+)(?:\\.\\d+)\\s+.*?gene:(\\w+)/){print "$1\\t$2\\n"}elsif(/^>(ERCC\\S+)/){print"$1\\t$1-gene\\n"}' \\
153153
!{fasta} > trans_gene.txt
154-
if (( $(grep -c ENST trans_gene.txt) < 1000 )); then
154+
if (( $(grep -c "^ENS" trans_gene.txt) < 1000 )); then
155155
echo 'Not enough Ensembl transcripts. This test is present because this script section is ugly.'
156156
echo 'Currently it makes an effort to recognise ERCC information.'
157157
echo 'If you want to run a gencode genome, update this section, make this file aware of gencode/Ensembl distinction'

0 commit comments

Comments
 (0)