Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • GC content of duplicate reads only

    I have a duplicate level of approx 10% in a DNA resequencing experiment and my read depth drops off in a linear fashion towards approx 70% GC.
    The assay is a SureSelect, done on an Illumina analyser.
    I'm fairly confident that the poor coverage in highGC is due to some sort of PCR amplification bias, but was interested in looking specifically at the duplicate reads re their GC content.
    How do I go about splitting the duplicates out of my bwa-aligned bam into their own bam so that I can look into this specifically?
    Thanks, Chris

  • #2
    if you use picard MarkDuplicates (REMOVE_DUPLICATES=FALSE) you can filter using samtools view -b -f 1024 to obtain only reads which are marked as duplicates in the flag field.

    Comment


    • #3
      I think at the moment my pipeline is converting sam to bam, then marking duplicates. Presumably I can split the dups out of a bam file rather than having to convert back to sam, splitting out the dups then converting to bam for all the other downstream analyses?
      Chris

      Comment


      • #4
        volks suggestion works on bam (hence the -b option).

        Comment


        • #5
          yes, from bam to bam (-b).
          there is hardly ever a reason to write a sam file to disk.

          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...
            04-22-2024, 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, Yesterday, 11:49 AM
          0 responses
          15 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-24-2024, 08:47 AM
          0 responses
          16 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-11-2024, 12:08 PM
          0 responses
          61 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 10:19 PM
          0 responses
          60 views
          0 likes
          Last Post seqadmin  
          Working...
          X