Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bug in featureCounts?

    We are using featureCounts for both and single-end and pair-end data. In the command we always have the '-p' flag on, assuming it will be ignored for single-end data. However, we notice that the function of '-M' is affected because of '-p' flag.

    When we run single-end data with '-M' flag, all the reads with NH > 1 are ignored, as expected. However, when we run single-end data with both '-M' and '-p' flag, the multi-mapping reads (NH > 1) are no longer ignored.

    The '-p' should have an effect on filtering multi-mapping reads. I'm attaching a sample sam and gtf file where this may be observed.
    Attached Files

  • #2
    Originally posted by brinda View Post
    We are using featureCounts for both and single-end and pair-end data. In the command we always have the '-p' flag on, assuming it will be ignored for single-end data. However, we notice that the function of '-M' is affected because of '-p' flag.

    When we run single-end data with '-M' flag, all the reads with NH > 1 are ignored, as expected. However, when we run single-end data with both '-M' and '-p' flag, the multi-mapping reads (NH > 1) are no longer ignored.

    The '-p' should have an effect on filtering multi-mapping reads. I'm attaching a sample sam and gtf file where this may be observed.
    As stated in the program usage and users guide, '-p' should only be used for paired-end reads. So don't use it for single end reads.

    When '-M' is specified, reads with NH>1 will NOT be ignored. They will be counted. Again, please read the program usage and users guide.

    Wei

    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