Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

"bcf_synced_reader.h:108 operator()" error wen using vt sort #110

Open
estolle opened this issue Apr 2, 2020 · 0 comments
Open

"bcf_synced_reader.h:108 operator()" error wen using vt sort #110

estolle opened this issue Apr 2, 2020 · 0 comments

Comments

@estolle
Copy link

estolle commented Apr 2, 2020

Hi there

I was running an older script which previously worked, but now threw an error at the step when I use vt sort on a VCF file:

vt sort $TF/test2.vcf
[E:bcf_synced_reader.h:108 operator()] alleles not defined in VCF record.
There is nothing I found being wrong with the VCF, I changed headers, sorted, checked alleles and genotypes, tested on subsets and other VCFs. It always comes up for a reason unclear to me.

Seems its similar to this issue:
DecodeGenetics/graphtyper-pipelines#3 (comment)
This describes that an older version works (and that commit fixed some htslib stuff). Indeed I never changed the previously working vt installation (vt v0.57721) but I updated htslib. So perhaps the problem is related to that.

Other tools work:
bcftools sort -Ov $TF/test2.vcf
no error

vcfstreamsort -a $TF/test2.vcf
no error

vt validate $TF/test2.vcf or vt peek $TF/test2.vcf work as expected too.

In the meantime I am using bcftools or vcflib to sort the vcf, but perhaps the issue in vt can be fixed.

ning-y added a commit to ning-y/docker-vt that referenced this issue Jun 23, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant