Показать все варианты
У меня есть файл ДНК vcf объемом 2 ГБ, и я пытаюсь использовать vcf_to_zarr(), чтобы распечатать все варианты со всеми фиксированными полями, но я получаю сообщение об ошибке KeyError: 'options / *'
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)
1 ответ
Решение
Чтобы перебрать все поля вариантов, выполните:
for a in callset['variants']:
print(a)
Zarr не понимает подстановочные знаки ('*') в путях иерархии.