Update logic used in how to track novel site when no writer is used (#…

…286)
BimberLab · Sep 26, 2023 · 361e163 · 361e163
1 parent af1902e
commit 361e163
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/main/java/com/github/discvrseq/walkers/AnnotateNovelSites.java b/src/main/java/com/github/discvrseq/walkers/AnnotateNovelSites.java
@@ -127,10 +127,10 @@ public void apply(List<VariantContext> variantContexts, ReferenceContext referen
             // Site is novel:
             if (refVariants == null || refVariants.isEmpty()) {
                 vcb.attribute(novelSiteAnnotationName, novelSiteAnnotationValue);
+                novelSites++;
                 writer.add(vcb.make());
 
                 if (novelSitesWriter != null) {
-                    novelSites++;
                     novelSitesWriter.add(vcb.genotypes().make());
                 }
             }
@@ -162,13 +162,13 @@ public void apply(List<VariantContext> variantContexts, ReferenceContext referen
                 if (!foundMatchingRef) {
                     vcb.attribute(novelSiteAnnotationName, novelSiteAnnotationValue);
                     isNovel = true;
+                    novelSites++;
                 }
 
                 vc = vcb.make();
                 writer.add(vc);
 
                 if (isNovel && novelSitesWriter != null) {
-                    novelSites++;
                     novelSitesWriter.add(vcb.genotypes().make());
                 }
             }