Holiday Notice:
The Frontline Support team will be slow to respond December 17-18 due to an institute-wide retreat and offline December 22- January 1, while the institute is closed. Thank you for your patience during these next few weeks. Happy Holidays!

High-depth filter for Mutect calls

akrishnaakrishna Member
edited December 2014 in Ask the GATK team

Hi,

We are working on a cancer genome project and we use Mutect to call somatic mutations.
We have a question related to the high-coverage depth filtering. Below is a list of nearby positions output by Mutect
where the coverage in the tumor sample is very high but the coverage in the normal sample is not.

Partial output from *call_stats.out file -

contig position t_q20_count n_q20_count failure_reasons judgement

2 89849000 622 108 KEEP

2 89850117 508 152 KEEP

2 89850498 713 732 KEEP

2 89850649 583 403 KEEP

2 89850993 849 142 KEEP

2 89872002 540 286 KEEP

2 89877959 607 259 KEEP

These positions are not filtered out by Mutect.
I wonder why is a high-depth filter not used?
Would it make sense to filter these variants and which high-depth threshold should we choose?

Regards,

Abhimanyu Krishna

Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    Hi there,

    Unfortunately, it is very difficult to give one-size-fits-all recommendations, and at this time we cannot give you detailed guidelines for further filtering. You will need to determine which filtering parameters are appropriate for your dataset. Good luck.

Sign In or Register to comment.