-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmain.nf
37 lines (30 loc) · 939 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
37
nextflow.preview.output = true
include { HOMER_GETMAPPABLEREGIONS } from './modules/local/homer/getmappableregions'
include { HOMER_CREATEUNIQMAP } from './modules/local/homer/createuniqmap'
workflow {
main:
ch_versions = Channel.empty()
// Split FASTA by chromosome
split_fastas = file(params.fasta).splitFasta(by: 1, file: true)
// Generate mappable regions
HOMER_GETMAPPABLEREGIONS(
split_fastas,
params.parallel_sequences_analyzed,
params.read_length
)
ch_versions = ch_versions.mix(HOMER_GETMAPPABLEREGIONS.out.versions)
// Create uniqmap directory
HOMER_CREATEUNIQMAP(
HOMER_GETMAPPABLEREGIONS.out.txt,
params.genome_name,
params.read_length
)
ch_versions = ch_versions.mix(HOMER_CREATEUNIQMAP.out.versions)
publish:
HOMER_CREATEUNIQMAP.out.uniqmap_zip >> 'homer'
}
output {
'homer' {
mode 'copy'
}
}