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
            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
          • seqadmin
            Investigating the Gut Microbiome Through Diet and Spatial Biology
            by seqadmin




            The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
            02-24-2025, 06:31 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 03-03-2025, 01:15 PM
          0 responses
          177 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-28-2025, 12:58 PM
          0 responses
          267 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-24-2025, 02:48 PM
          0 responses
          652 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-21-2025, 02:46 PM
          0 responses
          267 views
          0 likes
          Last Post seqadmin  
          Working...
          X