Update logic for VcfToLuceneIndexer and fractionHet to include just s…

…amples with the allele
BimberLab · Dec 14, 2023 · dbd365e · dbd365e
1 parent d749224
commit dbd365e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/main/java/com/github/discvrseq/walkers/VcfToLuceneIndexer.java b/src/main/java/com/github/discvrseq/walkers/VcfToLuceneIndexer.java
@@ -264,7 +264,7 @@ else if (line.getCountType() == VCFHeaderLineCount.INTEGER || line.getCountType(
                     doc.add(new IntPoint("nCalled", (int)nCalled));
                     doc.add(new StoredField("nCalled", (int)nCalled));
 
-                    float fractionHet = (float) nHet / (float) nCalled;
+                    float fractionHet = (float) nHet / (float) (nHet + nHomVar);
                     doc.add(new FloatPoint("fractionHet", fractionHet));
                     doc.add(new StoredField("fractionHet", fractionHet));