This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Variant calling step order question: base recalibration & mark duplicates, which is first?

Hi all,

We're going through & revising our variant calling pipeline on NGS data from cancer patients and a question came up:

Which step should be done first (and why), base recalibration or mark duplicates?

Currently we recalibrate bases first and then mark duplicates.

The reason I'm asking this is that we originally based part of our pipeline on the following article, which said that you recalibrate bases and then mark duplicates: http://www.htslib.org/workflow/#mapping_to_variant

However, in the following Broad Institute best practices page it says the opposite, you mark duplicates and then recalibrate bases, saw it in another paper as well: https://software.broadinstitute.org/gatk/best-practices/bp_3step.php?case=GermShortWGS

Thanks in advance!

Alon

ngs variant calling cancer pipeline
• 3,452 views • link
updated by Brian Bushnell Brian Bushnell • written by alons alons

2 answers

I'd probably remove duplicates first, since BSRC is generating some sort of covariation model with all of the supplied reads. I'm assuming that having a bunch of clonal artifacts in your dataset might throw this off a little. But honestly, you should ask the GATK people as they have a better understanding of the underlying model.

• 0 views • link
written by mforde84 mforde84

Recalibrating bases should not really improve (or affect) duplicate detection. But duplicate removal can improve recalibration, so I'd do that first. And the earlier you remove duplicates, the faster everything else becomes.

• 0 views • link
written by Brian Bushnell Brian Bushnell

Log in to answer this question.

More posts like this