Hi, I am new to pacbio and have 2 sets of .h5 files as output from pacbio. I am planning to use celera assembler and for that i need fastq files from .h5 files.
1) Is there any way to convert .h5 to fastq.
2) Is there any specific method to filter pacbio reads based on quality?
3) Do we combine both sets of data and then work on it for assembly?
Thanks!
3 answers
1) and 2) Use bash5tools.py
bash5tools.py --minLength 500 --readType subreads --minReadScore 0.8 --outType fastq
Depends on your dataset, but if you just sequenced 2 SMRT cells to get more coverage, then you can merge them prior to assembly.
https://github.com/PacificBiosciences/pbh5tools/blob/master/doc/index.rst
•
written
by
Biomonika (Noolean)
Have a look at dextract. It's very quick and lets you set a score cutoff. However, I think it only generates FASTA.
I've written a perl wrapper for the hdf5 library that you might find useful. It hasn't been tested on pacbio files though I have no reason to think it wouldn't read them.
•
written
by
Jean-Karim Heriche
Log in to answer this question.
-
Pacbio SRA data format error
written by shounak.chakraborty1990 •
-
Using Pacbio fastq file for Arrow polishing of Canu assembly
written by rongji •
-
Filtering low quality score bases from BAM file
written by tayyabaalvi82 •
-
Is there any written script to assemble multiple fastq file from accession id?
written by benkhedhermaryem044 •
-
PacBio: Subread Extraction with bash5tools.py
written by mehmetgoktay1989 •
-
Keep quality score when correcting with Canu assembler?
written by lizgzara •
-
Questions about using the Celera Genome Assembler for HGAP
written by tptacek3050
-
PacBio data assembly with Celera
written by mhasa006
-
Simulating PacBio long reads and generating consensus
written by bdeonovic
-
Best Software For 454/Pacbio Hybrid Assembly
written by mgalactus
The quality values are sufficiently low that reads may be artificially trimmed by celera. I've found it's best to just fake fastq from fasta with high enough quality value that reads are retained. The assembly quality needs to be improved later using quiver.
Could you please tell me If bash5tools.py also removes adapter sequences?
I have just used it and got subreads from raw data but I am not sure whether subreads still contains adapter sequences?
You need to post this as another question, also please refer to the manual