Question: Fisher Strand in GATK, a good measure to detect Strand Bias?
gravatar for cristianrohr768
4.9 years ago by
cristianrohr76830 wrote:


I follow the GATK good practices to call variants on 11 samples from a custom truseq design. I do Joint genotyping and apply Hard Filters to my data

In a second test, i skip the joint genotyping step, and just apply hard filters to my individuals VCFs.

With boot strategies i found the same variant in the same individual wich is a pathogenic variant for a disease i was interested in.

For this variant i have a read depth of 180, and i see A:103 and G:77, G is the reference, there aren't homopolymers near. But the problem is that this position is ONLY covered by forward reads.

The FS calculated by GATK for this position is FS=0.000

FS is Phred-scaled p-value using Fisher’s Exact Test to detect strand bias (the variation being seen on only the forward or only the reverse strand) in the reads. More bias is indicative of false positive calls.

So a value of 0 must be indicative of no bias, but this position is only covered by forward reads, so i undoubtedly had strand bias in this position.

So what's going on?? i can have confidence on this variant?





fisher strand strand bias gatk • 3.5k views
ADD COMMENTlink modified 12 months ago by Pablo Marin-Garcia1.8k • written 4.9 years ago by cristianrohr76830
gravatar for Zaag
4.9 years ago by
Zaag800 wrote:

It tests the difference between variant and ref alleles so probably the ref alleles on that position are also on the forward strand, so the bias is the same for ref and var reads so there is no reason there to discard the variant.

ADD COMMENTlink written 4.9 years ago by Zaag800
gravatar for Pablo Marin-Garcia
12 months ago by
Pablo Marin-Garcia1.8k wrote:

One nice explanation could be found in this samtools mailing list post

If you have only reads in one strand you can not "have" strand bias, or at least you can not measure it, since the other strand has 0 reads.

The strand bias measure if the proportion of alleles is statistically different between strands. As you don't have reads in the other strand, you can not measure it. You can only measure allele bias (departure of the 0.5 if you are expecting germline heterozygotes)

ADD COMMENTlink written 12 months ago by Pablo Marin-Garcia1.8k
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: 706 users visited in the last hour