seqMerge error

[davidroberson]

I split a large bcf (bcftools view) and converted each shard to a gds. I am now trying to merge the shards back together but not all of the smaller gds will merge. Occasionally I get this error

Error in .append_node_variant(gfile, "position", "int32", storage.option, :
Invalid number of variants in 'position'.
Calls: seqMerge -> .append_node_variant

I have been careful to not split near variant positions. Any ideas? Thanks.

[davidroberson]

If the original file had multi-alleles on separate records would this cause the error shown above?

[davidroberson]

just a kind ping..please help :)

[zhengxwen]

Not sure what the issue is.
Did you split the bcf file by variants or samples?
If the variants are split, then seqMerge() could just simply concatenate the data.

Also, show me the session information: sessionInfo().