Don't try to phase triple-hets either.

This commit is contained in:
Kiran V Garimella 2011-07-23 01:38:58 -04:00
parent 45f2ca8d99
commit 57e3d136eb
1 changed files with 1 additions and 1 deletions

View File

@ -234,7 +234,7 @@ public class PhaseByTransmission extends RodWalker<Integer, Integer> {
finalGenotypes.add(father);
finalGenotypes.add(child);
if (mother.isCalled() && father.isCalled() && child.isCalled()) {
if (mother.isCalled() && father.isCalled() && child.isCalled() && !(mother.isHet() && father.isHet() && child.isHet())) {
ArrayList<Genotype> possibleMotherGenotypes = createAllThreeGenotypes(ref, alt, mother);
ArrayList<Genotype> possibleFatherGenotypes = createAllThreeGenotypes(ref, alt, father);
ArrayList<Genotype> possibleChildGenotypes = createAllThreeGenotypes(ref, alt, child);