How to merge multiple patient's vcf files (indel and snv) with different IDs?

Hi all, I have some VCF files for my patients, each patient has 2 files( indel.vcf , snv.vcf) and I want to merge these file by the script bellow:

java -jar gatk-package-4.2.0.0-local.jar  MergeVcfs -I /PATH_TO_patient1_ID_indel.vcf  -I /PATH_TO_patient1_ID_snv.vcf -O /PATH_TO_patient1_ID_merged.vcf

I don't know how to do this for multiple patients with different IDs.

Thanks for any help.

vcf gatk

• 1,262 views • link

updated 5.0 years ago by German.M.Demidov

• written 5.0 years ago by Zahra •

bcftools concat may help

• 0 views • link

written 5.0 years ago by German.M.Demidov

0 answers

No answers yet.

More posts like this

ASEReadCounter output wrong number of coverage
written 2.7 years ago by junhuili •

Hi, I am using ASEReadCounter to count the number of reads per variant in a BAM file. For some positions, it will report 1 read …
How to analyze whole-genome data using Picard and GATK?
written 4.6 years ago by Dan •

Hello, I am new to whole-genome sequencing data analysis. I am not sure whether I am correctly using Picard and GATK for my whole-genome analysis. …
How to give multiple input files to the gatk ?
written 5.0 years ago by Zahra •

Hi all, I want to use gatk-package-4.2.0.0 for updating my VCF dictionary. I use the script below for just one VCF file, and it has …
Reference file for base recalibration
written 5.0 years ago by priya.bmg •

Hello I need help in the Base recalibration step in variant analysis. I downloaded the fasta reference file for the recalibration step and there seems …
Merge multiple VCF files together with different numbers of samples
written 5.2 years ago by hemr3 •

I am aware that there are other questions on this topic. However, the answers aren't working for me. I have 22 VCF files that I …
About shell script
written 5.5 years ago by r00628112 •

Good afternoon everyone: I saved the script for each sam file in a A.txt file as below: java -jar /data/share/tools/gatk-4.1.7.0/gatk-package-4.1.7.0-local.jar SortSam -I DRR212437.sam -O DRR212437.sort.bam …
GATK tool to merge INDEL with SNPs with the same set of samples
written 5.5 years ago by MAPK

Hi All, I am trying to merge two VCF files, one with SNPs and the other with INDELS. I was looking at these three methods, …
GATK variant calling on multiple bam files
written 5.8 years ago by User000

Hello, I am trying to do variant calling using GATK on a list of several hundreds of bam files like this: bob.dup.bam smith.dup.bam will.dup.bam This …
Error while merging multiple vcf files using Picard MergeVcfs
written 5.9 years ago by nitinra •

Hello all, I am trying to merge multiple vcf files using MergeVcfs from Picard tools. Here is my command: java -jar ./picard.jar I=input1.vcf I=input2.vcf I=inputx.vcf …
How to get vcf not from g.vcf?
written 7.1 years ago by woojoy14 •

java -jar gatk-package-4.x.x.x-local.jar HaplotypeCaller -R ucsc.hg19.fasta -I SRR000982.mapped.sorted.markdup.recal.bam -O SRR000982.g.vcf -ERC GVCF I've been following the pipeline based on GATK BestPractice and I was told …