Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 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!

  • #2
    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.

    Comment


    • #3
      You can find useful PacBio software here (including blasr mentioned by mchaisso): https://github.com/PacificBiosciences

      Comment


      • #4
        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, 02:51 AM.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Genetic Variation in Immunogenetics and Antibody Diversity
          by seqadmin



          The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
          11-06-2024, 07:24 PM
        • seqadmin
          Choosing Between NGS and qPCR
          by seqadmin



          Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
          10-18-2024, 07:11 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 11-08-2024, 11:09 AM
        0 responses
        205 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 11-08-2024, 06:13 AM
        0 responses
        151 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 11-01-2024, 06:09 AM
        0 responses
        80 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-30-2024, 05:31 AM
        0 responses
        26 views
        0 likes
        Last Post seqadmin  
        Working...
        X