Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: MarkDuplicates output file in GATK pipeline
1
gravatar for gprashant17
11 months ago by
gprashant17 • 70
gprashant17 • 70 wrote:

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?
sequencing rna-seq bam alignment gatk • 231 views
ADD COMMENT • link •
modified 11 months ago by Pierre Lindenbaum ♦ 128k • written 11 months ago by gprashant17 • 70
1
gravatar for Pierre Lindenbaum
11 months ago by
Pierre Lindenbaum ♦ 128k
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum ♦ 128k wrote:

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.
ADD COMMENT • link written 11 months ago by Pierre Lindenbaum ♦ 128k

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?

ADD REPLY • link written 11 months ago by gprashant17 • 70

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

ADD REPLY • link written 11 months ago by Pierre Lindenbaum ♦ 128k
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 952 users visited in the last hour