Question: HaplotypeCaller filter out variants covered only in one direction
gravatar for stanedav
2.3 years ago by
stanedav40 wrote:


I would like to ask you, if anyone has experience with filtering out variants, that are not covered in both directions by reads (WES and Gene panels). I am using GATK workflow with haplotype caller and hard filtering, but after all steps it remains even the variants, which are covered only 2+ 0- or 3+ 0- and I want these variants with zero in any direction not to be called.

Any ideas? Thank you very much

Here is syntax of mine HaplotypeCaller command:

java -jar $gatk -T HaplotypeCaller -R hg19.fasta -L bedfile -I sample.bam -o output.g.vcf -nct 16 --emitRefConfidence GVCF --variant_index_type LINEAR --variant_index_parameter 128000 --dontUseSoftClippedBases

ADD COMMENTlink modified 2.3 years ago by pfs270 • written 2.3 years ago by stanedav40
gravatar for pfs
2.3 years ago by
pfs270 wrote:

This link may help: . This option will probably require post processing of the VCF output file. I would recommend going back in your pipeline and removing any steps that removes duplicates. Then calculating a strand bias ratio. Then using this ratio value as a filter.

ADD COMMENTlink written 2.3 years ago by pfs270

You do want to calculate StrandBias, unless you are doing amplicon sequencing though you absolutely do not want to remove duplicate filtering. That is just adding noise and bias into your data.

ADD REPLYlink written 2.3 years ago by Dan Gaston7.1k

Good afternoon Dan, It is unclear to me why you would not calculate a strand bias ratio for amplicon sequencing. There should be reads associated with both the forward and reverse primers. With respect to removing duplicates, if you have computational/time constraints I understand removing duplicates otherwise I would argue you are introducing bias and removing signal by removing duplicates. Given the low read depth the poster gave why would you remove signal? How would you calculate an accurate strand bias ratios if you are removing duplicates?

ADD REPLYlink written 2.3 years ago by pfs270

I think you misinterpreted my meaning. In all cases the OP wants to use StrandBias as a filter. However, if you are working with any form of hybridization data you absolutely should be doing duplicate removal. If you are working with amplicon data you should skip doing duplicate removal because you will be throwing out legitimate data. The original poster stated they are working with both WES and Gene panel data, the gene panel data could be coming from amplicon or hybrid capture approaches. If they are not doing amplicon sequencing they should never skip PCR duplicate removal.

ADD REPLYlink written 2.3 years ago by Dan Gaston7.1k
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: 1227 users visited in the last hour