File tree Expand file tree Collapse file tree 2 files changed +12
-3
lines changed
subworkflows/local/prepare_genome Expand file tree Collapse file tree 2 files changed +12
-3
lines changed Original file line number Diff line number Diff line change @@ -35,6 +35,7 @@ include { GTF2BED } from '../../../modules/local/gt
3535include { GTF_FILTER } from ' ../../../modules/local/gtf_filter'
3636include { STAR_GENOMEGENERATE_IGENOMES } from ' ../../../modules/local/star_genomegenerate_igenomes'
3737
38+
3839workflow PREPARE_GENOME {
3940 take :
4041 fasta // file: /path/to/genome.fasta
Original file line number Diff line number Diff line change @@ -278,18 +278,26 @@ workflow LNCPIPE {
278278/*
279279* Step 5: Transcript assembly using Stringtie and merge gtf into one
280280*/
281- ch_genome_bam. view()
282281 STRINGTIE_WORKFLOW (
283282 ch_genome_bam,
284283 ch_gtf
285284 )
286285 ch_versions = ch_versions. mix(STRINGTIE_WORKFLOW . out. versions)
287286 ch_merged_gtf = STRINGTIE_WORKFLOW . out. stringtie_gtf_merged
288- ch_merged_gtf. view()
289287/*
290288* Step 6: Compare assembled gtf with known annotations
291289*/
292- GFFCOMPARE (ch_merged_gtf, ch_fasta, ch_gtf)
290+ ch_fasta_meta_fai = ch_fasta. combine(ch_fai)
291+ .map { fasta, fai ->
292+ def meta2 = [ id : fasta. getBaseName(), description : ' Genome FASTA with index' ]
293+ return [ meta2, fasta, fai ]
294+ }
295+
296+ ch_reference_gtf = ch_gtf. map { gtf_file ->
297+ def meta3 = [ id : gtf_file. getBaseName(), description : ' Reference GTF file' ]
298+ return [ meta3, gtf_file ]
299+ }
300+ GFFCOMPARE (ch_merged_gtf, ch_fasta_meta_fai, ch_reference_gtf)
293301 ch_versions = ch_versions. mix(GFFCOMPARE . out. versions)
294302
295303 ch_multiqc_files = ch_multiqc_files
You can’t perform that action at this time.
0 commit comments