Fixing genotypes from split vcf

I have split a multi-sample vcf into each sample, and I was wondering if there was a simple method for fixing the genotype for variants to 0/1 for each single-sample vcf, and changing the alternate allele to match the new genotype. For instance, if the vcf contains a read like

chr1 945122 . C A,T . PASS <INFO STRING> GT:ABQ:AD:ADF:ADR:DP:FREQ:GQ:PVAL:RBQ:RD:RDF:RDR:SDP 0/2:57:4:1:3:584:0.68%:0:9.8E-1:50:578:449:129:584

then I would like it to read

chr1 945122 . C T . PASS <INFO STRING> GT:ABQ:AD:ADF:ADR:DP:FREQ:GQ:PVAL:RBQ:RD:RDF:RDR:SDP 0/1:57:4:1:3:584:0.68%:0:9.8E-1:50:578:449:129:584

so that the alt allele and the genotype still match, just no other alternate alleles are in the VCF.

Is there a tool for tidying a vcf up like this?

vcf

• 1,448 views • link

written 6.1 years ago by graeme.thorn

Isn't this just a bcftools norm -m-any solution?

• 0 views • link

written 6.1 years ago by Kevin Blighe

Yeah, if done before splitting. Also, I'd recommend vt decompose over bcftools norm -m-any - vt retains information on the variants it "transforms".

• 0 views • link

written 6.1 years ago by Ram

You should split multi-allelic sites before splitting VCFs into sample-specific VCFs, if such "clean" genotypes are a necessity.

Use vt decompose to split multi-allelics, then any tool of your choice to get single sample VCFs.

By the way, an entry in a VCF file is a site/location/variant, not a "read".

• 0 views • link

written 6.1 years ago by Ram

0 answers

No answers yet.

More posts like this

No valid chromosomes found! on Michigan Imputation Server
written 2.8 years ago by graeme.thorn

I have two vcf files - one hg19 and one hg38, analysing data from the same participants on two slightly different SNP platforms. Both files …
Calculating Variant Allele Frequency
written 3.3 years ago by jason •

I got a VCF which I need to calculate variant allele frequency for each variant at each position. My understanding is that variant allele frequency …
bedtools intersect VCF with BED Invalid record error.
written 3.8 years ago by nicole.zeltser •

I am attempting to intersect a VCF file with coordinates from a BED file. e.g. bedtools intersect -a /path/to/vcf -b /path/to/.bed The command results in …
Merging VCFs from same sample created by different variant callers and adding new INFO field for va…
written 4.4 years ago by Jordi •

Hi, I am generating 2 distinct SNVs/small indels for each sample from WGS data. The two VCFs are generated by `gatk HaplotypeCaller` and `deepvariant`. After …
Extract vcf file header and save it in Key-Value csv
written 7.7 years ago by tinyskinn •

Hi everyone, If I have a vcf file with header like this: ``` ##fileformat=VCFv4.0 ##fileDate=20090805 ##source=myImputationProgramV3.1 ##reference=1000GenomesPilot-NCBI36 ##phasing=partial ##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data"> ##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth"> …
Get genotype GT from multi-sample VCF (per chromosome)
written 7.7 years ago by zx8754

Example input: multi-sample VCF (adapted from [www.internationalgenome.org][1]): **Note:** my actual file is bgzipped and tabixed with ~2Mln variants (rows) and ~1000 samples (columns). ##fileformat=VCFv4.0 ##fileDate=20090805 …
About structure variants simulation
written 8.7 years ago by Darren Lin •

Guys, I am trying to use DWGSIM to generate some simulation paired-end reads with a VCF file of candidate mutations. I use the command like …
add coordinates to line of output extracted from vcf
written 8.7 years ago by bioguy24

I am trying to extrapolate/convert variants from a vcf file to genotype. I am close but can not seem to add the chr and start …
Convert VCF dosage file into dosage only file
written 9.0 years ago by dp0b •

Hi, Is it possibly to only output dosages from a generated vcf dosage file from minimac? Input format is ##fileformat=VCFv4.1 ##filedate=2017.7.5 ##source=Minimac3 ##contig=<ID=29> ##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> ##FORMAT=<ID=DS,Number=1,Type=Float,Description="Estimated …
Calculating Variant Allele Frequency
written 9.7 years ago by novice

I'm reading an interesting paper, [Malachi et al., Cell Systems 2015][1], that talks about Variant Allele Frequency (VAF). Could someone please help me understand how …