This is a test version of Biostars. For the public version, visit https://www.biostars.org.
vcf file format

I used plink2 to make a vcf file from my .bed, .bim, and .fam files. This produced a vcf v4.3 file. Just realized that vcftools only supports 4.0, 4.1, and 4.2. How can I convert my file to a lower version?

vcf
• 5,213 views • link
updated by richyanicky • richyanicky • written by johnja • johnja

Oh dear, VCF versions rarely have tangible change logs or backwards compatible tools. I had to go from VCF4.2 to a rough 4.1 by removing all problematic loci I found, no clue what's changed in 4.3 :-(

• 0 views • link
written by Ram Ram

Thank you. I actually just downloaded plink version 1.9 instead of 2.0 and then made the vcf file again from the original data.

• 0 views • link
written by johnja • johnja

Have you actually tried whatever you want to do? It may still very well work. The changes in v4.3 aren't all that extreme and mostly consist of a few new tags and making hard lines about how certain tags should be formatted and named.

I'd try it first. At worst, you might have to change some of the header lines.

• 0 views • link
written by jared.andrews07 jared.andrews07

2 answers

In plink2 there is a vcf-4.2 option .

  • 'vcf', 'vcf-4.2': VCF (default version 4.3). If PAR1 and PAR2 are present, they are automatically merged with chrX, with proper handling of chromosome codes and male ploidy. When the 'bgz' modifier is present, the VCF file is block-gzipped. The 'id-paste' modifier controls which .psam columns are used to construct sample IDs (choices are maybefid, fid, iid, maybesid, and sid; default is maybefid,iid,maybesid), while the 'id-delim' modifier sets the character between the ID pieces (default '_'). Dosages are not exported unless the 'vcf-dosage=' modifier is present. The following five dosage export modes are supported: 'GP': genotype posterior probabilities (v4.3 only). 'DS': Minimac3-style dosages, omitted for hardcalls. 'DS-force': Minimac3-style dosages, never omit. 'HDS': Minimac3-style phased dosages, omitted for hardcalls and unphased calls. Also includes 'DS' output. 'HDS-force': Always report DS and HDS.
• 0 views • link
written by richyanicky • richyanicky

EDIT: As genomax said, the VCFv4.3 doc seems to include the changes from 4.2 - I guess they did not make the same omission twice :-)

The tip mentioned here might be of use if you wish to find the difference between v4.3 and v4.2: https://bioinformatics.stackexchange.com/questions/344/whats-the-difference-between-vcf-spec-versions-4-1-and-4-2

you can download both specs in .tex format and do diff.

• 0 views • link
written by Ram Ram

Log in to answer this question.

More posts like this