Hello support,

I am merging VCF files 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

How can I get one row of output per position when using bcftools 1.2 or higher?

Thanks,

Steve Davis
Office of Analytics and Outreach
FDA Center for Food Safety and Applied Nutrition
5100 Paint Branch Pkwy
College Park, MD 20740
Office: 240-402-4834

------------------------------------------------------------------------------
_______________________________________________
Samtools-help mailing list
Samtools-help@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/samtools-help

Reply via email to