Compare the amino acid changes of supplied mutations to cosmic and omim databases.
This document describes genome music cosmic-omim (2013-05-14 at 16:03:04)
genome music cosmic-omim --maf-file=? --output-file=? --reference-build=? [--omimaa-dir=?] [--cosmic-dir=?] [--verbose] [--wu-annotation-headers] [--aa-range=?] [--nuc-range=?] [--show-known-hits]
... music cosmic-omim \
--maf-file input_dir/myMAF.tsv \
--output-file output_dir/myMAF_output.tsv \
--no-verbose
... music cosmic-omim \
--maf-file input_dir/myMAF.tsv \
--output-file output_dir/myMAF_output.tsv \
--omimaa-dir omim_dir/ \
--cosmic-dir cosmic_dir/ \
--no-verbose
list of annotated mutations in \s-1MAF\s0 format (or any file with MAF+annotation headers)
Output file contains the input file with two columns appended to the end, corresponding to cosmic and omim mutation comparisons, respectively
Put either \*(L"Build36\*(R" or \*(L"Build37\*(R" Default value 'Build37' if not specified
omim amino acid mutation database folder
cosmic amino acid mutation database folder
Use this to display the larger working output Default value 'false' (--noverbose) if not specified
Use this if input \s-1MAF\s0 contains \s-1WUSTL\s0 annotation format headers Default value 'false' (--nowu-annotation-headers) if not specified
Set how close a 'near' match is when searching for amino acid near hits Default value '2' if not specified
Set how close a 'near' match is when searching for nucleotide position near hits Default value '5' if not specified
When a finding is novel, show known \s-1AA\s0 in that gene Default value 'true' if not specified
This tool looks at the amino acid changes for the given set of mutations and compares the genomic coordinates as well as the affected amino acid to the coordinates and amino acids of all cancer-specific mutations listed in the Cosmic and \s-1OMIM\s0 databases. The database files are specially prepared for this task and provided with the MuSiC suite. The tool reports various types of matches, including matches within \*(L"near proximity\*(R", where \*(L"near proximity\*(R" is currently defined as a linear \s-1DNA\s0 distance of 5 bases or 2 amino acids. (This type of matching helps to account for the possibility of subtle differences in reported positions for variants due to differences in transcript definitions or other things of this nature.) Any site without a match in a particular databases is reported as \*(L"novel\*(R" with respect to that database.
The output of this script returns each row the original input \s-1MAF\s0 file with two columns appended to the end of each, one column for each of the databases. Also included is a \s-1STDOUT\s0 printout of a summary of what was found in the input \s-1MAF\s0. Neither output can be suppressed in the current version. The --verbose option is used to display working notes that are useful for various purposes in debugging potential \s-1MAF\s0 problems. The Omim and Cosmic directories must point to the output of the downloader, named appropriately, as they don't recognize the \s-1OMIM\s0 database in the raw download format.
This tool only compares build 36 or build 37 coordinates that are specified in Cosmic to the coordinates in your maf file. This is a weakness of the Cosmic database (not all position entries are currently available for both builds) that is out of our control.
In addition to the standard version 2.3 \s-1MAF\s0 headers, there needs to be 3 columns appended. These column headers in the \s-1MAF\s0 must have these names in the header in order for the tool to find them:
transcript_name - the transcript name, such as \s-1NM_000028\s0
amino_acid_change - the amino acid change, such as p.R290H
c_position - the nucleotide position changed, such as c.869
Copyright (C) 2010-2011 Washington University in St. Louis.
It is released under the Lesser \s-1GNU\s0 Public License (\s-1LGPL\s0) version 3. See the associated \s-1LICENSE\s0 file in this distribution.
Brian Dunford-Shore David Larson, Ph.D. Michael C. Wendl, Ph.D. William Schierding, M.S.
This tool depends on copies of data from the following databases, packaged in a form useable for quick analysis:
* COSMIC - http://www.sanger.ac.uk/genetics/CGP/cosmic/ * OMIM - http://www.ncbi.nlm.nih.gov/omim
genome-music(1), genome(1)