Question: Lower GATK's HaplotypeCaller threshold for allele frequency as part of variant calling pipeline
gravatar for alons
3.9 years ago by
alons270 wrote:

Hi all!

As part of a variant calling pipeline for cancer i'm interested in lowering the threshold for allele frequency tolerance in GATK's HaplotypeCaller variant caller to 0.01 (1%).

Background: I know for a fact that there's at least one variant in my sample that doesn't appear in my VCF files if allele frequencies below 0.1 (10%) are filtered out during the variant calling as it's the default in some of the programs. I can see the variant when I inspect the corresponding bam file with samtools tview and Agilent's SureCall program called it itself. 

Now, i've managed to get that variant using freebayes by executing the following, essentially lowering the allele frequency threshold to 0.01:

freebayes -f ref.fa -F 0.01 -C 1 --pooled_continuous
aligned_sorted.bam > var.vcf

I want to do the same for GATK's HaplotypeCaller.  

Thanks in advance, 



ADD COMMENTlink written 3.9 years ago by alons270
gravatar for Daniel Swan
3.9 years ago by
Daniel Swan13k
Aberdeen, UK
Daniel Swan13k wrote:

Use MuTect if you're insistent on using something aligned with the GATK suite. HaploType caller isn't (AFAIK) designed with cancer calling in mind.  FreeBayes is an excellent choice for this however.

ADD COMMENTlink written 3.9 years ago by Daniel Swan13k

Daniel Swan, thank you! I will look into it. Actually i'm not very insistent on GATK, if you know other variant callers that are suitable for my purpose i'd love to hear about them.

ADD REPLYlink modified 3.9 years ago • written 3.9 years ago by alons270

You're using SureCall from Agilent - is this a SureSelect or Haloplex panel?  If it's Haloplex make sure you're not removing duplicates from the data, because this will ruin your chances of detecting a variant.

ADD REPLYlink written 3.9 years ago by Daniel Swan13k

I'm using SureSelect with removing duplicates. Actually i'm looking to use SureCall as more of a validation tool, that's why i'm interested in different variant callers.

ADD REPLYlink written 3.9 years ago by alons270
Please log in to add an answer.


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