This is a test version of Biostars. For the public version, visit https://www.biostars.org.
merging multiple vcfs by column in Picard

Hello everyone,

I used picard MergeVcfs to combine individual vcf files with the following command:

java -jar /picard.jar MergeVcfs -I input.vcf -I input2.vcf -I inputs.vcf  O= output.vcf.gz

The resulting output file looks like this: 

#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  20
Chromosome_1    517     .       CTT     C       79.60   .       AC=1;AF=0.500;AN=2;BaseQRankSum=-0.080;DP=16;ExcessHet=3.0103;FS=0.000;MLEAC=1;MLEAF=0.500;MQ=57.36;MQRankSum=-3.750;QD=4.98;ReadPosRankSum=0.068;SOR=0.446     GT:AD:DP:GQ:PL  0/1:13,3:16:87:87,0,537
Chromosome_1    562     .       CATTTCTCTA      C       46.60   .       AC=1;AF=0.500;AN=2;BaseQRankSum=1.180;DP=14;ExcessHet=3.0103;FS=0.000;MLEAC=1;MLEAF=0.500;MQ=53.19;MQRankSum=0.558;QD=3.88;ReadPosRankSum=0.493;SOR=1.179       GT:AD:DP:GQ:PL  0/1:10,2:12:54:54,0,404

The resulting file is 20.1gb and samples seemed to be combined by rows and not columns. How do I change it so that I can reduce the file size of my vcf?

Thank you!

vcf snp genome picard
• 3,011 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by nitinra • nitinra

2 answers

The resulting file is 20.1gb and samples seemed to be combined by rows and not columns.

I suppose all your VCF have the same sample named 20. They should all have some different names.

use bcftools reheader to rename the samples or GATK3.8 CombineVariants with --genotypeMergeOptions UNIQUIFY

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

"Merge" can mean (at least) two operations with VCF data. Using bcftools notations, you can either concat VCF data that describe variants occurring in different regions of the same sample, or you can merge data that describe variants occurring in the same region in diiferent samples. I think your input data determines what Picard's MergeVcfs does, and in your case, the sample names could be the same, resulting in the tool working in a concat fashion.

• 0 views • link
written by Ram Ram

Log in to answer this question.

More posts like this
  • Annotate genes from gtf file to vcf file
    written by QX • QX
    Hi all, Is there any tool that can help to annotate the genes from GTF file to a (columns) in the vcf file? ``` #CHROM …
  • why CombineGCVFs in gatk not showing all the samples name?
    written by rj.rezwan • rj.rezwan
    Hi, I combined the 64 `.vcf` files using the `CombineGVCFs` in `gatk`. The command was completed successfully but its showing the output with only one …
  • HaplotypeCaller GENOTYPE GIVEN ALLELES doesn't genotype given alleles
    written by kamanovae • kamanovae
    Hi! I am trying to run the Gatk HaplotypeCaller (human data): ``` ./gatk-4.1.2.0/gatk HaplotypeCaller\ --reference ref.fa \ --input file.bam \ --genotyping-mode GENOTYPE_GIVEN_ALLELES \ --alleles allele_chunk_file.vcf …
  • PGT only available for some variants in GATK .vcf
    written by michael.flower.14 • michael.flower.14
    I've got a vcf file someone else prepared using GATK. I'm interested in the phasing information in the PGT tag e.g. 0|1. This information seems …
  • Extract the DP from VCF file along with the chromosome postion and alteration
    written by manojkumarbioinfo manojkumarbioinfo
    Hi, I would like to extract the DP from my VCF file, along with the chromosome position and alteration Example VCF file CHROM POS ID …
  • How is allelle frequency calculated ?
    written by DareDevil DareDevil
    How is allelle frequency calculated ? This a variation where Allelle Frequency (AF) is 0.500 and AD = 57,36 and DP = 93 My understanding …
  • Strange Allele Frequency in vcf file
    written by GP_PG • GP_PG
    Hi everyone. I performed a snp calling using GATK and following all the suggestions in the BestPractices sections. I performed single snp calling for two …
  • why is there "." in Filter?
    written by arshiya.g02 • arshiya.g02
    This was the command I used `sudo ./gatk HaplotypeCaller -R /home/arshiya/efs/FastQC/hg38.fa -I /home/arshiya/efs/FastQC/sample_sort_reads.bam -O /home/arshiya/efs/FastQC/sample_haplo.vcf.gz` and the following was the output I got chr1 7820343 …
  • Diploid calls on chrY?
    written by dec986 dec986
    Hello, I recently got my genome sequenced. chrY 59030922 . A G 1244.77 PASS AC=1;AF=0.500;AN=2;BaseQRankSum=-6.580e-01;ClippingRankSum=0.00;DP=41;ExcessHet=3.0103;FS=9.230;MLEAC=1;MLEAF=0.500;MQ=37.33;MQRankSum=-4.563e+00;QD=30.36;ReadPosRankSum=-7.570e-01;SOR=0.076;VQSLOD=0.452;culprit=FS GT:AD:DP:GQ:PL 0/1:9,32:41:99:1273,0,180 chrY 59032665 . A G 418.77 PASS AC=1;AF=0.500;AN=2;BaseQRankSum=-4.182e+00;ClippingRankSum=0.00;DP=52;ExcessHet=3.0103;FS=0.000;MLEAC=1;MLEAF=0.500;MQ=45.59;MQRankSum=-6.082e+00;QD=8.21;ReadPosRankSum=-1.528e+00;SOR=0.756;VQSLOD=2.78;culprit=DP …
  • Large VCF QUAL Values
    written by pld pld
    I'm currently working on calling variants for a virus that was sequenced. For lack of a better solution, I'm using the GATK UnifiedGenotyper pipeline. I'm …