try filtering

bpblanken · bpblanken · commit 62125650aea9 · 2024-11-30T20:10:27.000-05:00
diff --git a/v03_pipeline/lib/reference_datasets/topmed.py b/v03_pipeline/lib/reference_datasets/topmed.py
@@ -15,6 +15,7 @@
 
 def get_ht(path: str, reference_genome: ReferenceGenome) -> hl.Table:
     ht = vcf_to_ht(path, reference_genome)
+    ht = ht.filter(ht.locus.position == hl.int(ht.info.OriginalStart))
     return ht.select(
         **{k: parse_nested_field(ht, v) for k, v in SELECT.items()},
     )

Original file line number	Diff line number	Diff line change
`@@ -15,6 +15,7 @@`
`15`	`15`
`16`	`16`	`def get_ht(path: str, reference_genome: ReferenceGenome) -> hl.Table:`
`17`	`17`	`ht = vcf_to_ht(path, reference_genome)`
	`18`	`+ ht = ht.filter(ht.locus.position == hl.int(ht.info.OriginalStart))`
`18`	`19`	`return ht.select(`
`19`	`20`	`**{k: parse_nested_field(ht, v) for k, v in SELECT.items()},`
`20`	`21`	`)`