vcf file has an empty sample id
0
0
Entering edit mode
4.2 years ago
yh362 ▴ 50

It seems like my VCF's header contains at least an empty smaple ID. I am not sure how to figure out where it goes wrong because, if I run awk 'print NF', all lines have the same number of fields.

When I ran bcftools, it told me "[E::bcf_hdr_add_sample] Empty sample name: trailing spaces/tabs in the header line?"

And when I ran refinedIBD, it gives an error of the same reason

> > java -jar ../../../../tools/refined-ibd.17Jan20.102.jar gt=test.vcf out=test ERROR: java.lang.IllegalArgumentException: id.isEmpty()
> java.lang.IllegalArgumentException: id.isEmpty()
>         at beagleutil.SampleIds.getIndices(SampleIds.java:84)
>         at beagleutil.Samples.fromIds(Samples.java:83)
>         at vcf.VcfHeader.samples(VcfHeader.java:172)
>         at vcf.VcfHeader.<init>(VcfHeader.java:127)
>         at vcf.RefIt.<init>(RefIt.java:120)
>         at vcf.RefIt.create(RefIt.java:107)
>         at refined.RefinedIbd.lambda$windowIt$0(RefinedIbd.java:140)
>         at vcf.WindowIt$Reader.run(WindowIt.java:244)
>         at java.lang.Thread.run(Thread.java:745) terminating program.

Yet I couldn't locate which column in the #header line goes wrong. Can anyone help? Thanks a lot!

vcf bcftools SNP • 1.3k views
ADD COMMENT
0
Entering edit mode

How you prepared this vcf file? BAM file or from internet source?

ADD REPLY
0
Entering edit mode

Could you attach header in your question?

ADD REPLY

Login before adding your answer.

Traffic: 2673 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6