OnWorks favicon

gmt-music-bmr-calc-bmrp - Online in the Cloud

Run gmt-music-bmr-calc-bmrp in OnWorks free hosting provider over Ubuntu Online, Fedora Online, Windows online emulator or MAC OS online emulator

This is the command gmt-music-bmr-calc-bmrp that can be run in the OnWorks free hosting provider using one of our multiple free online workstations such as Ubuntu Online, Fedora Online, Windows online emulator or MAC OS online emulator



gmt music bmr calc-bmr - Calculates mutation rates given per-gene coverage (from "music
bmr calc-covg"), and a mutation list


This document describes gmt music bmr calc-bmr version 0.04 (2016-01-01 at 23:10:19)


gmt music bmr calc-bmr --bmr-output=? --roi-file=? --gene-mr-file=? --reference-sequence=?
--bam-list=? --output-dir=? --maf-file=? [--skip-non-coding] [--skip-silent]
[--bmr-groups=?] [--show-skipped] [--separate-truncations] [--merge-concurrent-muts]

... music bmr calc-bmr \
--bam-list input_dir/bam_list \
--maf-file input_dir/myMAF.tsv \
--output-dir output_dir/ \
--reference-sequence input_dir/all_sequences.fa \
--roi-file input_dir/all_coding_exons.tsv

... music bmr calc-bmr \
--bam-list input_dir/bam_list \
--maf-file input_dir/myMAF.tsv \
--output-dir output_dir/ \
--reference-sequence input_dir/all_sequences.fa \
--roi-file input_dir/all_coding_exons.tsv \
--genes-to-ignore GENE1,GENE2


bmr-output Number

roi-file Text
Tab delimited list of ROIs [chr start stop gene_name] (See DESCRIPTION)

gene-mr-file Text

reference-sequence Text
Path to reference sequence in FASTA format

bam-list Text
Tab delimited list of BAM files [sample_name normal_bam tumor_bam] (See DESCRIPTION)

output-dir Text
Directory where output files will be written (Use the same one used with calc-covg)

maf-file Text
List of mutations using TCGA MAF specification v2.3


skip-non-coding Boolean
Skip non-coding mutations from the provided MAF file

Default value 'true' if not specified

noskip-non-coding Boolean
Make skip-non-coding 'false'

skip-silent Boolean
Skip silent mutations from the provided MAF file

Default value 'true' if not specified

noskip-silent Boolean
Make skip-silent 'false'

bmr-groups Integer
Number of clusters of samples with comparable BMRs (See DESCRIPTION)

Default value '1' if not specified

show-skipped Boolean
Report each skipped mutation, not just how many

Default value 'false' (--noshow-skipped) if not specified

noshow-skipped Boolean
Make show-skipped 'false'

separate-truncations Boolean
Group truncational mutations as a separate category

Default value 'false' (--noseparate-truncations) if not specified

noseparate-truncations Boolean
Make separate-truncations 'false'

merge-concurrent-muts Boolean
Multiple mutations of a gene in the same sample are treated as 1

Default value 'false' (--nomerge-concurrent-muts) if not specified

nomerge-concurrent-muts Boolean
Make merge-concurrent-muts 'false'

genes-to-ignore Text
Comma-delimited list of genes to ignore for background mutation rates


Given a mutation list (MAF), and per-gene coverage data calculated using "music bmr calc-
covg"), this script calculates overall Background Mutation Rate (BMR) and BMRs in the
categories of AT/CG/CpG Transitions, AT/CG/CpG Transversions, and Indels. An optional
category for truncational mutations can also be specified. The script generates a file
with per-gene mutation rates that can be used with the tool that tests for significantly
mutated genes (music smg).


The regions of interest (ROIs) of each gene are typically regions targeted for
sequencing or are merged exon loci (from multiple transcripts) of genes with 2-bp
flanks (splice junctions). ROIs from the same chromosome must be listed adjacent to
each other in this file. This allows the underlying C-based code to run much more
efficiently and avoid re-counting bases seen in overlapping ROIs (for overall covered
base counts). For per-gene base counts, an overlapping base will be counted each time
it appears in an ROI of the same gene. To avoid this, be sure to merge together
overlapping ROIs of the same gene. BEDtools' mergeBed can help if used per gene.
The reference sequence in FASTA format. If a reference sequence index is not found
next to this file (a .fai file), it will be created.
Provide a file containing sample names and normal/tumor BAM locations for each. Use
the tab- delimited format [sample_name normal_bam tumor_bam] per line. Additional
columns like clinical data are allowed, but ignored. The sample_name must be the same
as the tumor sample names used in the MAF file (16th column, with the header
Ideally, we want to test the mutation rate (MR) of a gene in a sample against the
background mutation rate (BMR) across that sample. But if the BMRs of some samples are
comparable, we can instead test the MR of a gene across a group of samples with
comparable BMR, against the overall BMR of that group. This argument specifies how
many such groups you want to cluster samples into. By default, it is assumed that all
samples have comparable BMRs (bmr-groups = 1).
This should be the same output directory used when running "music bmr calc-covg". The
following outputs of this script will also be created/written: overall_bmrs: File
containing categorized overall background mutation rates. gene_mrs: File containing
categorized per-gene mutation rates.
A comma-delimited list of genes to ignore for overall BMR calculations. List genes
that are known factors in this disease and whose mutations should not be classified as

Use gmt-music-bmr-calc-bmrp online using onworks.net services

Free Servers & Workstations

Download Windows & Linux apps

  • 1
    Strikr Free Software project. Artifacts
    released under a 'intent based'
    dual license: AGPLv3 (community) and
    CC-BY-NC-ND 4.0 international
    Download strikr
  • 3
    giflib is a library for reading and
    writing gif images. It is API and ABI
    compatible with libungif which was in
    wide use while the LZW compression
    algorithm was...
    Download GIFLIB
  • 4
    Alt-F provides a free and open source
    alternative firmware for the DLINK
    DNS-320/320L/321/323/325/327L and
    DNR-322L. Alt-F has Samba and NFS;
    supports ext2/3/4...
    Download Alt-F
  • 5
    Usm is a unified slackware package
    manager that handles automatic
    dependency resolution. It unifies
    various package repositories including
    slackware, slacky, p...
    Download usm
  • 6
    Chart.js is a Javascript library that
    allows designers and developers to draw
    all kinds of charts using the HTML5
    canvas element. Chart js offers a great
    array ...
    Download Chart.js
  • More »

Linux commands

  • 1
    a2j - Wrapper script to simulate
    a2jmidid's non-DBUS behaviour though
    a2jmidid actually being in DBUS mode ...
    Run a2j
  • 2
    a2jmidid - JACK MIDI daemon for ALSA
    MIDI ...
    Run a2jmidid
  • 3
    corebird - Native Gtk+ Twitter client
    for the Linux desktop. ...
    Run corebird
  • 4
    coredumpctl - Retrieve coredumps from
    the journal ...
    Run coredumpctl
  • 5
    ganyremote - GTK+ frontend for
    anyRemote ...
    Run ganyremotex
  • 6
    GAP - Groups, Algorithms and
    Programming DESCRIPTION: GAP is a system
    for computational discrete algebra with
    particular emphasis on computational
    group theory...
    Run gap
  • More »