Hey can you post the full command you used to solve this problem?
• 0 views
•
link
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?
Log in to answer this question.