HaplotypeCaller filter out variants covered only in one direction

Question: HaplotypeCaller filter out variants covered only in one direction

12 months ago by

stanedav • 20 wrote:

Hello,

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

haplotype caller gatk variant calling • 718 views

ADD COMMENT • link •

modified 12 months ago by pfs • 250 • written 12 months ago by stanedav • 20

12 months ago by

pfs • 250

USA/Boston

pfs • 250 wrote:

This link may help: https://gatkforums.broadinstitute.org/gatk/discussion/4939/haplotypecaller-strandbiasbysample-annotation . 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 COMMENT • link written 12 months ago by pfs • 250

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 REPLY • link written 12 months ago by Dan Gaston ♦ 7.0k

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 REPLY • link written 12 months ago by pfs • 250

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 REPLY • link written 12 months ago by Dan Gaston ♦ 7.0k

Similar posts • Search »