Hello, I am wondering why the error pattern of PacBio raw data is dominated by InDel. Could someone help to explain it ?
Thanks in advance.
Thanks in advance.
-
-
I think it's due to both the 'single molecule' and 'real time' nature of the system. The RS is essentially taking a video of the polymerase adding nucleotides in real time. As this happens very quickly, the imaging system might miss an addition (or maybe an unlabeled nucleotide slipped into the system), making it look like a small deletion. I'm not sure what would lead to a small insert (but I'm also not sure that error model really happens with PacBio).Originally posted by juckdnarocks View Post -
The insertion errors result from a nucleotide entering the detection zone for a significant amount of time, but without being incorporated.Comment
-
-
Indeed. The measured signal is based on the residence time of the tagged nucleotide rather than from the tagged end of a terminal incorporation. You can also have problems with pulses merging together when you hit a long homopolymerComment
-
It means that SMRT has the problem with a long homopolymer like 454? I thought it is a random error model.Originally posted by ELoomis View PostComment
-
"Problem" is relative. You'll have a harder time calling the exact number of nt's in the homopolymer run, but polymerase keeps running through it, so you can still get accuracy improvement with CCS (if you really want to know how many nt's are there) or flanking sequence on either side (if you just want to map your read). If the exact number of nt's in a homopolymer run is your thing, you could also delve into the basecalling parameters to improve/optimize since this isn't a major priority for the default user...Originally posted by juckdnarocks View Post
In my experience, the SMRT error profile is remarkably stable through very extreme sequence compositions (100% GC, trinucleotide repeats) and all the way to the end of the raw read.Comment
-
Homopolymers
PacBio's sequencing errors in homopolymers are still stochastic (random), just at a higher rate. With enough coverage, the consensus across homopolymers approaches 100%, just like in non-homopolymer regions. This is the case for de novo assembly, resequencing, and also single molecules (circular consensus).
By contrast, systematic errors don't go away with coverage, and limit the ultimate consensus accuracy. That's why many sequencing experiments plateau at Q40 or so. By contrast, because of the randomness of the errors, PacBio has demonstrated results greater than Q50 for a range of bacterial genomes and BACs.
See this blog from PacBio for a better explanation:
Comment
-
Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?
Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating...10-14-2026, 09:55 AM -
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...
Comment