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
            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
          • seqadmin
            Strategies for Sequencing Challenging Samples
            by seqadmin


            Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
            03-22-2024, 06:39 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 04-11-2024, 12:08 PM
          0 responses
          29 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 10:19 PM
          0 responses
          31 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 09:21 AM
          0 responses
          28 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-04-2024, 09:00 AM
          0 responses
          52 views
          0 likes
          Last Post seqadmin  
          Working...
          X