@@ -172,10 +172,6 @@ workflow DIRECTRNA{
172172 // INPUT_CHECK
173173 INPUT_CHECK ( ch_input )
174174 .set { ch_sample }
175- // playing around with channel transformations
176- // /ch_sample
177- // /.map { it -> [ it[0], it[1] ] } // take sample, replicate, reads
178- // /.set { ch_fastq }
179175
180176 // QC of fastq files
181177 // / MODULES: NANOQ & SEQUALI
@@ -196,7 +192,7 @@ workflow DIRECTRNA{
196192
197193 // Prepare the reference files
198194 // / SUBWORKFLOW: PREPARE_REFERENCE
199- /* if (!params.skip_prepare_reference) {
195+ if (! params. skip_prepare_reference) {
200196 PREPARE_REFERENCE (
201197 params. genome_fasta,
202198 params. genome_fasta_index,
@@ -223,8 +219,11 @@ workflow DIRECTRNA{
223219 ch_transcriptome_fasta = PREPARE_REFERENCE . out. transcriptome_fasta
224220 ch_annotation_gtf = PREPARE_REFERENCE . out. annotation_gtf
225221 ch_jaffal_reference_dir = PREPARE_REFERENCE . out. jaffal_reference
222+ // Combine genome fasta with genome fasta index into single channel -
223+ // some software expect both files in a single path/channel
224+ ch_genome_fasta_with_index = ch_genome_fasta. combine(ch_genome_fasta_index)
226225
227- // initialize sqanti qc references
226+ /* // initialize sqanti qc references
228227 if (!params.skip_sqanti_qc) {
229228 if (params.sqanti_qc_cage) {
230229 ch_sqanti_qc_cage_bed = PREPARE_REFERENCE.out.sqanti_qc_cage_bed
@@ -242,9 +241,6 @@ workflow DIRECTRNA{
242241*/
243242
244243
245- // Combine genome fasta with genome fasta index into single channel
246- ch_genome_fasta_with_index = ch_genome_fasta. combine(ch_genome_fasta_index)
247-
248244 // Mapping and sorting
249245 // SUBWORKFLOW: MAPPING
250246 //
@@ -347,6 +343,7 @@ workflow DIRECTRNA{
347343 ch_versions = ch_versions. mix(GFFREAD_GETFASTA . out. versions. first())
348344 }
349345
346+ // TODO
350347 // ISOQUANT
351348 // if (!params.skip_isoquant) {
352349 // ISOQUANT( ch_mixed_bam, ch_annotation_gtf, ch_genome_fasta)
@@ -424,13 +421,12 @@ workflow DIRECTRNA{
424421 }
425422
426423 //
427- // Transcript quantification
424+ // Transcript quantification will be added in first update to pipeline
428425 // TransSigner
429426 //if (!params.skip_quantification && !params.skip_mapping)
430427 // TRANSIGNER_MAP
431428 // TRANSIGNER_
432429 // TRANSIGNER_QUANT
433-
434430*/
435431
436432 //
0 commit comments