I am trying to replicate the results of the following paper: https://genome.cshlp.org/content/28/3/396
The data is long reads from Pacbio RS II and is available in SRA in the following location : https://www.ncbi.nlm.nih.gov/Traces/study/?acc=SAMN06328616&o=acc_s%3Aa
I am getting an error "ERROR, could not read quality metrics for FASTQ Sequence" when I try to use the program bax2bam to convert the legacy pacbio files to pacbio bam files.
Any suggestions how to convert the files from SRA to pacbio bam files to be used with the IsoSeq3 pipeline?
Are you using the original data files found under the
Data Accesstab of a the samples? A representative example from your data here. If not, try those
Yes I am trying out the bax files from the Data Access tab. Here is an example file name.
As I mentioned I got the error with these bax files. I have also tried to remove the .1 after the .h5 but still getting the same error. "ERROR, could not read quality metrics for FASTQ Sequence"
Have you tried multiple files and do they all get the same error?
I am getting the error for multiple files.
I have the same problem with bax.h5 files we generated in 2017. I'm seeing this enough to think it's a bug, rather than a problem with input files.
Interesting. Perhaps contacting PacBio tech support may be the only option then.