Non-left-aligned variant in sparse_split_multi Error

Hail Query & hailctl
1

Hi,

I am trying to export only non-reference rows from the MT merging 5k gvcfs.
My pipeline is as below.

read.table
hl.experimental.densify
hl.len(densified_mt.alleles) > 1
hl.experimental.sparse_split_multi
hl.vep
MT.write

But in the MT.write() function execution, I got an error at the very last minute.

image
image2058×292 129 KB

Please let me know how to correct it or recommend another pipeline if you have.

-Jina

2

Hi Hail team,

I am looking forward to your advice eagerly. If you need any more information, please let me know. Thank you always.

-Jina

3

This is an unexpected input – Hail’s sparse matrix utilities work against left-aligned variants. What caller generated your GVCFs?

4

Hi @tpoterba,

Thank you for your response. We applied the same variant calling pipeline including Haplotypecaller (GATK) for all genomes. By the way, this error occurred not in the “hl.experimental.sparse_split_multi” step, but in the “MT.write” step. Why does it happen?

Please suggest possible solutions you think. Thank you.

-Jina