2021-02-12 14:24:09 +00:00
|
|
|
#!/usr/bin/env nextflow
|
|
|
|
|
|
|
|
nextflow.enable.dsl = 2
|
|
|
|
|
2021-11-26 07:58:40 +00:00
|
|
|
include { BWAMEM2_INDEX } from '../../../../modules/bwamem2/index/main.nf'
|
|
|
|
include { BWAMEM2_MEM } from '../../../../modules/bwamem2/mem/main.nf'
|
2021-02-12 14:24:09 +00:00
|
|
|
|
2021-05-12 13:56:46 +00:00
|
|
|
//
|
|
|
|
// Test with single-end data
|
|
|
|
//
|
2021-02-12 14:24:09 +00:00
|
|
|
workflow test_bwamem2_mem_single_end {
|
2021-11-26 07:58:40 +00:00
|
|
|
input = [
|
|
|
|
[ id:'test', single_end:true ], // meta map
|
|
|
|
[
|
|
|
|
file(params.test_data['sarscov2']['illumina']['test_1_fastq_gz'], checkIfExists: true)
|
|
|
|
]
|
|
|
|
]
|
2021-03-25 09:24:21 +00:00
|
|
|
fasta = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
|
2021-02-12 14:24:09 +00:00
|
|
|
|
2022-07-26 20:48:47 +00:00
|
|
|
BWAMEM2_INDEX ( [ [:], fasta ] )
|
2021-12-08 14:31:27 +00:00
|
|
|
BWAMEM2_MEM ( input, BWAMEM2_INDEX.out.index, false )
|
2021-02-12 14:24:09 +00:00
|
|
|
}
|
|
|
|
|
2021-12-08 14:31:27 +00:00
|
|
|
//
|
|
|
|
// Test with single-end data and sort
|
|
|
|
//
|
|
|
|
workflow test_bwamem2_mem_single_end_sort {
|
|
|
|
input = [
|
|
|
|
[ id:'test', single_end:true ], // meta map
|
|
|
|
[
|
|
|
|
file(params.test_data['sarscov2']['illumina']['test_1_fastq_gz'], checkIfExists: true)
|
|
|
|
]
|
|
|
|
]
|
|
|
|
fasta = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
|
|
|
|
|
2022-07-26 20:48:47 +00:00
|
|
|
BWAMEM2_INDEX ( [ [:], fasta ] )
|
2021-12-08 14:31:27 +00:00
|
|
|
BWAMEM2_MEM ( input, BWAMEM2_INDEX.out.index, true )
|
|
|
|
}
|
|
|
|
|
|
|
|
|
2021-05-12 13:56:46 +00:00
|
|
|
//
|
|
|
|
// Test with paired-end data
|
|
|
|
//
|
2021-02-12 14:24:09 +00:00
|
|
|
workflow test_bwamem2_mem_paired_end {
|
2021-11-26 07:58:40 +00:00
|
|
|
input = [
|
|
|
|
[ id:'test', single_end:false ], // meta map
|
|
|
|
[
|
|
|
|
file(params.test_data['sarscov2']['illumina']['test_1_fastq_gz'], checkIfExists: true),
|
|
|
|
file(params.test_data['sarscov2']['illumina']['test_2_fastq_gz'], checkIfExists: true)
|
|
|
|
]
|
|
|
|
]
|
2021-03-25 09:24:21 +00:00
|
|
|
fasta = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
|
2021-02-12 14:24:09 +00:00
|
|
|
|
2022-07-26 20:48:47 +00:00
|
|
|
BWAMEM2_INDEX ( [ [:], fasta ] )
|
2021-12-08 14:31:27 +00:00
|
|
|
BWAMEM2_MEM ( input, BWAMEM2_INDEX.out.index, false )
|
|
|
|
}
|
|
|
|
|
|
|
|
//
|
|
|
|
// Test with paired-end data and sort
|
|
|
|
//
|
|
|
|
workflow test_bwamem2_mem_paired_end_sort {
|
|
|
|
input = [
|
|
|
|
[ id:'test', single_end:false ], // meta map
|
|
|
|
[
|
|
|
|
file(params.test_data['sarscov2']['illumina']['test_1_fastq_gz'], checkIfExists: true),
|
|
|
|
file(params.test_data['sarscov2']['illumina']['test_2_fastq_gz'], checkIfExists: true)
|
|
|
|
]
|
|
|
|
]
|
|
|
|
fasta = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
|
|
|
|
|
2022-07-26 20:48:47 +00:00
|
|
|
BWAMEM2_INDEX ( [ [:], fasta ] )
|
2021-12-08 14:31:27 +00:00
|
|
|
BWAMEM2_MEM ( input, BWAMEM2_INDEX.out.index, true )
|
2021-02-12 14:24:09 +00:00
|
|
|
}
|