Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • swNGS
    Member
    • Nov 2011
    • 83

    #1

    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
  • volks
    Member
    • Jun 2010
    • 80

    #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

    • swNGS
      Member
      • Nov 2011
      • 83

      #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

      • jflowers
        Member
        • Oct 2011
        • 42

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

        Comment

        • volks
          Member
          • Jun 2010
          • 80

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

          Comment

          Latest Articles

          Collapse

          • SEQadmin2
            Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
            by SEQadmin2



            Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

            Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
            ...
            10-14-2026, 09:55 AM
          • SEQadmin2
            New Genomics Technologies Take Aim at Long-Standing Limits
            by SEQadmin2


            Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

            We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
            ...
            09-28-2026, 10:25 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, 10-15-2026, 10:52 AM
          0 responses
          23 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-14-2026, 05:39 AM
          0 responses
          27 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-13-2026, 09:56 AM
          0 responses
          23 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-09-2026, 11:14 AM
          0 responses
          28 views
          0 reactions
          Last Post SEQadmin2  
          Working...