-
Notifications
You must be signed in to change notification settings - Fork 8
/
trimming_aligning4.nf
58 lines (43 loc) · 1.08 KB
/
trimming_aligning4.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
#!/usr/bin/env nextflow
/* Testing for NCBI SpEW project, trimming module
* Zhou (Ark) Fang [email protected]
*/
inPath = "/zfs1/ncbi-workshop/AP/nextflow/"
inFiles = "$inPath/*.fastq.gz"
singleEnd = true
adapter1 = "ADATPER_FWD"
adapter2 = "ADATPER_REV"
inFiles
params.in = "/zfs1/ncbi-workshop/AP/RNAseq/*.fastq"
sequences = file(params.in)
/* records = "/zfs1/ncbi-workshop/AP/RNAseq/AlignedReads/" */
process modules{
script"
"""
bash module load bowtie2
bash module load tophat
bash module load cutadapt
"""
}
process trimming{
input:
file inFiles
output:
file 'trimmed_*' /* into records */
script:
if (singleEnd==true)
"""
bash trimming.sh -i inPath -s -a1 adapter1 -a2 adapter2
"""
else
"""
bash /zfs1/nci-workshop/AP/nextflow/trimming.sh -i inPath -a1 adapter1 -a2 adapter2
"""
}
process alignment {
input:
file '*.fastq.gz'
output:
file 'aligned_*' into /* records */
"""
bash /zfs1/ncbi-workshop/AP/nextflow/align2.sh