From 510b3efcc203a6e1788cdae3b8e4eff008d7d859 Mon Sep 17 00:00:00 2001
From: kiran <kiran@348d0f76-0448-11de-a6fe-93d51630548a>
Date: Mon, 3 May 2010 18:46:33 +0000
Subject: [PATCH] Fixed an issue where asking for the alternate alleles at
 hom-ref sites would result in an array out-of-bounds exception.

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@3292 348d0f76-0448-11de-a6fe-93d51630548a
---
 .../gatk/walkers/vcftools/VCFSubsetWalker.java         | 10 ++++++----
 1 file changed, 6 insertions(+), 4 deletions(-)

diff --git a/java/src/org/broadinstitute/sting/playground/gatk/walkers/vcftools/VCFSubsetWalker.java b/java/src/org/broadinstitute/sting/playground/gatk/walkers/vcftools/VCFSubsetWalker.java
index 55ceb8cd6..0dff1a011 100755
--- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/vcftools/VCFSubsetWalker.java
+++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/vcftools/VCFSubsetWalker.java
@@ -147,13 +147,15 @@ public class VCFSubsetWalker extends RodWalker<ArrayList<VCFRecord>, VCFWriter>
             VCFRecord subset = subsetRecord(record);
 
             boolean isVariant = false;
-            for ( VCFGenotypeEncoding ge : subset.getVCFGenotypeRecords().get(0).getAlleles() ) {
-                if (!record.getReference().equals(ge.getBases())) {
-                    isVariant = true;
+
+            if (subset.getVCFGenotypeRecords().size() > 0) {
+                for ( VCFGenotypeEncoding ge : subset.getVCFGenotypeRecords().get(0).getAlleles() ) {
+                    if (!record.getReference().equals(ge.getBases())) {
+                        isVariant = true;
+                    }
                 }
             }
 
-            //if (isVariant && !subset.isFiltered()) {
             if ((isVariant || INCLUDE_NON_VARIANTS) && (!subset.isFiltered() || INCLUDE_FILTERED)) {
                 if (vwriter != null) {
                     vwriter.addRecord(subset);