Created a new tool, SiblingIBD, which finds Identical-By-Descent regions in two siblings.
-When parental genotypes are available, implements an HMM on genotype observations in the quartet. -Outputs IBD regions as well as per-site posterior probabilities of being in each IBD state. -Includes an experimental heuristic based mode for when parental genotypes are not available. -Made a method in MendelianViolation public static to reuse code. -Added the mockito library to private/gatk-tools-private/pom.xml
This commit is contained in:
parent
fbc75c8e14
commit
ba1d23e535
|
|
@ -373,7 +373,15 @@ public class MendelianViolation {
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
private boolean isViolation(Genotype gMom, Genotype gDad, Genotype gChild) {
|
/**
|
||||||
|
* Evaluate the genotypes of mom, dad, and child to detect Mendelian violations
|
||||||
|
*
|
||||||
|
* @param gMom
|
||||||
|
* @param gDad
|
||||||
|
* @param gChild
|
||||||
|
* @return true if the three genotypes represent a Mendelian violation; false otherwise
|
||||||
|
*/
|
||||||
|
public static boolean isViolation(final Genotype gMom, final Genotype gDad, final Genotype gChild) {
|
||||||
//1 parent is no "call
|
//1 parent is no "call
|
||||||
if(!gMom.isCalled()){
|
if(!gMom.isCalled()){
|
||||||
return (gDad.isHomRef() && gChild.isHomVar()) || (gDad.isHomVar() && gChild.isHomRef());
|
return (gDad.isHomRef() && gChild.isHomVar()) || (gDad.isHomVar() && gChild.isHomRef());
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue