Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Issues with the countOverlaps function

    Hi,

    I am busy with an RNA analysis, and following the steps outlined in Next Generation Sequencing Wiki on http://en.wikibooks.org/wiki/Next_Ge...cing_(NGS)/RNA, I have stumbled upon a problem I cannot solve: the counOverlaps function seems to count hits for a gene even if the gene wholy falls into an area between two reads. As an example, consider the IGV.pdf file attached below. It shows the area occupied by the PPP1R26P2 pseudogene, which is located on chr22:20,499,890-20503484. PPP1R26P2 is fully contained between two reads in case of 43 pairs, and countOverlaps(tx_by_gene, reads) yields the result 43 for this gene.
    The result should be 0, because "reads" does not contain any information about read pairs, only information about each read's location separately, and locations of PPP1R26P2 and neighbouring reads are disjoint.

    Do you have any idea how to make this function work properly or what to use instead?

    Best regards,
    Marcin.
    Attached Files

  • #2
    You may try the featureCounts function in the bioconductor package package Rsubread.

    Comment


    • #3
      Thank you very much for the suggestion, shi. I will give it a try.

      Best regards,
      Marcin.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Quality Control Essentials for Next-Generation Sequencing Workflows
        by seqadmin




        Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

        Nucleic Acid Quality Control
        Preparing for NGS starts with isolating the...
        Yesterday, 01:58 PM
      • seqadmin
        An Introduction to the Technologies Transforming Precision Medicine
        by seqadmin


        In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
        01-27-2025, 07:46 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 02-07-2025, 09:30 AM
      0 responses
      48 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-05-2025, 10:34 AM
      0 responses
      66 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-03-2025, 09:07 AM
      0 responses
      54 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 01-31-2025, 08:31 AM
      0 responses
      43 views
      0 likes
      Last Post seqadmin  
      Working...
      X