Is it possible to combine single sample VCF files to multi sample VCF files? The samples are different strains from the same species, mapped against the same reference.
The single sample vcf files have been made with GATK and the multi sample vcf file should preferentially also be compatible with GATK for downstream processing.
We chose not to do multi-sample calling in GATK because we don't want the genotype calls and their quality from one sample to be influenced by the genotype calls and their quality from other samples (which GATK does by default when doing multi-sample calling ).
We do however want to to process the samples the and their genotypes as one big matrix.
•
written
by
William
3 answers
The problem is that unless your vcf's have a line for every single locus, samples that match consensus at a position will not have an entry for that position, so if you have a SNP in one sample at a given position, other samples will have no vcf entry for that position, and you won't know if your other samples really match consensus there, or if they have low coverage there, and a variant couldn't be safely called.
That's why calling them all together is preferable. If one variant is iffy in one sample, the fact that that variant exists in other samples makes it more likely that the the iffy variant is real.
use vcf-merge
•
written
by
Pierre Lindenbaum
Documentation migrated to https://vcftools.github.io/perl_module.html#vcf-merge
•
written
by
setempler •
You can also try Gene-Talk.de, upload the single files, and create a collection from all. afterwards you can filter and reduce your data to clinically relevant variants.
Up to now GeneTalk supports only human hg19 referenced VCF files.
•
written
by
alexej.knaus •
Log in to answer this question.
-
quality check of vcf files
written by rthapa •
-
How to combine a multi sample VCF from multiple sample VCF and interpret a multi sample vcf of bact…
written by S AR
-
compare and merge VCF files
written by qwzhang0601
-
How does samtools mpileup handle multiple .bam inputs?
written by traviata •
-
Multi-sample variant calling
written by sumudu_rangika •
-
Do freeBayes and platypus support gVCF (ie incremental variant calling for large sets of genomes)
written by William
-
Multi-Sample Vcf To Phylogenetic Tree.
written by William
-
Gatk, Variant Quality Score Recalibration, How To Work With Custom Truth And Training Datasets?
written by William
-
Gatk Multi-Sample Vcf Variantfiltration
written by William
-
Extract Private Snps From Multi-Sample Vcf File
written by William
