Question: HaplotypeCaller filter out variants covered only in one direction
0
gravatar for stanedav
3 months ago by
stanedav10
stanedav10 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

ADD COMMENTlink modified 3 months ago by pfs170 • written 3 months ago by stanedav10
0
gravatar for pfs
3 months ago by
pfs170
USA/Boston
pfs170 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 COMMENTlink written 3 months ago by pfs170
1

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 3 months ago by Dan Gaston6.9k

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 3 months ago by pfs170

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 3 months ago by Dan Gaston6.9k
Please log in to add an answer.

Help
Access

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