diff --git a/java/src/org/broadinstitute/sting/gatk/walkers/indels/IndelGenotyperV2Walker.java b/java/src/org/broadinstitute/sting/gatk/walkers/indels/IndelGenotyperV2Walker.java index ecee83bb5..5929eea38 100644 --- a/java/src/org/broadinstitute/sting/gatk/walkers/indels/IndelGenotyperV2Walker.java +++ b/java/src/org/broadinstitute/sting/gatk/walkers/indels/IndelGenotyperV2Walker.java @@ -624,6 +624,7 @@ public class IndelGenotyperV2Walker extends ReadWalker { //if ( pos == 3534096 ) System.out.println("pos="+pos +" total reads: "+context.getReads().size()); Iterator read_iter = context.getReads().iterator(); + while ( read_iter.hasNext() ) { ExpandedSAMRecord rec = read_iter.next(); SAMRecord read = rec.getSAMRecord(); @@ -631,6 +632,7 @@ public class IndelGenotyperV2Walker extends ReadWalker { byte[] quals = rec.getExpandedQuals(); int mm = rec.getMMCount(); + if( read.getAlignmentStart() > pos || read.getAlignmentEnd() < pos ) continue; long local_right = right; // end of nqs window for this particular read. May need to be advanced further right @@ -638,9 +640,9 @@ public class IndelGenotyperV2Walker extends ReadWalker { // automatically since flags/quals are set to -1 there boolean read_has_a_variant = false; - boolean read_has_consensus = ( consensus_indel!= null && consensus_indel.getReadSet().contains(read) ); + boolean read_has_consensus = ( consensus_indel!= null && consensus_indel.getReadSet().contains(rec) ); for ( IndelVariant v : variants ) { - if ( v.getReadSet().contains(read) ) { + if ( v.getReadSet().contains(rec) ) { read_has_a_variant = true; local_right += v.lengthOnRef(); break;