Question: VCF output when merging with bcftools merge
1
gravatar for steven.davis
3.0 years ago by
steven.davis10
steven.davis10 wrote:

I am merging multiple single-sample VCF files into one multi-sample VCF file with bcftools merge using this command:

bcftools merge --info-rules NS:sum -o outFilePath dir/*.gz

Previously, using bcftools 1.1, the output VCF file contained one row per position, like this:

#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO
barref  1094    .       A       G       .       PASS    NS=24
barref  1221    .       G       A       .       PASS    NS=24
barref  2015    .       G       A       .       PASS    NS=24

Using newer version of bcftools 1.2 or 1.3.1, I am seeing two rows per position, like this:

#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO
barref  1094    .       A       .       .       PASS    NS=20
barref  1094    .       A       G       .       PASS    NS=4
barref  1221    .       G       .       .       PASS    NS=20
barref  1221    .       G       A       .       PASS    NS=4
barref  2015    .       G       .       .       PASS    NS=20
barref  2015    .       G       A       .       PASS    NS=4

First, should I be concerned about the change in output?

Secondly, I would prefer the old style with a single row per position. Does anyone know how to get one row of output per position when using bcftools 1.2 or higher?

bcftools vcf • 2.0k views
ADD COMMENTlink modified 3.0 years ago • written 3.0 years ago by steven.davis10
0
gravatar for steven.davis
3.0 years ago by
steven.davis10
steven.davis10 wrote:

Solved it with help from a colleague. Problem solved with "--merge all" command line option.

ADD COMMENTlink written 3.0 years ago by steven.davis10
Please log in to add an answer.

Help
Access

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