This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Does Markduplicates Work If Mate Pair Information Is Not 'Fixed'?

I would like to use the picard's MarkDuplicates to count the number of duplicates in a bam file. However, I have not yet run samtools fixmate to ensure that all mate-pair information is correct. Will this effect the output from MarkDuplicates? Does MarkDuplicates work even if mate pair information is not 'fixed'?

The bam file was produced with bwa and it has been sorted use picard SortSam.

bam duplicates markduplicates picard
• 4,135 views • link
updated by Biostar Biostar • written by ian.streeter • ian.streeter

1 answer

You should be able to run MarkDuplicates if you use the VALIDATION_STRINGENCY=LENIENT option.

• 1 views • link
written by Johan Johan

Log in to answer this question.

More posts like this
  • How To Sort Sam File
    written by Chen Sun Chen Sun
    <p>I want to process the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file created by <a href='http://bio-bwa.sourceforge.net/'>BWA</a>, if it is sorted by the sequence coordinate on reference, it will be …
  • Sam To Bam - Loss Of Data Or Just Great Compression?
    written by Clare Clare
    <p>I use picard tools to take a <a href='http://bio-bwa.sourceforge.net/'>BWA</a> alignment <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and convert it into a sorted <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. Normally this works …
  • Remove Duplicates Workflow
    written by Anjan Anjan
    <p>Hello all, Am using <a href='http://samtools.sourceforge.net/'>samtools</a> to remove duplicates. Here is my workflow for removing duplicates from a *.bam file. 1. Sort <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file …
  • Removing Phix Data From Ngs Data For Analysis
    written by crivenster crivenster
    <p>I have a NGS data containing 95-98% of PHIX data,which i dont need. So after some reading, I used <a href='http://bio-bwa.sourceforge.net/'>BWA</a> to align the reads …
  • Questions About Pair-End Reads Which Are Not Match For Bwa
    written by Tonyzeng • Tonyzeng
    <p>My question is that my original data or pair end read data are consisted of Read1 and Read2 and they are corresponding to each other. …
  • Can You Make A Bam File Picard Compatible? I Forgot -M With Bwa Mem
    written by ajc8 • ajc8
    <p>Hi,</p> <p>I just aligned several exome samples with <a href='http://bio-bwa.sourceforge.net/'>bwa</a> mem and then processed the output through several pipeline steps. This is the first time …
  • Problem Running Pindel With Bam Files
    written by xenophiliuslovegood • xenophiliuslovegood
    <p>Hi, I am running <strong>Pindel version 0.2.4p, March 26 2012</strong>, on a set of <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files. It chokes on them, printing out:</p> <p><strong>No currentState.Reads …
  • Markduplicates Error
    written by Leandro Batista • Leandro Batista
    <p>Hello,</p> <p>I'm processing whole genome <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files. Since I'm specifically interested on chromosome 11, I have split my files and I'm working only on …
  • Soap Result Converted To Bam Not Working In Samtool
    written by Sarwar • Sarwar
    <p>I converted the output of soap2 which is map file to <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file using map2sam.pl and then converted this <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> to <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> …
  • Samtools + Picard Markduplicates
    written by Nery • Nery
    <p>Hello,</p> <p>I'm using <a href='http://samtools.sourceforge.net/'>Samtools</a> to call variants and I am using Picard MarkDuplicates to mark duplicates in my <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. My question is: …