我有一个 2GB 的 vcf DNA 文件,我正在尝试使用 vcf_to_zarr() 打印出具有所有固定字段的所有变体,但我收到错误 KeyError: 'variants/*'
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)