Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • seqqeq
    Junior Member
    • Nov 2009
    • 3

    Coverage requirement for denovo sequencing

    Is there any study to explore the relationship between read length, raw read accuracy, coverage, genome size and the resulting number of contigs for denovo assembly? I am particularly interested in human genome. For example, if we have 10k reads, do we need less throughput to achieve the same quality assembly?

    Appreciate someone give me some pointers. Simulation is fine too.
  • krobison
    Senior Member
    • Nov 2007
    • 734

    #2
    A place to start would be Lander-Waterman statistics


    however, these don't incorporate the actual structure of the human genome -- notably the repeats. These will favor long reads even more.

    Comment

    • seqqeq
      Junior Member
      • Nov 2009
      • 3

      #3
      Originally posted by krobison View Post
      A place to start would be Lander-Waterman statistics


      however, these don't incorporate the actual structure of the human genome -- notably the repeats. These will favor long reads even more.
      Thanks for the pointer. Its a very useful set of notes. In the NIH estimate for human genome sequencing cost:

      The following 'sequence coverage' values were used in calculating the cost per genome:

      Sanger-based sequencing (average read length=500-600 bases): 6-fold coverage
      454 sequencing (average read length=300-400 bases): 10-fold coverage
      Illumina and SOLiD sequencing (average read length=50-100 bases): 30-fold coverage

      (http://www.genome.gov/sequencingcosts/)

      They listed different coverage numbers from that estimated from Lander-Waterman statistics. If I use genome size 3Gb, 1 contig, read lengths 600, 400 and 100 bps, the formula gives coverage of 20.2, 18.7 and 18.3. (The numbers above are 30, 10 and 6). Any idea why its way off, esp for Sanger?

      Comment

      • jebowers
        Member
        • Apr 2013
        • 19

        #4
        I would assume the NIH estimates take into account repeats and error rates as well as sequencing biases.

        Different sequencing technologies favor different GC% sequences, and to get the 99% coverage some technologies must sequence more to get the needed coverage of GC% extremes.

        Also repeats complicate matters, as short sequences containing repeats can not be effectively anchored, requiring higher coverage.

        Finally error rates result in some loss of coverage, again depending on the tech used.

        So in summary, yes a formula would be great, but guidelines based on experimental results are probably going to be more accurate as several different factors besides read length strongly impact the results.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Pathogen Surveillance with Advanced Genomic Tools
          by seqadmin




          The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
          03-24-2025, 11:48 AM
        • seqadmin
          New Genomics Tools and Methods Shared at AGBT 2025
          by seqadmin


          This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

          The Headliner
          The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
          03-03-2025, 01:39 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 03-20-2025, 05:03 AM
        0 responses
        49 views
        0 reactions
        Last Post seqadmin  
        Started by seqadmin, 03-19-2025, 07:27 AM
        0 responses
        57 views
        0 reactions
        Last Post seqadmin  
        Started by seqadmin, 03-18-2025, 12:50 PM
        0 responses
        49 views
        0 reactions
        Last Post seqadmin  
        Started by seqadmin, 03-03-2025, 01:15 PM
        0 responses
        200 views
        0 reactions
        Last Post seqadmin  
        Working...