From 97caf94807a0bf5ba9d3c35e5d9069c7922c0b22 Mon Sep 17 00:00:00 2001 From: Laura Gauthier Date: Thu, 23 Apr 2015 13:01:47 -0400 Subject: [PATCH] Fix implementation of allowNonUniqueKmersInRef so that it applies to all kmer sizes --- .../haplotypecaller/readthreading/ReadThreadingAssembler.java | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/haplotypecaller/readthreading/ReadThreadingAssembler.java b/protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/haplotypecaller/readthreading/ReadThreadingAssembler.java index fa54a4fda..f9de3173a 100644 --- a/protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/haplotypecaller/readthreading/ReadThreadingAssembler.java +++ b/protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/haplotypecaller/readthreading/ReadThreadingAssembler.java @@ -124,7 +124,7 @@ public class ReadThreadingAssembler extends LocalAssemblyEngine { while ( results.isEmpty() && numIterations <= MAX_KMER_ITERATIONS_TO_ATTEMPT ) { // on the last attempt we will allow low complexity graphs final boolean lastAttempt = numIterations == MAX_KMER_ITERATIONS_TO_ATTEMPT; - addResult(results, createGraph(reads, refHaplotype, kmerSize, givenHaplotypes, lastAttempt, lastAttempt)); + addResult(results, createGraph(reads, refHaplotype, kmerSize, givenHaplotypes, lastAttempt, allowNonUniqueKmersInRef || lastAttempt)); kmerSize += KMER_SIZE_ITERATION_INCREASE; numIterations++; }