Combine Single Sample Vcf To Multi Sample Vcf
3
5
Entering edit mode
11.2 years ago
William ★ 5.3k

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.

gatk snp vcf • 18k views
ADD COMMENT
6
Entering edit mode
11.1 years ago

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.

ADD COMMENT
1
Entering edit mode

I did use EMIT_ALL_CONFIDENT_SITES for the output mode of the unified genotyper. So each sample vcf also has reference calls.

ADD REPLY
5
Entering edit mode
11.2 years ago

use vcf-merge

ADD COMMENT
1
Entering edit mode
ADD REPLY
1
Entering edit mode
10.9 years ago
alexej.knaus ▴ 130

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.

ADD COMMENT

Login before adding your answer.

Traffic: 2905 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6