Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • How to evaluate the Polymerase read and subread statistics from the raw data?

    Hi

    How can I check before assembly that the PacBio data provided by service provider is good. What I have received from my service provider is the raw data and two types of statistics:
    Polymerase read statistics
    Subread statistics

    Thanks
    Last edited by swatie; 01-10-2018, 06:36 PM. Reason: Want to add type of data

  • #2
    The subread statistics are the important numbers for assembly. You should have enough coverage in the subread bases, and good subread mean / N50 readlingth.

    The subread stats will depend on the library quality (how much long DNA was in the sample) as well as generally sequencing variables.

    Comment


    • #3
      Hi, Thank you for the reply.
      So the subread statistics are;
      Subread bases Subreads Subread N50 Avg. subread length
      sample1_cell1 1,420,065,796 163,879 12166 8665
      sample1_cell2 1,314,563,766 150,564 12,399 8,730

      so the average subread length/ N50 subread for cell1 is 8665/12166 (0.71)
      and for cell2 is 8730/12,300 (0.70)
      Are these good? What is the threshold to decide these values good or bad ? (average genome size of such strains is 36Mb)

      In addition, I ran the RS_Subreads protocol from the SMRT portal from the raw data. The loading P1 values for cell is 77.45% (P0 is 11.06% and P2 is 11.49%) and for cell2 is 73.59% (P0 is 9.48% and P2 is 16.94%). In general, I guess a good P1 is between 30-40%, greater than 40% value of P1 gives a large number of shorter contigs, right ?

      ( The Chemistry used was P6-C4)

      Comment


      • #4
        The P1 is high, probably a sign of overloading. This results in low quality data, and limited subread lengths.
        It's always worth assembling the data, even if the raw data isn't of the absolute highest quality. Simply run HGAP with 36Mb as the estimated genome size.
        With the preassembly stats and assembly results it should be possible to estimate how much the data quality has effected the results.

        Comment


        • #5
          Dear Dr. Hall,
          We have a similar problem than swatisinha. We got sequel data for a de novo genome assembly, in total we have 30Gb of sub-reads with mean size 8.7kb. Our lab outsourced a bioinformatics service company to run Falcon, but after trying all the parameters the assembly is not good, with N50 around 70kb. We have discarded heterozygosity and contamination. We believe our raw data is flawed and want to test this, Would you, by any chance, have any advice on how to measure the quality of our reads?

          Here are some of the stats from the sequencing facility. Also, Well A01 was done on a different Run than the rest.
          Samples -- Yield (GB) -- Pol_RDL_mean -- sub-RDL mean -- ZMW_0 % -- ZMW_1 % -- ZMW_2 % -- Pippin cutoff
          Well A01 -- 2.55 -- 13,144 -- 10,763 -- 70.32 -- 19.44 -- 10.24 -- 15kb
          Well B01 -- 6.08 -- 8,585 -- 7,640 -- 13.10 -- 69.40 -- 17.50 -- 15kb
          Well C01 -- 7.13 -- 9,759 -- 8,372 -- 13.70 -- 71.70 -- 14.60 -- 15kb
          Well D01 -- 7.55 -- 10,546 -- 8,873 -- 17.40 -- 70.20 -- 12.40 -- 15kb
          Well E01 -- 6.87 -- 9,046 -- 7,908 -- 11.30 -- 74.50 -- 143 -- 15kb

          Thanks so much for any insight you may have.
          Last edited by am33567; 03-19-2018, 06:07 PM. Reason: Spaces

          Comment


          • #6
            Assembly stats is not indicator of sequence quality as it will depend on input DNA heterozygosity, genome composition, assembly software and coverage among other factors.

            Well A01 is underloaded and others overloaded. But your sequencing QC looks Sequel average. In PacBio systems read length and yield generally have inverse correlation. If they were loaded according to PacBio recommendation of P1~35-45% you would get around 1Kb longer subreads for B-E wells in the expense of yield and 1kb shorter for A01 with increased yield.

            Comment


            • #7
              The data looks a little overloaded, but for sequel these number do not look unreasonable.
              The mean subread size does look low, but given a 15kb size selection, it's within expectations. Was there issues with the initial DNA quality? For denovo assembly a much higher size selection would be recommended.
              How did you
              discarded heterozygosity and contamination
              ?
              The first place I would be looking in this case is the assembly, not the raw data. What do the assembly parameters look like, preassembeld seed cut, yield, N50, coverage, Draft assembly size, repeat content etc.

              Comment

              Latest Articles

              Collapse

              • 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
              • seqadmin
                The Impact of AI in Genomic Medicine
                by seqadmin



                Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                02-26-2024, 02:07 PM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 03-14-2024, 06:13 AM
              0 responses
              34 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-08-2024, 08:03 AM
              0 responses
              72 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-07-2024, 08:13 AM
              0 responses
              82 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-06-2024, 09:51 AM
              0 responses
              68 views
              0 likes
              Last Post seqadmin  
              Working...
              X