forked from nextflow-io/patterns
-
Notifications
You must be signed in to change notification settings - Fork 0
/
multiple_inputs.nf
82 lines (69 loc) · 2.01 KB
/
multiple_inputs.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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
#!/usr/bin/env nextflow
/*
*
* Copyright (c) 2013-2018, Centre for Genomic Regulation (CRG).
* Copyright (c) 2013-2018, Paolo Di Tommaso and the respective authors.
*
* This file is part of 'Nextflow'.
*
* Nextflow is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* Nextflow is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with Nextflow. If not, see <http://www.gnu.org/licenses/>.
*/
/*
*
*
* Try to run this example by entering the following command line:
*
* nextflow run multiple_inputs.nf --query data/prot_\*.fa
*
* It will execute a blast search for each file that matches the
* the wildcard path matcher
*
*/
/*
* Define the pipeline default parameters.
*/
params.query = "$baseDir/data/sample.fa"
params.db = "$baseDir/blast-db/pdb/tiny"
params.out = 'blast_result.txt'
/*
* the path where the BLAST DB is located
*/
db_name = file(params.db).name
db_path = file(params.db).parent
/*
* A channel emitting the fasta files specified as parameter
*/
seq = Channel.fromPath(params.query)
/*
* Execute a BLAST job for each chunk for the provided sequences
*/
process blast {
tag "$x.baseName"
input:
file x from seq
file db_path
output:
file 'out' into blast_result
"""
blastp -db $db_path/$db_name -query $x -outfmt 6 > out
"""
}
/*
* Collect all the outputs produced by the `blast` process
* executions to a single file, whose name is defined by
* the `params.out` parameters
*/
blast_result
.collectFile(name: file(params.out))
.println { "Result saved to file: $it" }