Error executing process > 'TAILFINDR_ESTIMATE_TAIL:estimateTailSize (0)

[PraddyumnaR]

Hello Luca,
I tried running the MOP_TAIL module, but I got this error message:

ERROR ~ Error executing process > 'TAILFINDR_ESTIMATE_TAIL:estimateTailSize (0)'
Caused by:
Process TAILFINDR_ESTIMATE_TAIL:estimateTailSize (0) terminated with an error exit status (1)
Command executed:
R --vanilla --slave -e "library(tailfindr); find_tails(fast5_dir = './' , save_dir = './', , csv_filename = '0_findr.csv', num_cores = 1)"
gzip *_findr.csv
Command exit status:
1
Command output:
────────────────────────────────────────────────────────────────────────────────
── Started tailfindr (version 0.1.0) ───────────────────────────────────────────
────────────────────────────────────────────────────────────────────────────────
☰ You have configured tailfindr as following:
❯ fast5_dir: ./
❯ save_dir: ./
❯ csv_filename: 0_findr.csv
❯ num_cores: 1
❯ basecall_group: Basecall_1D_000
❯ save_plots: FALSE
❯ plot_debug_traces: FALSE
❯ plotting_library: rbokeh
── Processing started at 2024-08-07 04:20:51 ───────────────────────────────────
• Searching for all Fast5 files...
Done! Found 1 Fast5 files.
• Analyzing a single Fast5 file to assess if your data
is in an acceptable format...
✔ The data has been basecalled using Guppy.
✔ Flipflop model was used during basecalling.
✔ Every read is in a single fast5 file of its own.
✔ The experiment type is RNA, so we will search
for poly(A) tails.
✔ The reads are 1D reads.
• Starting a parallel compute cluster...
Done!
• Searching for Poly(A) tails...
Processing chunk 1 of 1

|                                                                            
|                                                                      |   0%
|                                                                            
|======================================================================| 100%

• Formatting the tail data...

Command error:
────────────────────────────────────────────────────────────────────────────────
── Started tailfindr (version 0.1.0) ───────────────────────────────────────────
────────────────────────────────────────────────────────────────────────────────
☰ You have configured tailfindr as following:
❯ fast5_dir: ./
❯ save_dir: ./
❯ csv_filename: 0_findr.csv
❯ num_cores: 1
❯ basecall_group: Basecall_1D_000
❯ save_plots: FALSE
❯ plot_debug_traces: FALSE
❯ plotting_library: rbokeh
── Processing started at 2024-08-07 04:20:51 ───────────────────────────────────
• Searching for all Fast5 files...
Done! Found 1 Fast5 files.
• Analyzing a single Fast5 file to assess if your data
is in an acceptable format...
✔ The data has been basecalled using Guppy.
✔ Flipflop model was used during basecalling.
✔ Every read is in a single fast5 file of its own.
✔ The experiment type is RNA, so we will search
for poly(A) tails.
✔ The reads are 1D reads.
• Starting a parallel compute cluster...
Done!
• Searching for Poly(A) tails...
Processing chunk 1 of 1

|                                                                            
|                                                                      |   0%
|                                                                            
|======================================================================| 100%

• Formatting the tail data...
Error in round(result$tail_length, digits = 2) :
non-numeric argument to mathematical function
Calls: find_tails
In addition: Warning messages:
1: cols is now required when using unnest().
Please use cols = c()
2: In rm(polya_fastq) : object 'polya_fastq' not found
3: Unknown or uninitialised column: tail_length.
Execution halted

Work dir:
/mnt/datalake/ngsa/analysis/bhavesh/Priyanka/project_data/work/0b/faf4ed36232de5dfcfb66fe47e1119

Tip: view the complete command output by changing to the process work dir and entering the command cat .command.out

This is my the param.config file -
params {

input_path         = "/mnt/datalake/ngsa/analysis/bhavesh/Priyanka/project_data/MOP_output/"
reference          = "/mnt/datalake/ngsa/analysis/bhavesh/Priyanka/project_data/test/GCF_000002985.6_WBcel235_genomic.fa.gz"

pars_tools         = "/mnt/datalake/ngsa/analysis/bhavesh/Priyanka/project_data/test/tools_opt.tsv"

output             = "/mnt/datalake/ngsa/analysis/bhavesh/Priyanka/project_data/MOP_tail_output"

tailfindr          = "YES"
nanopolish         = "YES"

email              = ""

}