Показать все варианты

У меня есть файл ДНК vcf объемом 2 ГБ, и я пытаюсь использовать vcf_to_zarr(), чтобы распечатать все варианты со всеми фиксированными полями, но я получаю сообщение об ошибке KeyError: 'options / *'

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)

1 ответ

Решение

Чтобы перебрать все поля вариантов, выполните:

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

Zarr не понимает подстановочные знаки ('*') в путях иерархии.

Другие вопросы по тегам