That's a great tool, thanks! I added a command line parameter for the name of snpeff vcf, feel free to use it if interested. (https://github.com/dakl/vcf2maf)
Vcf To Maf (Mutation Annotation Format) Conversion ?
2 answers
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 script 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 canonical "best" transcript.
• 115 views
•
link
written
by
Cyriac Kandoth
Yea that makes sense - to give the user the option to run snpEff themselves. Actually, the first version of this script was a "converter of a pre-annotated VCF" :) Then I wanted to package it all-in-one.
Update: I released vcf2maf v1.1 that allows you to use a VCF that is already annotated with snpEff or Ensembl's VEP.
• 225 views
•
link
written
by
Cyriac Kandoth
FYI, I recently started getting ERROR: Unrecognized biotype "non_coding". Please update your hashes! at vcf2maf.pl line 287, <GEN0> line 171..I added it with priority 3 which had other non-coding RNAs in it. Just so you know.
Info on the biotype from here: http://vega.sanger.ac.uk/info/about/gene_and_transcript_types.html
Thanks. Which transcript database are you using? I don't see non_coding as a valid transcript biotype in the Ensembl 74 GTF, but I do see it listed in the GENCODE specs. I have now updated the script to handle all the GENCODE biotypes.
• 109 views
•
link
written
by
Cyriac Kandoth
If a variant locus maps to multiple genes/transcripts, which biotype is most well defined and/or more likely to be disease associated.
• 0 views
•
link
written
by
Cyriac Kandoth
Please see fork of the code mentioned above by @Danielk. Alternatively, my script skips snpEff annotation for an input VCF named file.vcf if it finds an annotated VCF in the same folder named file.anno.vcf.
Update: I released vcf2maf v1.1 that allows you to use a VCF that is already annotated with snpEff or Ensembl's VEP.
• 64 views
•
link
written
by
Cyriac Kandoth
MAF contains annotation about the variant effects on transcripts/proteins while VCF typically does not. You might find that using tools like annovar, snpeff, and the Ensembl Variant Effect Predictor get you pretty close. I'm not aware of a script that applies one or more of the tools to a VCF file to produce MAF directly.
• 1 views
•
link
written
by
Sean Davis
But you'll probably still need to run annovar or snpeff or something like that (unless you are into reinventing wheels). The output of annovar or snpeff is what gets fed to your script is what I would envision.
• 0 views
•
link
written
by
Sean Davis
I should comment here that MAF is not really considered a "standard" format, so you may want to make sure that the output of one of the software packages mentioned above would not suffice for your final purpose.
• 251 views
•
link
written
by
Sean Davis
FWIW, MAF is a "standard" format within the TCGA project. Here's documentation: https://wiki.nci.nih.gov/display/TCGA/Mutation+Annotation+Format+(MAF)+Specification
• 129 views
•
link
written
by
Chris Miller
Thanks Sean. The problem started when I wanted to use MuSiC gmt.genome.wustl.edu/genome-music/0.2/index.html) This requires the mutations in MAF format and I have a bunch of vcfs. You are right, that I first need to extract information from the vcf through annovar/snpeff.
If there are two variant alleles, then you will find that in the ALT column of the VCF file as a comma-separated value. In most cases, there will not be a second variant allele present, I do not think.
• 0 views
•
link
written
by
Sean Davis
Log in to answer this question.
More posts like this
-
Is there any method to convert a text file to VCF format?
written by Abbas.M •
-
VCF file (annotated with snpeff) to MAF
written by JJ
-
Whole exomes - More than 3 time points
written by Kasthuri
-
Variant Control Format (VCF) TO Mutation Annotation Format (.MAF)
written by genomics Newbie •
-
How to convert from a .txt to a .maf file?
written by ksanchez1296 •
-
Converting Vcf File To Maf
written by anshupa.vssut •
-
Tool or Script to convert MAF file to VCF ?
written by ttom
-
How to get the mutated protein sequence from the DNA mutation information like the VCF or MAF file?
written by bioinflix •
-
Visualization Of Mutect Output For Point Mutation Detection
written by jeremy
-
Whole Exome From Whole Genome
written by Kasthuri
See www.biostars.org/p/74822/ and seqanswers.com/forums/showthread.php?t=16740
I'm afraid your pointers are not useful here:
I have snpeff annotated vcf files and I am converting these to maf format. When I run vcf2maf I get the rerror
Can you please point out the reason for this error.
Please open a new question, and use tags and keywords like vcf, maf, vcf2maf... so the relevant folks can find it.