Thank you very much for the suggestion, shi. I will give it a try.
Best regards,
Marcin.
-
You may try the featureCounts function in the bioconductor package package Rsubread.Leave a comment:
-
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.
-
Spatial biology is an exciting field that encompasses a wide range of techniques and technologies aimed at mapping the organization and interactions of various biomolecules in their native environments. As this area of research progresses, new tools and methodologies are being introduced, accompanied by efforts to establish benchmarking standards and drive technological innovation.
3D Genomics
While spatial biology often involves studying proteins and RNAs in their...01-01-2025, 07:30 PM -
Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...12-16-2024, 07:57 AM
Leave a comment: