显示所有变体

Question

我有一个 2GB 的 vcf DNA 文件，我正在尝试使用 vcf_to_zarr() 打印出所有固定字段的变体，但我收到错误 KeyError: 'variants/*'

allel.vcf_to_zarr

import allel
import numcodecs
import zarr

def readVcf():

    allel.vcf_to_zarr('actual.vcf', 'example.zarr', fields='*', overwrite=True)
    callset = zarr.open_group('example.zarr', mode='r')
    allfield=callset['variants/*']

    for a in allfield:
         print(a)

Answer 1

要遍历所有变体字段，请执行：

for a in callset['variants']:
    print(a)

Zarr 不理解层次结构路径中的通配符 ('*')。