From 3dabb90eb0502bcce352d3a34e99e0782e683d20 Mon Sep 17 00:00:00 2001 From: Ryan Poplin Date: Mon, 30 Jul 2012 21:26:16 -0400 Subject: [PATCH] Updating example active region walker integration test. --- .../sting/gatk/walkers/haplotypecaller/GenotypingEngine.java | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/protected/java/src/org/broadinstitute/sting/gatk/walkers/haplotypecaller/GenotypingEngine.java b/protected/java/src/org/broadinstitute/sting/gatk/walkers/haplotypecaller/GenotypingEngine.java index 5fc466336..2787689b5 100644 --- a/protected/java/src/org/broadinstitute/sting/gatk/walkers/haplotypecaller/GenotypingEngine.java +++ b/protected/java/src/org/broadinstitute/sting/gatk/walkers/haplotypecaller/GenotypingEngine.java @@ -379,7 +379,6 @@ public class GenotypingEngine { } if( thisVC == null || nextVC == null ) { continue; - //throw new ReviewedStingException("StartPos TreeSet has an entry for an event that is found on no haplotype. start pos = " + thisStart + ", next pos = " + nextStart); } if( isBiallelic ) { final double R2 = calculateR2LD( Math.pow(10.0, x11), Math.pow(10.0, x12), Math.pow(10.0, x21), Math.pow(10.0, x22) ); @@ -441,7 +440,7 @@ public class GenotypingEngine { altBases = ArrayUtils.addAll(altBases, nextVC.getAlternateAllele(0).getBases()); int iii = 0; - if( refBases.length == altBases.length ) { // special case of insertion + deletion of same length creates an MNP --> trim padding bases off the allele + if( refBases.length == altBases.length ) { // insertion + deletion of same length creates an MNP --> trim common prefix bases off the beginning of the allele while( iii < refBases.length && refBases[iii] == altBases[iii] ) { iii++; } } final ArrayList mergedAlleles = new ArrayList();