zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Merge pull request #123 from broadinstitute/rp_event_mapper_bug_GSA-877 

Bug fix in HC gga mode.
This commit is contained in:
Mark DePristo 2013-03-21 11:42:20 -07:00
parent 1a95ce5dcf b9c331c2fa
commit b30deda7ee
1 changed files with 4 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
protected/java/src/org/broadinstitute/sting/gatk/walkers/haplotypecaller/GenotypingEngine.java
View File

@ -631,6 +631,10 @@ public class GenotypingEngine {
if( eventToTest.getKey().equals(new Event(null)) ) if( eventToTest.getKey().equals(new Event(null)) )
continue; continue;
// only try to disambiguate for alleles that have had haplotypes previously assigned above
if( eventToTest.getValue().isEmpty() )
continue;
final Haplotype artificialHaplotype = eventToTest.getValue().get(0); final Haplotype artificialHaplotype = eventToTest.getValue().get(0);
if( isSubSetOf(artificialHaplotype.getEventMap(), h.getEventMap(), true) ) { if( isSubSetOf(artificialHaplotype.getEventMap(), h.getEventMap(), true) ) {
matchingEvent = eventToTest.getKey(); matchingEvent = eventToTest.getKey();