merging tetraploid vcf files through bcftools
1
1
Entering edit mode
19 days ago
analyst ▴ 20

I have 7 tetraploid vcf files (generated through gatk Hapotypecaller) that I tried to merge through bcftools but It gave me following error:

Incorrect number of FORMAT/PL values at LK031787:1508, cannot merge. The tag is defined as Number=G, but found 5 values and 2 alleles.

This is what I have in position 1508:

enter image description here

Whereas my comment section includes following information for TAG PL:

##FORMAT=<ID=PL,Number=G,Type=Integer,Description="Normalized, Phred-scaled likelihoods for genotypes as defined in the VCF specification">

vcf bcftools • 338 views
ADD COMMENT
0
Entering edit mode

I replaced G with . but now got this error:

[W::bcf_hdr_check_sanity] PL should be declared as Number=G

ADD REPLY
0
Entering edit mode

Is it necessary to perform imputation for GBS data ?

ADD REPLY
2
Entering edit mode
19 days ago
Axzd ▴ 70

Can you show your merging command?

Also check you are using last version of bcftools

Finally normalise your VCF before merging.

ADD COMMENT
0
Entering edit mode

Thanks Axzd !

bcftools version 1.15.1

I first bgzipped vcf files using command:

bgzip -c file1.vcf > file1.vcf.gz

Then I used following command for indexing:

tabix -C -h -p vcf file1.vcf.gz

Finally performed merging as follows:

bcftools merge -l input_file.txt -O v -o merged_snps.vcf

ADD REPLY
1
Entering edit mode

I have fixed issue by removing PL TAG. I hope it is not required for filtration and analysis.

ADD REPLY

Login before adding your answer.

Traffic: 1128 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