2021-03-23 13:15:26 +00:00
|
|
|
#!/usr/bin/env nextflow
|
|
|
|
|
|
|
|
nextflow.enable.dsl = 2
|
|
|
|
|
2021-11-26 07:58:40 +00:00
|
|
|
include { GATK4_APPLYBQSR } from '../../../../modules/gatk4/applybqsr/main.nf'
|
2021-03-23 13:15:26 +00:00
|
|
|
|
|
|
|
workflow test_gatk4_applybqsr {
|
2021-03-24 09:53:41 +00:00
|
|
|
input = [ [ id:'test' ], // meta map
|
2021-03-26 11:46:46 +00:00
|
|
|
file(params.test_data['sarscov2']['illumina']['test_paired_end_sorted_bam'], checkIfExists: true),
|
2021-04-13 17:03:09 +00:00
|
|
|
file(params.test_data['sarscov2']['illumina']['test_paired_end_sorted_bam_bai'], checkIfExists: true),
|
2021-03-26 11:46:46 +00:00
|
|
|
file(params.test_data['sarscov2']['illumina']['test_baserecalibrator_table'], checkIfExists: true)
|
2021-03-24 09:53:41 +00:00
|
|
|
]
|
2021-03-26 11:46:46 +00:00
|
|
|
fasta = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
|
|
|
|
fai = file(params.test_data['sarscov2']['genome']['genome_fasta_fai'], checkIfExists: true)
|
|
|
|
dict = file(params.test_data['sarscov2']['genome']['genome_dict'], checkIfExists: true)
|
2021-03-24 09:53:41 +00:00
|
|
|
|
|
|
|
GATK4_APPLYBQSR ( input, fasta, fai, dict, [] )
|
2021-03-23 13:15:26 +00:00
|
|
|
}
|
|
|
|
|
|
|
|
workflow test_gatk4_applybqsr_intervals {
|
2021-03-26 11:46:46 +00:00
|
|
|
input = [ [ id:'test' ], // meta map
|
|
|
|
file(params.test_data['sarscov2']['illumina']['test_paired_end_sorted_bam'], checkIfExists: true),
|
2021-04-13 17:03:09 +00:00
|
|
|
file(params.test_data['sarscov2']['illumina']['test_paired_end_sorted_bam_bai'], checkIfExists: true),
|
2021-03-26 11:46:46 +00:00
|
|
|
file(params.test_data['sarscov2']['illumina']['test_baserecalibrator_table'], checkIfExists: true)
|
|
|
|
]
|
|
|
|
fasta = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
|
|
|
|
fai = file(params.test_data['sarscov2']['genome']['genome_fasta_fai'], checkIfExists: true)
|
|
|
|
dict = file(params.test_data['sarscov2']['genome']['genome_dict'], checkIfExists: true)
|
|
|
|
intervals = file(params.test_data['sarscov2']['genome']['test_bed'], checkIfExists: true)
|
2021-03-24 09:53:41 +00:00
|
|
|
|
|
|
|
GATK4_APPLYBQSR ( input, fasta, fai, dict, intervals )
|
2021-03-23 13:15:26 +00:00
|
|
|
}
|
2021-10-29 11:01:05 +00:00
|
|
|
|
|
|
|
workflow test_gatk4_applybqsr_cram {
|
|
|
|
input = [ [ id:'test' ], // meta map
|
|
|
|
file(params.test_data['homo_sapiens']['illumina']['test_paired_end_sorted_cram'], checkIfExists: true),
|
|
|
|
file(params.test_data['homo_sapiens']['illumina']['test_paired_end_sorted_cram_crai'], checkIfExists: true),
|
|
|
|
file(params.test_data['homo_sapiens']['illumina']['test_baserecalibrator_table'], checkIfExists: true)
|
|
|
|
]
|
|
|
|
fasta = file(params.test_data['homo_sapiens']['genome']['genome_fasta'], checkIfExists: true)
|
|
|
|
fai = file(params.test_data['homo_sapiens']['genome']['genome_fasta_fai'], checkIfExists: true)
|
|
|
|
dict = file(params.test_data['homo_sapiens']['genome']['genome_dict'], checkIfExists: true)
|
|
|
|
intervals = file(params.test_data['homo_sapiens']['genome']['genome_bed'], checkIfExists: true)
|
|
|
|
|
|
|
|
GATK4_APPLYBQSR ( input, fasta, fai, dict, intervals )
|
|
|
|
}
|