This is a test version of Biostars. For the public version, visit https://www.biostars.org.
MarkDuplicates output file in GATK pipeline

I have used GATK's MarkDuplicates on a BAM file I obtained after alignment, which resulted in another file marked_duplicates.bam. So should I proceed with this marked_duplicates.bam file for analysis (converting to VCF), or this is just a file containing duplicates? In the latter case, is it possible to obtain a BAM file, with all the duplicates removed?

gatk rna-seq alignment bam sequencing
• 1,465 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by gprashant17 gprashant17

1 answer

So should I proceed with this marked_duplicates.bam file for analysis (converting to VCF),

yes. As a proof, test both files wth samtools flagstats

is it possible to obtain a BAM file, with all the duplicates removed?

in the manual : https://software.broadinstitute.org/gatk/documentation/tooldocs/4.0.4.0/picard_sam_markduplicates_MarkDuplicates.php#--REMOVE_DUPLICATES

--REMOVE_DUPLICATES / NA

If true do not write duplicates to the output file instead of writing them with appropriate flags set.
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

So if I did not use --REMOVE_DUPLICATES, the duplicate reads will still be present in the marked_duplicates.bam but they would have been flagged as duplicates right?

• 0 views • link
written by gprashant17 gprashant17

So if I did not use --REMOVE_DUPLICATES, the duplicate reads will still be present in the marked_duplicates.bam but they would have been flagged as duplicates right?

yes

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this