0
|
1 #!/usr/bin/env nextflow
|
|
2
|
|
3 python3 = "python3"
|
|
4 resfinder = "/home/projects/cge/apps/resfinder/resfinder/run_resfinder.py"
|
|
5
|
|
6 params.indir = './'
|
|
7 params.ext = '.fq.gz'
|
|
8 params.outdir = '.'
|
|
9 params.species = 'other'
|
|
10
|
|
11 println("Search pattern: $params.indir*{1,2}$params.ext")
|
|
12
|
|
13 Channel
|
|
14 .fromFilePairs("$params.indir*{1,2}$params.ext", followLinks: true)
|
|
15 .set{ infile_ch }
|
|
16
|
|
17 process resfinder{
|
|
18
|
|
19 cpus 5
|
|
20 time '30m'
|
|
21 memory '1 GB'
|
|
22 clusterOptions '-V -W group_list=cge -A cge'
|
|
23 executor "PBS"
|
|
24
|
|
25 input:
|
|
26 set sampleID, file(datasetFile) from infile_ch
|
|
27
|
|
28 output:
|
|
29 stdout result
|
|
30
|
|
31 """
|
|
32 set +u
|
|
33 module unload mgmapper metabat fastqc
|
|
34 module unload ncbi-blast perl
|
|
35 source /home/projects/cge/apps/env/rf4_env/bin/activate
|
|
36 module load perl
|
|
37 module load ncbi-blast/2.8.1+
|
|
38 if [ $params.species = 'other' ]
|
|
39 then
|
|
40 $python3 $resfinder -acq -ifq $datasetFile -o '$params.outdir/$sampleID' -s '$params.species'
|
|
41 else
|
|
42 $python3 $resfinder -acq -ifq $datasetFile -o '$params.outdir/$sampleID' -s '$params.species' --point
|
|
43 fi
|
|
44 """
|
|
45 }
|
|
46
|
|
47 /*
|
|
48 result.subscribe {
|
|
49 println it
|
|
50 }
|
|
51 */
|