Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • KlausNZ
    Junior Member
    • Feb 2010
    • 2

    #1

    SAMtools rmdup - saving the duplicates?

    Is there a way to save the duplicates (write to a separate file) that are removed by SAMtools' rmdup? I'd like to have a look at them....

    Many thanks in advance,
    K
  • volks
    Member
    • Jun 2010
    • 80

    #2
    try picard MarkDuplicates and then using samtools view -f 1024

    Comment

    • KlausNZ
      Junior Member
      • Feb 2010
      • 2

      #3
      Cheers, that makes sense. I wish there was a way of doing this in a single step with samtools (I'm not so fond of Picard at the moment... )
      K

      Comment

      • bruce01
        Senior Member
        • Mar 2011
        • 160

        #4
        Is it too obtuse to say that you could grep out the non-matching lines from the rmdup SAM and the original SAM?

        Also whats wrong with Picard? I've always found it more consistent.

        Comment

        • robinpsmith
          Junior Member
          • Apr 2010
          • 5

          #5
          Picard MarkDuplicates doesn't work with SE reads, as far as I can tell.

          We're estimating duplicates in ChIP-seq data.

          Comment

          • volks
            Member
            • Jun 2010
            • 80

            #6
            Originally posted by robinpsmith View Post
            Picard MarkDuplicates doesn't work with SE reads, as far as I can tell.
            i am not sure of this. why wouldnt it?

            anyway, depending on your alignment file size you could also consider running a unix diff on the original and samtools-rmdup'ed SAM file.

            Comment

            • naragam
              Member
              • Apr 2012
              • 21

              #7
              Hi,

              Am using samtools rmdup ... but it doesn't seem to remove them in the out file? What's going on?

              Nash

              Comment

              • robinpsmith
                Junior Member
                • Apr 2010
                • 5

                #8
                @volks you are right, it does work on SE. My bad.

                @naragam. I ended up using Picard, it does a much better job.

                Comment

                Latest Articles

                Collapse

                • SEQadmin2
                  Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
                  by SEQadmin2



                  CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

                  Despite this, “CRISPR helped turn genome editing from a specialized technique into
                  ...
                  07-31-2026, 11:01 AM
                • SEQadmin2
                  Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
                  by SEQadmin2


                  Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

                  The systematic characterization of the human proteome has
                  ...
                  07-20-2026, 11:48 AM
                • SEQadmin2
                  Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
                  by SEQadmin2



                  Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
                  ...
                  07-09-2026, 11:10 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by SEQadmin2, 07-31-2026, 02:55 AM
                0 responses
                18 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 07-24-2026, 12:17 PM
                0 responses
                16 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 07-23-2026, 11:41 AM
                0 responses
                15 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 07-20-2026, 11:10 AM
                0 responses
                26 views
                0 reactions
                Last Post SEQadmin2  
                Working...