Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Try to find the strand_bias_counts definition? #12

Open
ipstone opened this issue Jun 26, 2018 · 0 comments
Open

Try to find the strand_bias_counts definition? #12

ipstone opened this issue Jun 26, 2018 · 0 comments

Comments

@ipstone
Copy link

ipstone commented Jun 26, 2018

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!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant