2021-03-22 12:38:08 +00:00
#!/usr/bin/env nextflow
nextflow.enable.dsl = 2
2021-03-24 09:53:41 +00:00
include { BISMARK_GENOMEPREPARATION } from '../../../../software/bismark/genomepreparation/main.nf' addParams( options: [:] )
include { BISMARK_ALIGN as BISMARK_ALIGN_SE } from '../../../../software/bismark/align/main.nf' addParams( options: [ publish_dir:'test_single_end' ] )
include { BISMARK_ALIGN as BISMARK_ALIGN_PE } from '../../../../software/bismark/align/main.nf' addParams( options: [ publish_dir:'test_paired_end' ] )
2021-03-22 12:38:08 +00:00
2021-05-12 13:56:46 +00:00
//
// Test with single-end data
//
2021-03-22 12:38:08 +00:00
workflow test_bismark_align_single_end {
input = [ [ id:'test', single_end:true ], // meta map
2021-03-25 09:12:08 +00:00
[ file(params.test_data['sarscov2']['illumina']['test_methylated_1_fastq_gz'], checkIfExists: true) ]
2021-03-24 09:53:41 +00:00
]
2021-03-25 09:12:08 +00:00
fasta = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
2021-03-22 12:38:08 +00:00
2021-03-24 09:53:41 +00:00
BISMARK_GENOMEPREPARATION ( fasta )
BISMARK_ALIGN_SE ( input, BISMARK_GENOMEPREPARATION.out.index )
2021-03-22 12:38:08 +00:00
}
2021-05-12 13:56:46 +00:00
//
// Test with paired-end data
//
2021-03-22 12:38:08 +00:00
workflow test_bismark_align_paired_end {
input = [ [ id:'test', single_end:false ], // meta map
2021-03-25 09:12:08 +00:00
[ file(params.test_data['sarscov2']['illumina']['test_methylated_1_fastq_gz'], checkIfExists: true),
file(params.test_data['sarscov2']['illumina']['test_methylated_2_fastq_gz'], checkIfExists: true) ]
2021-03-24 09:53:41 +00:00
]
2021-03-25 09:12:08 +00:00
fasta = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
2021-03-22 12:38:08 +00:00
2021-03-24 09:53:41 +00:00
BISMARK_GENOMEPREPARATION ( fasta )
BISMARK_ALIGN_PE ( input, BISMARK_GENOMEPREPARATION.out.index )
2021-03-22 12:38:08 +00:00
}