The x axis represents the samples. The y axis represents the exons, one row per exon, and they are sorted by average coverage across samples. For exons with exactly the same average coverage, they are sorted next by the %GC of the exon. (The secondary sort is especially useful for the zero-coverage exons at the bottom). If the figure is unpopulated, then full coverage is assumed (e.g. MutSig CV doesn't use WIGs and assumes full coverage).

Figure 1. 

The mutation spectrum is depicted in the lego plots below in which the 96 possible mutation types are subdivided into six large blocks, color-coded to reflect the base substitution type. Each large block is further subdivided into the 16 possible pairs of 5' and 3' neighbors, as listed in the 4x4 trinucleotide context legend. The height of each block corresponds to the mutation frequency for that kind of mutation (counts of mutations normalized by the base coverage in a given bin). The shape of the spectrum is a signature for dominant mutational mechanisms in different tumor types.

Figure 3.  Get High-res Image SNV Mutation rate lego plot for entire set. Each bin is normalized by base coverage for that bin. Colors represent the six SNV types on the upper right. The three-base context for each mutation is labeled in the 4x4 legend on the lower right. The fractional breakdown of SNV counts is shown in the pie chart on the upper left. If this figure is blank, not enough information was provided in the MAF to generate it.

Figure 4.  Get High-res Image SNV Mutation rate lego plots for 4 slices of mutation allele fraction (0<=AF<0.1, 0.1<=AF<0.25, 0.25<=AF<0.5, & 0.5<=AF) . The color code and three-base context legends are the same as the previous figure. If this figure is blank, not enough information was provided in the MAF to generate it.

Column Descriptions:

• N = number of sequenced bases in this gene across the individual set

• n = number of (nonsilent) mutations in this gene across the individual set

• npat = number of patients (individuals) with at least one nonsilent mutation

• nsite = number of unique sites having a non-silent mutation

• nsil = number of silent mutations in this gene across the individual set

• n1 = number of nonsilent mutations of type: *CpG->(A/T)

• n2 = number of nonsilent mutations of type: *Cp(A/C/T)->(A/T)

• n3 = number of nonsilent mutations of type: A->(C/G)

• n4 = number of nonsilent mutations of type: flip

• n5 = number of nonsilent mutations of type: indel+null

• n6 = number of nonsilent mutations of type: double_null

• p_classic = p-value for the observed amount of nonsilent mutations being elevated in this gene

• p_ns_s = p-value for the observed nonsilent/silent ratio being elevated in this gene

• p_cons = p-value for enrichment of mutations at evolutionarily most-conserved sites in gene

• p_joint = p-value for clustering + conservation

• p = p-value (overall)

• q = q-value, False Discovery Rate (Benjamini-Hochberg procedure)

In this analysis, COSMIC is used as a filter to increase power by restricting the territory of each gene. Cosmic version: v48.

Note:

n - number of (nonsilent) mutations in this gene across the individual set.

cos = number of unique mutated sites in this gene in COSMIC

n_cos = overlap between n and cos.

N_cos = number of individuals times cos.

cos_ev = total evidence: number of reports in COSMIC for mutations seen in this gene.

p = p-value for seeing the observed amount of overlap in this gene)

q = q-value, False Discovery Rate (Benjamini-Hochberg procedure)