-
Notifications
You must be signed in to change notification settings - Fork 2
/
basic.nf
67 lines (54 loc) · 1.29 KB
/
basic.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
#!/usr/bin/env nextflow
/*
* Copyright 2020-2022, Seqera Labs
* Copyright 2013-2019, Centre for Genomic Regulation (CRG)
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
nextflow.enable.dsl=1
/*
* Command line input parameter
*/
params.in = "$baseDir/data/sample.fa"
/*
* Define the input file
*/
sequences = file(params.in)
/*
* split a fasta file in multiple files
*/
process splitSequences {
input:
file 'input.fa' from sequences
output:
file 'seq_*' into records
"""
awk '/^>/{f="seq_"++d} {print > f}' < input.fa
"""
}
/*
* Simple reverse the sequences
*/
process reverse {
input:
file x from records
output:
stdout result
"""
cat $x | rev
"""
}
/*
* print the channel content
*/
result.subscribe { println it }