Question: bcftools merge is not merging ??
gravatar for abdulbenissa
2.2 years ago by
Saudi Arabia
abdulbenissa0 wrote:

Hello, Everyone !!

I have an easy one , I am sure I am missing something here. I am trying to merge a number of vcf files to one large vcf, utilising bcfttools merge.

the command I used is following

where is Samples is the list of vcf files as A.vcf.gz B.vcf.gz ,, etc)

$> bcftools merge --force-samples -l Sample > Meged.vcf.gz

the error I got is :

[W::bcf_hdr_check_sanity] GL should be declared as Number=G Warning: trying to combine "GQ" tag definitions of different types Error at chrM:309: wrong number of fields in CIGAR?


bcftools next-gen vcf • 2.0k views
ADD COMMENTlink modified 5 months ago by jaybee10 • written 2.2 years ago by abdulbenissa0
gravatar for Macspider
2.2 years ago by
Vienna - BOKU
Macspider2.8k wrote:

The bcftools merge help section says:

Usage:   bcftools merge [options] <A.vcf.gz> <B.vcf.gz> [...]

This means that you can specify file names either like above (in the last positional arguments of the command, one after each other, space separated) or in a file, which you can declare as -l.

-l, --file-list <file>             read file names from the file

You should write the --file-list as a file containing the file names you want to merge, including path, one per line. What you're doing is declaring them as a list in the command!

ADD COMMENTlink written 2.2 years ago by Macspider2.8k
gravatar for jaybee
5 months ago by
South Korea
jaybee10 wrote:

You can try /PATH/to/bcftools merge Home/data/*vcf.gz -Oz -o Merged.vcf.gz

[considering that your files A.vcf.gz , B.vcf.gz and so on are in the Home directory, under the folder data]

[If your bcftools is installed in the usr/bin directory then simply use:

/usr/bin/bcftools merge Home/data/*vcf.gz -Oz -o Merged.vcf.gz

This should work fine!

ADD COMMENTlink modified 5 months ago • written 5 months ago by jaybee10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 635 users visited in the last hour