diff --git a/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/SNPGenotypeLikelihoodsCalculationModel.java b/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/SNPGenotypeLikelihoodsCalculationModel.java index f57c8873e..54ce35bbd 100755 --- a/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/SNPGenotypeLikelihoodsCalculationModel.java +++ b/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/SNPGenotypeLikelihoodsCalculationModel.java @@ -81,10 +81,11 @@ public class SNPGenotypeLikelihoodsCalculationModel extends GenotypeLikelihoodsC return null; } if ( !vcInput.isBiallelic() ) { - logger.info("Record at position " + ref.getLocus() + " is not bi-allelic; choosing the first allele..."); - //return null; + // for multi-allelic sites go back to the reads and find the most likely alternate allele + initializeBestAlternateAllele(refBase, contexts); + } else { + bestAlternateAllele = vcInput.getAlternateAllele(0).getBases()[0]; } - bestAlternateAllele = vcInput.getAlternateAllele(0).getBases()[0]; } else { initializeBestAlternateAllele(refBase, contexts); }