Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Fad2012
    Member
    • Sep 2012
    • 62

    specific location coverage from average coverage?

    Hello there,

    I want to know if there is any way to calculate the coverage at a specific nucleotide across the amplicon using the average coverage, with the assumption that the reads are normally distributed across the ampliocn.

    Say I have sequenced a 1.3kb amplicon, with read length 300 nts, and number of reads 10000. using Lander/Waterman equation (LN/G), the average coverage will be ~2300x. This tells us that every nucleotide was covered at least with 2300 reads. My question is that, can I estimate roughly the actual number of reads cover a certain nucleotide across the ampliocn, say nucleotide number 9, or any other nucleotide. Is it doable in a theoretical way without using tools like bedtools or samtools?

    Many thanks
  • Brian Bushnell
    Super Moderator
    • Jan 2014
    • 2709

    #2
    You can't get a useful answer without mapping. Did you amplify some region then randomly shear it? If you're sequencing the amplicons directly, the coverage will all be in the same location. If you're sequencing the randomly sheared amplicons, there's no way to predict what the coverage will look like.

    Comment

    • Fad2012
      Member
      • Sep 2012
      • 62

      #3
      mmmm I see..yes this is the plan, to shear the amplicon.

      What If I use a simulated data, that is perfectly normally distributed, can I use this average coverage to calculate a specific point coverage?

      Many thanks

      Comment

      • Brian Bushnell
        Super Moderator
        • Jan 2014
        • 2709

        #4
        Certainly, if you use synthetic data with a normal distribution, it's possible to estimate the coverage at any point along the amplicon from the average.

        Comment

        • Fad2012
          Member
          • Sep 2012
          • 62

          #5
          Thanks very much Brian, could you please tell me a key though of how this can be done? Do I have to find first the actual coverage for one point, and use it as a reference? Or I can use the coverage average along with the mean and standard deviation? In both ways, I am not sure how to carry out this process...Confused

          Could please help me further.

          Much appreciated

          Comment

          • Brian Bushnell
            Super Moderator
            • Jan 2014
            • 2709

            #6
            Well, it's possible, once you know the constants, using the formula:



            Using 300bp reads instead of 1bp reads would make it harder; you'd need to modify the formula somehow. But it's an approximation either way, and not one that I think would be useful at all in the real world, so I don't see much point in putting in any effort to solve it.

            Comment

            • Fad2012
              Member
              • Sep 2012
              • 62

              #7
              Many thanks

              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...