Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Advanced coverage predictor

    Hi all,

    Is there some kind of a Poisson calculator out there that will predict not only mean coverage but coverage at various percents? For example, given a human genome, what mean coverage do I need to get 95% of my bases covered at greater than 30X?

    Thanks!!

    Seth

  • #2
    Use ppois in R

    If I'm not mistaken, what you're asking is essentially: if I assume a Poisson distribution, what mean rate do I need to have for at least 95% of the area under my curve above 29?

    The answer to that can be found using a simple lapply or for loop in R. Here's what the for loop would look like:

    Code:
    for(i in 30:75) {
     # ppois returns the density under the curve, lower.tail=F
     # specifies that you want to return the area under the curve
     # at or above 30
      cat(i,ppois(29,i,lower.tail=F),"\n")
    }
    This prints out the mean coverage, and the percentage of bases with at least 30X coverage. The answer to your question appears to be that you should aim for mean coverage of at least 40X.

    The code isn't that great, but it should give you an idea of how to answer the question.

    Comment


    • #3
      I made this. Does it look right?



      It suggests that to get >95% 30-fold coverage I need an average coverage of 41...can that be right??
      Last edited by scrosby; 01-18-2014, 12:54 PM.

      Comment


      • #4
        Yup! I think I got that you need at least 40X average coverage.

        I haven't done much WGS, mainly exome sequencing, but I'd also bear in mind that due to GC content, the genome you're sequencing, the instrument you're using, and other factors, you may not see 95% of the genome with at least 30X coverage.

        Comment


        • #5
          scrosby: Fixed it for you using imgur.com. Sent you a PM.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Essential Discoveries and Tools in Epitranscriptomics
            by seqadmin


            The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
            Yesterday, 07:01 AM
          • seqadmin
            Current Approaches to Protein Sequencing
            by seqadmin


            Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
            04-04-2024, 04:25 PM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 04-11-2024, 12:08 PM
          0 responses
          39 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 10:19 PM
          0 responses
          41 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 09:21 AM
          0 responses
          35 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-04-2024, 09:00 AM
          0 responses
          55 views
          0 likes
          Last Post seqadmin  
          Working...
          X