How many reads will be removed after markduplicates in general?
Hello! Now I'm dealing with bulk DNA-seq data.
Firstly, I aligned it to human genome by bwa-mem.
Before MarkDuplicates by picard with only required options(I, O, M), its mean depth of total region is 34.
However, after Markduplicates, its mean depth of total region is only 18! Moreover, after MarkDuplicates, its bam file is only 1G smaller. Why has the depth descreased so much? The depth is calculated by mosdepth (mosdepth -t 4 -n).
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
PopGenome lost individual information
written by Yixuan •
-
Bowtie2 Alignment with low percentage of properly paired reads (ATAC-seq)
written by zau saa •
-
MarkDuplicates gives too few reads
written by wormball •
-
Understanding MarkDuplicates and featureCounts results
written by xosugarxspicexo •
-
Two bams with similar flagstat output have very different file size
written by freeheart.steven •
-
Peak calling with MACS for paired end DNase-seq data
written by Bioinformatist Newbie
-
mirDeep2 using bowtie vs. bwa - why do more aligned reads yield less miRNA
written by Michael Dondrup
-
True or False, SNVs using bwa+gatk
written by J.F.Jiang
-
True or False, SNVs using bwa+gatk
written by J.F.Jiang
-
Error On Executing "Local Realign" Step Of Call For The Input Bam File
written by febz.babz •
run
samtools flagstatsto count the number of marked/total reads.Did you simply mark or actually removed the duplicates? BAM file could have shrunk if reads simply got rearranged.
or the compression level has changed ...
I realize that I only identified the duplicates.
If I forget to remove duplicates, will it influences the SNPs called by gatk HaplotypeCaller?
No. https://gatk.broadinstitute.org/hc/en-us/articles/360037225632-HaplotypeCaller
Thanks for your answer!
Returning to the original question， would you think it's strange that the depth of bulk DNA-seq data decreases too much after markduplicates?