-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmain.nf
36 lines (29 loc) · 775 Bytes
/
main.nf
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
// Nextflow DSL2
nextflow.enable.dsl=2
process kir_mapper_map {
tag "$sampleId"
publishDir 'output/map'
input:
tuple val(sampleId), file(r1), file(r2)
file config
file kir_mapper_db
output:
path "output/map/$sampleId"
script:
"""
wd=\$(pwd)/output
kir-mapper map \
-r1 $r1 \
-r2 $r2 \
-sample $sampleId \
-config $config \
-output \$wd \
--exome
"""
}
workflow {
fq = Channel.fromFilePairs('wes/*{1,2}.merged.fastq.gz', flat: true)
kir_mapper_db = file('kir-mapper_db_latest')
config_file = file("${workflow.projectDir}/assets/.kir-mapper_config")
kir_mapper_map(fq, config_file, kir_mapper_db)
}