Question: Converting Vcf File To Maf
7
gravatar for l.pirhaji
4.4 years ago by
l.pirhaji70
l.pirhaji70 wrote:

I'm going to use VarScan2 and MuSiC software to identify mutated genes from exome sequencing data.

at first, I used VarScan2 to identity somatic mutations. I will then use MuSiC to identify significantly mutated genes from VarScan2 results. However, the output file generated from VarScan2 is in VCF format, but MuSiC only accepts MAF file format as an input.

I was wondering if there is any available tools for converting VCF file format to MAF.

music • 9.8k views
ADD COMMENTlink modified 9 days ago by Biostar ♦♦ 20 • written 4.4 years ago by l.pirhaji70
3

Possible duplicate of VCF to MAF (Mutation Annotation Format) Conversion ?

ADD REPLYlink written 4.4 years ago by Eric Normandeau9.6k
1
gravatar for mdm-two
4.4 years ago by
mdm-two200
mdm-two200 wrote:

http://code.google.com/p/ngs-analysis/source/browse/modules/somatic/vcf2maf.py

Will read in an annotated vcf file outputted by SNPEff It will parse and output the data in tcga maf format (currently v 2.3 http://goo.gl/2Shtq)

By default, every effect and transcript annotated per variant will be outputted as a row in the maf file. (In other words you will need to filter to one row per sample per variant)

This included somatic caller support for varscan and GATK.

There are many other vcf2maf alternatives such as

https://github.com/Sage-Bionetworks/vcf2maf/blob/master/vcf2maf/vcf2maf.py

You might also want to research Ensembl VEP for VCF annotation:

http://useast.ensembl.org/info/docs/variation/vep/index.html

ADD COMMENTlink written 4.4 years ago by mdm-two200
1
gravatar for Cyriac Kandoth
4.0 years ago by
Cyriac Kandoth4.7k
Memorial Sloan Kettering, New York, USA
Cyriac Kandoth4.7k wrote:

I recently posted a VCF->MAF conversion script at: github.com/ckandoth/vcf2maf. It's plenty documented so that you understand what information is lost in translation.

Briefly - each VCF variant must be annotated to only one of all possible gene transcripts/isoforms that it might affect. This selection of a single affected transcript/isoform per variant, is often subjective. For now, the scripts tries to follow best-practices: it chooses the "worst" effect on the "best" transcript. If there are multiple such candidates, it annotates the variant effect on the longest "best" transcript.

ADD COMMENTlink written 4.0 years ago by Cyriac Kandoth4.7k

Thanks for the script but the current version still requires snpEff yet I have already annotated using snpEff. Could you please push a version that doesn't depend on snpEff? Thanks again!

ADD REPLYlink written 3.7 years ago by tayebwajb80
2

The script skips snpEff annotation for an input VCF named file.vcf if it finds an annotated VCF named file.anno.vcf. You can use this as a workaround. Otherwise, try this fork of the script: https://github.com/dakl/vcf2maf

Update: I released vcf2maf v1.1 that allows you to use a VCF that is already annotated with snpEff or Ensembl's VEP.

ADD REPLYlink modified 3.7 years ago • written 3.7 years ago by Cyriac Kandoth4.7k
0
gravatar for mdm-two
4.2 years ago by
mdm-two200
mdm-two200 wrote:

Same question answered here: http://seqanswers.com/forums/showthread.php?t=16740 converting VCF file to MAF

ADD COMMENTlink written 4.2 years ago by mdm-two200
0
gravatar for rse
19 months ago by
rse30
United States
rse30 wrote:

Hi Cyriac, i used vcf2maf perl script. Installed vcf2maf-master, vep. But getting an error - WARNING: Chromosome 1 not found in cache on line . What can be the reason?

ADD COMMENTlink written 19 months ago by rse30

Please don't write a question in a section that's meant for an answer. Either start a new question, or submit an issue at https://github.com/mskcc/vcf2maf/issues

ADD REPLYlink written 19 months ago by Cyriac Kandoth4.7k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 778 users visited in the last hour