Skip to content

Commit cb1317e

Browse files
committed
fix GFFREAD
1 parent 76e2382 commit cb1317e

File tree

1 file changed

+16
-7
lines changed

1 file changed

+16
-7
lines changed

subworkflows/local/prepare_genome.nf

Lines changed: 16 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -49,30 +49,39 @@ workflow PREPARE_GENOME {
4949
//
5050
ch_fasta = Channel.empty()
5151
if (fasta.endsWith('.gz')) {
52-
ch_fasta = GUNZIP_FASTA ( [ [:], fasta ] ).gunzip.map{ it[1] }
52+
ch_fasta = GUNZIP_FASTA([[:], fasta]).gunzip.map { it[1] }
5353
ch_versions = ch_versions.mix(GUNZIP_FASTA.out.versions)
54-
} else {
55-
ch_fasta = Channel.value(file(fasta))
54+
}
55+
else {
56+
ch_fasta = Channel.value(file(fasta, checkIfExists: true))
5657
}
5758

5859
//
5960
// Uncompress GTF annotation file or create from GFF3 if required
6061
//
6162
if (gtf) {
6263
if (gtf.endsWith('.gz')) {
63-
ch_gtf = GUNZIP_GTF ( [ [:], gtf ] ).gunzip.map{ it[1] }
64+
ch_gtf = GUNZIP_GTF([[:], gtf]).gunzip.map { it[1] }
6465
ch_versions = ch_versions.mix(GUNZIP_GTF.out.versions)
6566
} else {
6667
ch_gtf = Channel.value(file(gtf))
6768
}
68-
} else if (gff) {
69+
else {
70+
ch_gtf = Channel.value(file(gtf, checkIfExists: true))
71+
}
72+
}
73+
else if (gff) {
6974
if (gff.endsWith('.gz')) {
70-
ch_gff = GUNZIP_GFF ( [ [:], gff ] ).gunzip.map{ it[1] }
75+
ch_gff = GUNZIP_GFF([[:], gff]).gunzip.map { it[1] }
7176
ch_versions = ch_versions.mix(GUNZIP_GFF.out.versions)
7277
} else {
7378
ch_gff = Channel.value(file(gff))
7479
}
75-
ch_gtf = GFFREAD ( [ [:], ch_gff ] ).gtf
80+
else {
81+
ch_gff = Channel.value(file(gff, checkIfExists: true)).map { [[:], it] }
82+
}
83+
84+
ch_gtf = GFFREAD(ch_gff, []).gtf.map { it[1] }
7685
ch_versions = ch_versions.mix(GFFREAD.out.versions)
7786
}
7887

0 commit comments

Comments
 (0)