Skip to content

Commit 49fc182

Browse files
Fixup
1 parent 4accfc5 commit 49fc182

File tree

1 file changed

+2
-9
lines changed

1 file changed

+2
-9
lines changed

bio2zarr/vcf2zarr/vcz.py

Lines changed: 2 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -848,17 +848,10 @@ def encode_local_alleles_partition(self, partition_index):
848848
store=self.wip_partition_array_path(partition_index, "call_genotype"),
849849
mode="r",
850850
)
851-
alleles_array = zarr.open_array(
852-
store=self.wip_partition_array_path(partition_index, "variant_allele"),
853-
mode="r",
854-
)
855851
for chunk_index in range(gt_array.cdata_shape[0]):
856-
A = alleles_array.blocks[chunk_index]
857-
G = gt_array.blocks[chunk_index]
858-
for alleles, var in zip(A, G):
852+
for genotypes in gt_array.blocks[chunk_index]:
859853
j = call_LAA.next_buffer_row()
860-
# TODO we should probably compute LAAs by chunk for efficiency
861-
call_LAA.buff[j] = compute_laa_field(var, alleles)
854+
call_LAA.buff[j] = compute_laa_field(genotypes)
862855

863856
call_LAA.flush()
864857
self.finalise_partition_array(partition_index, "call_LAA")

0 commit comments

Comments
 (0)