This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Passing parameters to gyppy_duplex pipeline

Hi,

I am using ont_guppy_duplex_pipeline to re-base-call raw data from ONT sequencing. I am invoking full pipeline by simplex base-calling in fast mode, then duplex calling with sup mode and finally base-calling of all reads that were not in simplex base-calling with sup mode again. The problem is that it seems the duplex base-calling does not fetch the correct parameters and aborts with an error. Normally, the duplex base-calling requires two parameters:

--duplex_pairing_mode from_pair_list --duplex_pairing_file <text pair file>

to function but to my understanding there is no where in the command arguments to set that. Here is the output of help:

guppy_duplex -h usage: guppy_duplex [-h] [--basecaller_exe BASECALLER_EXE] [--duplex_basecaller_exe DUPLEX_BASECALLER_EXE] --input_path INPUT_PATH
                --save_path SAVE_PATH [--simplex_config SIMPLEX_CONFIG] [--duplex_config DUPLEX_CONFIG] [--disable_logging]
                [--skip_simplex] [--skip_duplex] [--call_non_duplex_reads] [--device DEVICE] [--do_read_splitting]
                [--duplex_chunks_per_runner DUPLEX_CHUNKS_PER_RUNNER] [--recursive]

And here is my code for running the pipe that breaks when duplex calling starts. I must add that simplex calling completes successfully.

 guppy_duplex -i ./pod5/ -s ./output2/ --simplex_config dna_r10.4.1_e8.2_400bps_fast.cfg --duplex_config  dna_r10.4.1_e8.2_400bps_sup.cfg --device 'cuda:0'

This is how the pipeline fails:

enter image description here

I appreciate if someone with the knowledge, helps with this.

base-calling ont oxford-nanopore guppy

0 answers

No answers yet.

Log in to answer this question.