Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • count duplicates in bam files

    Hi,

    I would like to count duplicates in bam files. I am comparing two mapping tools and looking at the total counts. Typically, it should be the same
    value. My assumption is that BWA mem counts duplicates.

    What I did:

    1) convert sam -> bam file
    2) sort bam file
    3) use the Picard tool MarkDuplicates.jar
    4) use the BuildBamIndex.jar

    Mapped Reads (CLC) 6,876,285
    Mapped Reads (BWA) 6,375,889
    Unmapped Reads (CLC) 231,927
    Unmapped Reads (BWA) 7,367,30
    Total count (CLC) 7,108,212
    Total count (BWA) 7,112,619

    What would be the next step ? I ve tried to use the .bai files ... but do they have information about the number of duplicates ?

    Do you have any suggestions ?

    Best,
    Flo

  • #2
    If you are just looking to get the total number of duplicates you could run Qualimap on your bam files. You will get all sorts of additional info as well.

    This thread has some other suggestions: http://seqanswers.com/forums/showthread.php?t=23493

    Comment


    • #3
      Doesn't the metrics file produced by MarkDuplicates.jar already contain the information the duplicate counts?

      Comment


      • #4
        Doh! A file name is required for Metrics file so @Flo89 must have provided it when running MarkDuplicates.

        Comment


        • #5
          Hi guys,

          Thank you very much. I am not familiar with this mapping stuff. So, I guess it was an easy question for you guys.

          Best,
          Flo

          Comment

          Latest Articles

          Collapse

          • seqadmin
            The Impact of AI in Genomic Medicine
            by seqadmin



            Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
            02-26-2024, 02:07 PM
          • seqadmin
            Multiomics Techniques Advancing Disease Research
            by seqadmin


            New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

            A major leap in the field has
            ...
            02-08-2024, 06:33 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, Yesterday, 06:12 AM
          0 responses
          20 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-23-2024, 04:11 PM
          0 responses
          67 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-21-2024, 08:52 AM
          0 responses
          75 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-20-2024, 08:57 AM
          0 responses
          66 views
          0 likes
          Last Post seqadmin  
          Working...
          X