I used GATK HaplotypeCaller for Variant calling of 2300 MTB strains. Now i want to make it a multi Sample VCF.
I used CombineVariants but im comfuse that either this is correct or not.
Each of my VCF is isolated means 1 VCF from one individual sample. But GATK says: CombineVariants can be used for combine variant calls that were produced from the same samples but using different methods, for comparison.
But my Variant files are from different samples. I just want to make a multi VCF for comparison purpose of each strains variants. The VCF file generated from CombineVariant is :
Mutltisample.vcf I'm not sure is it giving the proper results?
And how to interpret this union.vcf. Like in GT column it is giving ./. what does that means? and its a haploid genome so the GT should contain a single value but few columns are representing the 1|0, 0|1...
Can anyone help?
Thank you
•
written
by
S AR
2 answers
You could merge the files with bcftools:
bcftools merge -O z -o merged.vcf.gz sample1.vcf.gz sample2.vcf.gz
But GATK says: CombineVariants can be used for combine variant calls that were produced from the same samples but using different methods, for comparison.
yes CombineVariants can merge vcf with the same samples (--genotypeMergeOptions PRIORITIZE ) but in your case, you need to use the option --genotypeMergeOptions REQUIRE_UNIQUE
Require that all samples/genotypes be unique between all inputs.
Hey Pierre, do you have any advice for using GenotypeGVCF for multiple samples from viruses that are similar but have different reference strains. ie in order to improve read alignment slightly different reference strains were used. But I would like to compare them after even though the reference strains are slightly different. thank you
Log in to answer this question.
-
how to extract variants from the vcf.gz files linked below?
written by Mahad •
-
Error occured at the PHG imputation step: allele in genotype not in the variant context
written by xeod21 •
-
What happens to the genotype field in a g.vcf after CombineGVCFs
written by nhaus •
-
VCF files never get annotated with PG tag for CalculateGenotypePosteriors
written by lait
-
Annotating Mycobacterium tuberculosis VCF file using snpEFF and aNNOVAR
written by S AR
-
How to merge/combine multiple vcf from different samples to make one multi sample vcf file using GA…
written by S AR
-
SRA files bulk downloads
written by S AR
-
Tools to filter Mutect2 output
written by nanana
-
Multi-sample variant calling
written by sumudu_rangika •
-
VCF format conventions for multi-sample vs. single-sample data
written by bberns •
