Skip to content

Commit bf37590

Browse files
committed
Bugfig to lofreq merge: set stop for indels
1 parent 4ca1268 commit bf37590

File tree

1 file changed

+5
-2
lines changed

1 file changed

+5
-2
lines changed

SequenceAnalysis/src/org/labkey/sequenceanalysis/run/analysis/MergeLoFreqVcfHandler.java

Lines changed: 5 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -691,8 +691,11 @@ private void runSnpEff(JobContext ctx, Map<String, SiteAndAlleles> siteToAllele,
691691

692692
VariantContextBuilder b = new VariantContextBuilder();
693693
b.chr(tokens[0]);
694-
b.start(Integer.parseInt(tokens[1]));
695-
b.alleles(Arrays.asList(Allele.create(tokens[2], true), Allele.create(tokens[3])));
694+
int start = Integer.parseInt(tokens[1]);
695+
b.start(start);
696+
List<Allele> alleles = Arrays.asList(Allele.create(tokens[2], true), Allele.create(tokens[3]));
697+
b.alleles(alleles);
698+
b.computeEndFromAlleles(alleles, start);
696699
snpEffSites.add(b.make());
697700
}
698701

0 commit comments

Comments
 (0)