Is there a way to reformat vcf 4.2 down to 4.1?
1
0
Entering edit mode
8.6 years ago
umn_bist ▴ 390

I've been using MuTect2 (part of GATK toolkit) to call somatic variants for my RNAseq data set. The vcf file is in 4.2 format which I'm not familiar with. Is there any advantage to using 4.2 instead if I'm annotating with snpEff/annovar and filtering with snpSift?

The documentations for annotation/filtering tools have been based on 4.1 vcf format which has proven difficult to follow for my particular vcf files.

Any help will be greatly appreciated.

variant calling snpSift snpEff • 6.6k views
ADD COMMENT
1
Entering edit mode
8.6 years ago

Try this script (courtesy of Daniel Cook).

ADD COMMENT
0
Entering edit mode

How does one execute this script?

ADD REPLY
0
Entering edit mode

Put it into your .bashrc, source that file again (source ~/.bashrc) and execute as vcf_downgrade oldfile.vcf

ADD REPLY

Login before adding your answer.

Traffic: 2062 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