Question: HaplotypeCaller filter out variants covered only in one direction
gravatar for stanedav
20 months ago by
stanedav30 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 20 months ago by pfs260 • written 20 months ago by stanedav30
gravatar for pfs
20 months ago by
pfs260 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 20 months ago by pfs260

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 20 months 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 20 months ago by pfs260

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 20 months 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: 1514 users visited in the last hour