Hi there, I'm currently developing a pipeline for ...
# modules
c
Hi there, I'm currently developing a pipeline for chimeric reads detection. Trying to implement yacrd, I stumbled across an issue with minimap2: I need to map the reads back onto themselves. If I do this:
Copy code
MINIMAP2_ALIGN(reads, reads, [], [], [], [])
it fails because of a file name collision when staging. If I leave out the second "reads", it will fail because of a null path value. How can I solve this? I feel like a solution is already implemented because the "target" of minimap2 falls back to the reads from the first channel if the reference (the second parameter) is null.
f
I haven't used minimap in this way, but since the target uses the reads if reference is empty, could you not give it an empty channel?
MINIMAP2_ALIGN(reads, [[],[]], [], [], [], [])
c
Tried that, but that gives an error:
Path value cannot be null
which is weird because as I see it, it will fall back to the reads:
def target = reference ?: (bam_input ? error("BAM input requires reference") : reads)
c
I made a mistake with the brackets 🙃 Thanks! Now it works
👍 2