Hey, is it possible to code a single process with different channels?
nLines_per_WES_bed = ((file(params.bed_flank_100_file).readLines().size()+50)/30).toInteger()
Channel.fromPath( params.bed_flank_100_file ).splitText( by:nLines_per_WES_bed, file:params.yy_splits_dir)
.into { WESbedSplitsA; WESbedSplitsB; WESbedSplitsC }
Channel.fromPath( params.bed_wgs_file ).splitText( by:1, file:params.yy_splits_dir)
.into { WGSbedSplitsA; WGSbedSplitsB; WGSbedSplitsC }
wgs=true
process test {
echo true
input:
file(x) from ThisChannel
script:
if wgs:
ThisChannel = WGSbedSplitsA
else:
ThisChannel = WESbedSplitsA
"""
echo $x
"""
}
I want a HaplotypeCaller with a different .bed content, depends whether it is WGS or WES.
Thanks
Damian