-
Notifications
You must be signed in to change notification settings - Fork 55
Description
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!