Skip to content

tombo preprocess annotate_raw_with_fastqs==>Added sequences to a total of 0 reads./tombo resquiggle==>OSError: Can't read data (can't open directory: /xxx/anaconda3/envs/nano/lib/hdf5/plugin) #459

@AmisareWAN

Description

@AmisareWAN

Hi!This is my first analysis of nanopore sequencing.
I'm stuck at the tombo annotation methylation step when analyzing nanopore sequencing.
Here are my two attempts, one starts by combining the fastq file with the fast5 file, using "tombo preprocess annotate_raw_with_fastqs".Another one is direct tombo analysis, using“tombo resquiggle”.But both cases are reported as errors, what should I do now to fix it?
1st:
(nano) xxx:~/xxx/tombo/test$ tombo preprocess annotate_raw_with_fastqs --overwrite --fast5-basedir /xxx/20250116/tombo/test_barcode01/fast5 --fastq-filenames /xxx/20250116/new_way/1old_guppy_fastq1/cat_fastq/barcode01.fastq --sequencing-summary-filenames /xxx/20250116/new_way/1old_guppy_fastq1/fastq/sequencing_summary.txt
[20:47:07] Getting read filenames.
[20:47:07] Parsing sequencing summary files.
******************** WARNING ********************
Some FASTQ records from sequencing summaries do not appear to have a matching file.
[20:47:18] Annotating FAST5s with sequence from FASTQs.
****** WARNING ****** Some FASTQ records contain read identifiers not found in any FAST5 files or sequencing summary files.
0it [00:03, ?it/s]
[20:47:21] Added sequences to a total of 0 reads.

2nd:
(nano) wyx@luo-2022:~/xxx/tombo/test$ tombo resquiggle /xxx/20250116/tombo/test/fast5/ /xxx/20250116/nextpolish/polished_fasta/19_rename_barcode_Salmonella/barcode01.fasta --processes 16 --overwrite --num-most-common-errors 5
[20:49:31] Loading minimap2 reference.
[20:49:31] Getting file list.
[20:49:31] Loading default canonical ***** DNA ***** model.
[20:49:31] Re-squiggling reads (raw signal to genomic sequence alignment).
5 most common unsuccessful read types (approx. %):
-----
-----
-----
-----
-----
0%| | 0Exception in thread Thread-4:
Traceback (most recent call last):
File "/xxx/anaconda3/envs/nano/lib/python3.7/threading.py", line 926, in _bootstrap_inner
self.run()
File "/xxx/anaconda3/envs/nano/lib/python3.7/threading.py", line 870, in run
self._target(*self._args, **self._kwargs)
File "/xxx/anaconda3/envs/nano/lib/python3.7/site-packages/tombo/resquiggle.py", line 166read_worker
sig_len_rng, seq_len_rng)
File "/xxx/anaconda3/envs/nano/lib/python3.7/site-packages/tombo/resquiggle.py", line 139
all_raw_signal = th.get_raw_read_slot(fast5_data)['Signal'][:]
File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper
File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper
File "/xxx/anaconda3/envs/nano/lib/python3.7/site-packages/h5py/_hl/dataset.py", line 741
return self._fast_reader.read(args)
File "h5py/_selector.pyx", line 370, in h5py._selector.Reader.read
OSError: Can't read data (can't open directory: /xxx/anaconda3/envs/nano/lib/hdf5/plugin)
Exception in thread Thread-11:
Traceback (most recent call last):
File "/xxx/anaconda3/envs/nano/lib/python3.7/threading.py", line 926, in _bootstrap_inner
self.run()
File "/xxx/anaconda3/envs/nano/lib/python3.7/threading.py", line 870, in run
self._target(*self._args, **self._kwargs)
File "/xxx/anaconda3/envs/nano/lib/python3.7/site-packages/tombo/resquiggle.py", line 166read_worker
sig_len_rng, seq_len_rng)
File "/xxx/anaconda3/envs/nano/lib/python3.7/site-packages/tombo/resquiggle.py", line 139
all_raw_signal = th.get_raw_read_slot(fast5_data)['Signal'][:]
File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper
File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper
File "/xxx/anaconda3/envs/nano/lib/python3.7/site-packages/h5py/_hl/dataset.py", line 741
return self._fast_reader.read(args)
File "h5py/_selector.pyx", line 370, in h5py._selector.Reader.read
OSError: Can't read data (can't open directory: /xxx/anaconda3/envs/nano/lib/hdf5/plugin)

Is there a step here that I haven't noticed, and how should I fix this problem!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions