2020-11-25 00:16:56 +00:00
|
|
|
#!/usr/bin/env nextflow
|
|
|
|
|
|
|
|
nextflow.enable.dsl = 2
|
|
|
|
|
2021-11-26 07:58:40 +00:00
|
|
|
include { BWA_INDEX } from '../../../../modules/bwa/index/main.nf'
|
|
|
|
include { BWA_MEM } from '../../../../modules/bwa/mem/main.nf'
|
2020-11-25 00:16:56 +00:00
|
|
|
|
2021-05-12 13:56:46 +00:00
|
|
|
//
|
|
|
|
// Test with single-end data
|
|
|
|
//
|
2020-11-25 00:16:56 +00:00
|
|
|
workflow test_bwa_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)
|
2020-11-25 00:16:56 +00:00
|
|
|
|
2021-03-24 09:53:41 +00:00
|
|
|
BWA_INDEX ( fasta )
|
2021-12-08 14:31:27 +00:00
|
|
|
BWA_MEM ( input, BWA_INDEX.out.index, false )
|
|
|
|
}
|
|
|
|
|
|
|
|
//
|
|
|
|
// Test with single-end data and sort
|
|
|
|
//
|
|
|
|
workflow test_bwa_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)
|
|
|
|
|
|
|
|
BWA_INDEX ( fasta )
|
|
|
|
BWA_MEM ( input, BWA_INDEX.out.index, true )
|
2020-11-25 00:16:56 +00:00
|
|
|
}
|
|
|
|
|
2021-05-12 13:56:46 +00:00
|
|
|
//
|
|
|
|
// Test with paired-end data
|
|
|
|
//
|
2020-11-25 00:16:56 +00:00
|
|
|
workflow test_bwa_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)
|
2020-11-25 00:16:56 +00:00
|
|
|
|
2021-03-24 09:53:41 +00:00
|
|
|
BWA_INDEX ( fasta )
|
2021-12-08 14:31:27 +00:00
|
|
|
BWA_MEM ( input, BWA_INDEX.out.index, false )
|
|
|
|
}
|
|
|
|
|
|
|
|
//
|
|
|
|
// Test with paired-end data and sort
|
|
|
|
//
|
|
|
|
workflow test_bwa_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)
|
|
|
|
|
|
|
|
BWA_INDEX ( fasta )
|
|
|
|
BWA_MEM ( input, BWA_INDEX.out.index, true )
|
2020-11-25 00:16:56 +00:00
|
|
|
}
|