broadinstitute/mutect

Try to find the strand_bias_counts definition?

ipstone opened this issue · 0 comments

Hello,

We would like to find the stand_bias information in mutect output. Currently we are using mutect 1.1.7 version, and there is a stand_bias_count field in the output.

An example of stand_bias_counts is like:
(117,114,7,7)

We assumed that it might be:

  • the reference allele on the forward strand
  • the reference allele on the reverse strand
  • the alternate allele on the forward strand
  • the alternate allele on the reverse strand

But the numbers don't add up to read depth field.
It would be very helpful if someone can tell the definition of this field, or where we can find the strand information for the mutant allele call etc. Thanks

Searching the source code, we found something might be relevant:
candidate.setStrandContingencyTable(SequenceUtils.getStrandContingencyTable(forwardPileup, reversePileup, (byte) upRef, (byte) altAllele));
but what is the deinfition of the forwardPileup, reversePileup, (byte) upRef, (byte) altAllele ?

Really appreciate any help and suggestions!