diff --git a/java/src/org/broadinstitute/sting/gatk/walkers/filters/VECFisherStrand.java b/java/src/org/broadinstitute/sting/gatk/walkers/filters/VECFisherStrand.java index dd0b74f8d..cbb212ac0 100644 --- a/java/src/org/broadinstitute/sting/gatk/walkers/filters/VECFisherStrand.java +++ b/java/src/org/broadinstitute/sting/gatk/walkers/filters/VECFisherStrand.java @@ -186,6 +186,10 @@ public class VECFisherStrand implements VariantExclusionCriterion { SAMRecord read = reads.get(readIndex); int offset = offsets.get(readIndex); + // skip over deletion sites + if ( offset == -1 ) + continue; + int readAllele = BaseUtils.simpleBaseToBaseIndex(read.getReadString().charAt(offset)); boolean isFW = !read.getReadNegativeStrandFlag(); diff --git a/java/src/org/broadinstitute/sting/gatk/walkers/filters/VariantFiltrationWalker.java b/java/src/org/broadinstitute/sting/gatk/walkers/filters/VariantFiltrationWalker.java index 250e65cbb..8f93b80e5 100755 --- a/java/src/org/broadinstitute/sting/gatk/walkers/filters/VariantFiltrationWalker.java +++ b/java/src/org/broadinstitute/sting/gatk/walkers/filters/VariantFiltrationWalker.java @@ -197,6 +197,14 @@ public class VariantFiltrationWalker extends LocusWalker { */ public Integer reduceInit() { return 0; } + + /** + * We want reads that span deletions + * + * @return true + */ + public boolean includeReadsWithDeletionAtLoci() { return true; } + /** * For each site of interest, rescore the genotype likelihoods by applying the specified feature set. *