Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • minion error rate

    how can I know the error rate of my reads Minion?
    how to distinguish the rate of insertion, deletion and substitutions?
    thanks

  • #2
    You won't unless there is a reference available to compare the reads to.

    Comment


    • #3
      for now I am trying to look for the reference.
      If there is other solutions to determine the rate of error thank you for telling it to me.

      Comment


      • #4
        If there's no reference, a reasonable ballpark (especially for passed reads) is about 10% total error, distributed about 1/3 insertions, 1/3 deletions, and 1/3 SNPs.

        It's possible to remap reads against other reads, but you need to be very careful with that, particularly across homopolymer regions. There's a lot of systematic error in the base calling (not the signal) which causes problems in consensus sequences.

        Comment


        • #5
          according to my research the overall error rate in the Minion sequencer is about 25%.
          what's homopolymer regions?
          Thank you

          Comment


          • #6
            Homo-polymers are long stretches of identical nucleotides (e.g. AAAAAAAAAAAAAAA)

            This paper claims an error rate of 38%, but it is early days. This rate may be dependent on the sample and its composition too.

            Comment


            • #7
              I want to make a local alignment Minion reads but I do not know how to put the parameters of: insertion, deletion and match / mismatch.
              I want to use "Affine GAP PENALTIES" but this method penalizes insertions and deletions. If I'm going to use "Affine GAP PENALTIES" but I also need a function that penalizes subsitutions.
              Is this is a good idea friends?

              Comment


              • #8
                according to my research the overall error rate in the Minion sequencer is about 25%.
                When was the sequencing done? The f1000 paper has the most recent error analysis that I'm aware of (total error 10-15%), and that is from sequencing done a couple of months ago, prior to changing to a hexamer model for base calling:

                Read the original article in full on F1000Research: MinION Analysis and Reference Consortium: Phase 1 data release and analysis

                Comment


                • #9
                  How high is the rate of "phantasy sequences", that have no resemblance to the reference, with the latest versions? I actually can't find many data of interest for me in the f1000 article. The athors are writing often a bout "a run". Are they writing about an average run or the cherry picked best run they ever encountered?

                  Comment


                  • #10
                    How high is the rate of "phantasy sequences", that have no resemblance to the reference, with the latest versions?
                    Other people call this "mismatch rate", or alternatively refer to "mapping rate". From the 2D reads, mapping rate was pretty close to 100% for all except two runs (see Figure 6, or section with heading "Proportion of target and control sample"). I'm a bit hesitant to trust this fully, because there's a high chance of false positive matches when you have a high error rate.

                    The authors are writing often about "a run". Are they writing about an average run or the cherry picked best run they ever encountered?
                    There were 20 runs which each tried to stick to the same specific sample preparation and sequencing protocol. You can call those cherry picked runs if you like, but these were the only runs that the groups did as part of their Phase 1 experiments. Only six runs were able to keep to this without variation, with the deviations from the standard protocol specified in Table S5.

                    Comment


                    • #11
                      Originally posted by gringer View Post
                      When was the sequencing done? The f1000 paper has the most recent error analysis that I'm aware of (total error 10-15%), and that is from sequencing done a couple of months ago, prior to changing to a hexamer model for base calling:

                      http://f1000research.com/articles/4-1075/v1
                      So now is hexamer model? How much base calling accuracy does it gain over the previous pentamer model?

                      Comment


                      • #12
                        Originally posted by ymc View Post
                        So now is hexamer model? How much base calling accuracy does it gain over the previous pentamer model?
                        See Jared Simpson's post here. There's no raw sequence accuracy mentioned there, but in consensus it was a 0.4% improvement.

                        I don't have results from any comparative sequencing that we've done, because getting *any* sequence is a bit of a challenge for us.

                        Comment


                        • #13
                          Consensus improvement isn't the same as per-base improvement. I wonder what the difference is there.

                          Also Jared said it was mainly down to improvements in homopolymers or whether it is simply due to the frequency distribution of different homopolymer lengths. I wonder if this means they've changed their strategy.

                          Obviously any homopolymer longer than 5 (now 6) looks like a single event. Eg TAAAAAAAAG yields TAAAA AAAAA AAAAA AAAAA AAAAA AAAAG. Those AAAAA events all get joined together into a single longer event, which meant the longest homopolymer previously reported was 5. I assume it's now 6.

                          However an alternative (and obvious) strategy is to observe the event length distributions for non-homopolymer. They're largely random, centred around a particular amount. They'll either be gaussian or poisson I'd guess. Given a homopolymer signal of length L we can hypothesise 5, 6, 7, 8, ... lengths and derive the probability of the homopolymer being more than 5. The error rate would likely be horrendous, but in theory it ought to be better than just say "5, never more".

                          No idea if the signal is strong enough (ie the distribution tight enough) to extract sufficient accuracy to make it anything better than a wild guess.

                          Comment


                          • #14
                            Originally posted by jkbonfield View Post
                            However an alternative (and obvious) strategy is to observe the event length distributions for non-homopolymer. They're largely random, centred around a particular amount. They'll either be gaussian or poisson I'd guess. Given a homopolymer signal of length L we can hypothesise 5, 6, 7, 8, ... lengths and derive the probability of the homopolymer being more than 5. The error rate would likely be horrendous, but in theory it ought to be better than just say "5, never more".
                            This doesn't work so well in the way that ONT is modelling for base-calling, because the event length actually depends quite strongly on the bases that are found about 20bp upstream from the signal site (where the DNA is being unwound and split). In other words, it's not particularly random.

                            Comment


                            • #15
                              Originally posted by gringer View Post
                              This doesn't work so well in the way that ONT is modelling for base-calling, because the event length actually depends quite strongly on the bases that are found about 20bp upstream from the signal site (where the DNA is being unwound and split). In other words, it's not particularly random.
                              In that case, it makes it possible to be more accurate than just guessing based on the random distribution, albeit perhaps also (too?) complex to tease out all the correlations.

                              Either way, I'm curious to know what the per-sequence-base (not per consensus base) accuracy is like with a 6mer model vs a 5mer model. We have some data of our own, but fundamentally the variabiity from run to run is high enough that I think it would need a large project to average out that variability to get robust numbers.

                              Comment

                              Latest Articles

                              Collapse

                              • seqadmin
                                Strategies for Sequencing Challenging Samples
                                by seqadmin


                                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                                03-22-2024, 06:39 AM
                              • seqadmin
                                Techniques and Challenges in Conservation Genomics
                                by seqadmin



                                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                                Avian Conservation
                                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                                03-08-2024, 10:41 AM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by seqadmin, Yesterday, 06:37 PM
                              0 responses
                              12 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, Yesterday, 06:07 PM
                              0 responses
                              10 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-22-2024, 10:03 AM
                              0 responses
                              51 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-21-2024, 07:32 AM
                              0 responses
                              68 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X