SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Pacific Biosciences



Similar Threads
Thread Thread Starter Forum Replies Last Post
average quality score for fastq file jgibbons1 Bioinformatics 12 07-16-2013 01:13 PM
Detecting FASTQ Quality Score Format bhan Bioinformatics 2 01-02-2012 03:40 PM
454 quality score, z-score,.. nii 454 Pyrosequencing 2 10-14-2011 09:46 AM
ELAND output quality score kwtennis311 Bioinformatics 0 06-30-2010 09:51 AM
Fastq quliaty score and MAQ output quality score baohua100 Bioinformatics 1 02-19-2009 09:21 AM

Reply
 
Thread Tools
Old 07-12-2013, 03:57 AM   #1
vamiyn
Junior Member
 
Location: Finland

Join Date: Apr 2013
Posts: 2
Default PacBio RS quality score output definition (FASTQ)

I'm simulating reads that mimic PacBio's data (long reads, about 1k-20k bases). Other things are pretty well in order, but I'd like to know what kind of quality values there usually are and what their distribution and meaning is. I'd be very thankful if someone showed me to the right direction or explained at least some part of this!

- which characters are used? (as there seems to be many kinds of variations of FASTQ... I just can't seem to find anything about PacBio's FASTQ directly)
- is there a bigger probability for an indel (or a substitution) if the quality score is bad?
- how often do you generally see certain quality scores, where goes the line between "probably an error" and "most likely fine" ?
- how likely is it to have a bad quality score if the read itself is errorless?

Thank you in advance!
vamiyn is offline   Reply With Quote
Old 07-15-2013, 12:50 PM   #2
mchaisso
Member
 
Location: Seattle, WA

Join Date: Apr 2008
Posts: 84
Default

PacBio produces reads that have different quality values for deletion, insertion, and substitution, so the read are in HDF5 rather than FASTQ.

There is a read simluator called alchemy that is packaged with blasr that simulates HDF5 in a large number of ways using an empirical error model.
mchaisso is offline   Reply With Quote
Old 07-15-2013, 05:15 PM   #3
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,076
Default

You can find useful PacBio software here (including blasr mentioned by mchaisso): https://github.com/PacificBiosciences
GenoMax is offline   Reply With Quote
Old 07-16-2013, 02:44 AM   #4
vamiyn
Junior Member
 
Location: Finland

Join Date: Apr 2013
Posts: 2
Default

Okay, I see... I have to use FASTQ for this, so unfortunately the HDF5 reads alone don't help that much. (PacBio really doesn't produce FASTQ directly?) Thank you, though, I'll definitely check the link anyway!

But okay, now I'd need to know what happens if the reads are converted to from HDF5 to FASTQ format afterwards, what will happen to which quality value or is there a standard for that? Seems like at least some information would be lost then? Or does FASTQ's quality string have enough information (compared to HDF5) that it's worth mimicking? (Does blasr use the quality string in FASTQ or only in the HDF5?) Plus the questions in my first post, still, if FASTQ has enough info.

Last edited by vamiyn; 07-16-2013 at 02:51 AM.
vamiyn is offline   Reply With Quote
Reply

Tags
explanation, fastq, pacbio, quality scores, simulation

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:56 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO