Genomatix-Logo
Overview of Help-Pages
GenomeBrowser-Logo

Genomatix Genome Browser


[Introduction] [Getting started] [The Genome Browser Window] [Customizing Genome Browser]

Introduction

The Genomatix Genome Browser is our intuitive tool for visualizing regions of the genome and the associated annotation overlaid with multiple instances of your own data (e.g. from NGS experiments). You can

Whether you want a detailed overview of the various transcripts of a gene, combine ChIP-Seq and RNA-Seq data in genomic context or export a graphic for your next paper, these are all tasks you can achieve with Genome Browser.

Here's an example view:

track/coverage popup

Getting started

Genome Browser is part of the Genomatix Software Suite and can be accessed directly from the navigation bar under 'Genes & Genomes':

accessing GenomeBrowser from the navigation bar

via direct links (e.g. from Gene Overview result):

accessing GenomeBrowser via direct link

The Genome Browser view will start with your default species but with no specific region (you can start to view data by entering a region or gene in the toolbar). Clicking on the organism icon in the top right corner allows you to change the organism or genome build:

The GenomeBrowser view

Coming from ElDorado or Gene Overview pages, the region or gene you specified will be shown:

GenomeBrowser with chosen gene

The GenomeBrowser Window

The browser window consists of four parts - the toolbar on top, the general annotation tracks and the user data tracks in the middle and a scale on the bottom. While the user data is displayed below the annotation data per default you can move tracks around between the two. We will introduce these sections in more detail in the next few paragraphs.

GenomeBrowser window sections

The toolbar

Here is an overview of the toolbar and its elements:

GenomeBrowser toolbar overview

toolbar
Move and zoom: clicking any of the arrows will shift the display view to the left or right of the current view (single arrows: a half window size, double arrows: a full window size). You can also move the view by clicking in the window, holding and moving the mouse. If there are too many tracks to fit them into your browser window, you can also scroll down using this method.

Zooming in or out can be either achieved by clicking on the scaling bars (zoom levels range from 50 bp to 1Mb), by using the '+' and '-' buttons or using your mouse wheel. For the two highest levels the nucleotide sequence will be displayed.


Go back and forward: These buttons allow you to jump back and forward between chromosomal positions viewed before in the current session. The file selection is not affected!

track selection button
Track selection (select annotation and user data): This button allows you to individually select annotation tracks from the ElDorado database, your user data in BED or BAM format, variants in (g)VCF format, public chromatin state data in BED format and other tracks (like the sequence). If you upload or save BED/BAM files in the project management while working with the Genome Browser you can update your BED/BAM file list (for the selected organism) by clicking on the reload button (refresh tracklist button) :
General annotation track selectionBED file track selection
To see the source of the elements annotated in ElDorado just move your pointer over any of the elements or have a look at this ElDorado help page. Probes are only shown for the microarray selected on the ElDorado or the Genome Browser start page.

matrix view button
ENCODE track selection (select public ENCODE data): This button allows you to individually filter and select public ENCODE data in BED format. The dialog offers various filter categories on the left, a matrix representation in the middle, and a classical tree, for track selection, on the right. Filters can also be applied, by clicking on any of the matrix's cells.
MatrixView

show/hide position line button
Show / hide position line: will display a thin vertical position line indicating the exact genomic location of its position. It can also be used to visually align different features.


Select region for ElDorado / BED-file-generation: Click on this button and select a region of interest using your mouse to query the ElDorado database, or to generate a BED file from the selected region.

settings button
Settings: The settings window allows you to change some global parameters of the visualization. You can define the resolution for the coverage tracks and modify the appearance of the track legend. Any settings here affect all tracks shown:

settings window

export button
Export: allows you to export the current view as jpg or png image.


Save:
Opens a panel, which allows you to save your current view state. Reloading a view state will display the chromosomal region with the customized data tracks preloaded. This panel allows you also to save and manage multiple view states. Please note that only view states of the currently selected project will be displayed for loading. Moving the mouse pointer over the species icon at the top right shows the currently selected project. Clicking the icon, allows to change the project.

selection of genes or regions
Selecting a gene or chromosomal regions: Using this input you can jump to a different chromosomal region by either entering the chromosomal coordinates (e.g. 'chr1' or 'NC_00009' plus start and end position), a gene symbol or a gene id. If you enter coordinates, just hit the return key. If you enter a gene symbol, a drop down list will pop up, where you can select your gene from:

gene selection list

In the lower part of the toolbar you can find the chromosomal location bar and an icon for the currently displayed species:

location bar

The yellow indicator shows the relative chromosomal position of the current view. You can also drag it using the mouse to get to another chromosomal position. If cytogenetic G banding pattern information from UCSC is available, then cytobands are drawn. Darker bands are AT-rich, and lighter bands are GC-rich. Moving the mouse pointer over the cytobands will show the cytoband type. 'gneg' and 'gpos' define the stain intensity, 'acen' means acrocentric bands, 'gvar' stands for variable heterochromatic region and 'stalk' refers to tightly constricted regions on the short arms of acrocentric chromosomes. Cytobands are available for Homo sapiens, Mus musculus, Rattus norvegicus and Drosophila melanogaster. Moving the mouse pointer over the species icon will show the species name, ElDorado version, project and the currently selected chromosome.

The general annotation tracks

All tracks in the Genome Browser consist of a track button (a small triangle or square to the very left of the track), a track label indicating the type of annotation data or based on the name of user data and the visualization of the data. If a track can't be shown, e.g. because the selected range is too big to get a meaningful visualization of the data, this will be indicated by a yellow overlay. An animated loading circle indicates that data is still being loaded for that track. You can hide the track buttons and labels or choose to use the track legend window instead via the settings menu as described above. Here are screenshots showing the track buttons within the tracks and the track legend window:

track labelstrack legend window

Annotation tracks

All annotation data from Genomatix' ElDorado database can be displayed in this section. The tracks selected by default are microRNAs (from miRBase, pink), modules (known functional transcription factor binding site combinations, dark blue), promoter regions (yellow), repeat regions (light orange), transcription start sites (red) and transcript (black). If you started the Genome Browser via ElDorado, the user sequence track (dark orange) will also be shown, marking the location of the gene you entered or the aligned sequence if you submitted a sequence or accession number.

To add or remove tracks use the track selection button (track selection button) in the toolbar. (Please note that the 'transcripts' track is located in the 'other tracks' section.)

Clicking on a track button (the small triangles on the left side of the tracks: ) will expand the associated track. For most tracks this will result in plus and minus strand being displayed separately, but for the transcript track all transcript variants are shown. Strand orientation for these is indicated by the small grey arrows next to the transcript accession numbers. Tracks with a right pointing triangle can't be expanded. Clicking on a track button of an expanded annotation track () will show the annotations stacked among each other. This screenshot shows the promoter regions in an unexpanded track, the repeats in a strand-specific track and the modules (track with the button) are stacked on each other:

Transcript track

The unexpanded transcript track shows all transcripts in one row. The thin lines denote the introns and the thicker lines denote the exons. The thinner part of the exon is the UTR and the thicker part is the coding sequence. Non-coding transcripts are represented by white lines with black borders. The gene symbol for each locus is displayed as long as it does not overlap with the previous gene symbol.

The expanded transcript track shows each transcript stacked above each other. Transcripts from different loci are drawn with more space between them.

If you move your mouse pointer over a transcript, additional information is shown in a popup window:

transcript info popup

Furthermore, right-clicking the transcript opens a popup menu. Via the 'Links' option you can e.g. go to ElDorado's 'More Gene Info' or 'Transcript Details' pages for more detailed information on the chosen transcript:

element links popup

If you want to display a subset of transcripts, use the "Remove genes and transcripts" option from the context menu, to choose them individually.

Enhancer track

Right-clicking an enhancer, allows to jump to the locus of a linked gene via the 'Genes' menu, if available.

enhancer gene link

The 'Sequence' track

The sequence track can be used to display the underlying reference sequence of the region you're looking at. If you zoom in to the 250 bp level (second smallest scaling bar) the sequence will be displayed as color code:

sequence track colored nucleotides

If you zoom in further to the 50 bp zoom level (smallest scaling bar) the nucleotides will be displayed:

sequence track full nucleotides

For zoom levels where displaying either the sequence or the color coded nucleotides isn't possible, the sequence track will just show checkered grey and yellow squares.

The 'Position' track

The position track shows the exact chromosomal position of the region currently displayed.

The scale

The scale below the tracks helps you to quickly verify the approximate genomic distances you're looking at in the current view. It can be hidden via the settings menu.

The user data tracks

All BED and BAM files stored in the project management for the currently selected organism and genome version can be added to the user data section. To add data, use the track selection button () in the toolbar, then pick the 'BED file tracks' or 'BAM file tracks' option, open the desired project and tick those files that you want to add as coverage or as alignment track:

adding tracks

Coverage tracks

By default the coverage tracks are displayed as bar charts in black and with auto scaling. Each bar represents a distinct region of the underlying sequence, which we call 'bin'. For each bin a coverage value is calculated by dividing the number of 'covered' nucleotides by the total number of nucleotides in the bin. Multiple coverage is taken into account.

Example

                GATC   (region 1 - partial)
  ACCCTCCACA           (region 2)
    CCTCCACACTC        (region 3)
AGACCCTCCACACTCTGATC   (bin)

Here the bin is 20 nucleotides long. It is covered by 3 regions from the BED file, with a coverage value of (4 + 10 + 10)/20 = 1.2

Please note that the bins are recalculated dynamically whenever you change the zoom level or move the view. You can change the number of bins displayed via the 'Settings' menu.

Moving your mouse pointer over a bar will give you the name of the BED file associated with the track and coverage for that bin:

track/coverage popup

You can customize each user data track by setting different colors or choosing another chart type. Combining tracks is also possible. To learn more, please refer to the section 'Customizing Genome Browser'.

Coverage tracks with scores

The regions in BED files can contain scores in the fifth column. The official format allows a value between 0 and 1,000. The Genomatix Genome Browser does also handle values greater than 1,000 as well as negative values. The scores can be visualized in the coverage track by right clicking the track, selecting coverage settings and then selecting either Visualize positive scores or Visualize positive and negative scores. Please note that visualization of positive and negative scores in combination with a strand-specific view is not possible for the same coverage track.

The coverage of each bin is calculated by adding the products of region length and score and dividing it through the length of the bin.

Example

                GATC   (region 1 / score: 0.5)
  ACCCTCCACA           (region 2 / score: 1)
    CCTCCACACTC        (region 3 / score: 2)
AGACCCTCCACACTCTGATC   (bin)

Here the bin is 20 nucleotides long. It is covered by 3 regions from the BED file, with a coverage value of (4*0.5 + 10*1 + 10*2)/20 = 1.6

The next screenshot shows an alignment track and coverage track of a BED file with the visualization of positive scores. The outter left/right bins have a low score and the middle bins have a high score.

coverage score tracks

Strand-specific coverage tracks

Coverage-tracks can be expanded like most annotation tracks by clicking the triangle button (). Then one coverage plot is shown for the plus and one for the minus strand. The bin range on the right still displays the highest value for both strands together. The next screenshot displays two strand-specific coverage tracks. The first coverage track shows strand-unspecific RNA-seq data and the second coverage track shows strand-specific data:

Strand-specific coverage tracks with the plot type 'Region map' are displayed like strand-specific annotation tracks.

Alignment tracks

Alignment tracks display the segments of your BED files or aligned reads from your BAM files. These tracks are shown up to a view range of 25,000 bp, otherwse you have to zoom in.

If you zoom close enough to a segment/read, then the strand-orientation is displayed as arrowhead:

BAM files

Spliced reads and deletions are visualized with a thin line between the aligned read parts. At higher resolutions read bases that mismatch the reference are visualized with the same color code as in the sequence track (). Insertions are visualized with a yellow bar ().

If you move your mouse pointer over a read, additional information like the name, quality and the sequence is shown in a popup window. Mismatches are highlighted in lower-case and insertions are marked bold.

(g)VCF files

SNPs and indels are visualized by red / blue colored rectangles, while the portion of red indicates the alternative allele frequency. The rectangles are aligned where bases are actually modified, so that for e.g. deletions, the length and position of the rectangles correspond to the actual bases that are deleted.
Samples are shown on the extended track and are colored as follows:

hom refhomhetno call
whiteblackblack/whitegrey

Insertions have a little triangle cut out and are positioned between the two bases of insertion.

Looking at gVCFs, non-variant blocks are exclusively displayed on the extended track, as white blocks with black outline. However, additional corresponding tooltips are displayed on the main track.

Public data

Under 'BED files' you can find the group 'public' with more than 1,500 BED files. These BED files are available for every user.

Chromatin states

These BED files contain chromatin states of nine cell lines according to Jason Ernst et al., Nature 473: 43-49 (2011). Ernst et al. identified different chromatin states corresponding to regulatory elements such as active and inactive promoters or weakly and strongly transcribed regions.

In the next screenshot the regulatory elements for the cell line Gm12878 are displayed. The BED files have been added as coverage tracks with the plot type 'Region map'. More information about plot types can be found in the chapter 'Customizing Genome Browser' below.

ENCODE

The data of these tracks are from the ENCODE project (ENCODE Project Consortium. "A user's guide to the encyclopedia of DNA elements (ENCODE)." PLoS Biol 9.4 (2011): e1001046.). The tracks are divided into the categories DNase hypersensitivity site, Histone modification, TF binding site, Polymerase binding site and Other factor binding site.

The tracks in DNase hypersensitivity site show regions where the chromatin is hypersensitive to cutting by the DNase enzyme. The tracks can be loaded individually for each cell type or as a supertrack consisting of DNase hypersensitivity sites from all cell types. For supertracks, the number below the alignment track shows how many cell lines are hypersensitive in the region. At the functional level hypersensitive sites often represent transcriptional active regions.

The tracks in Histone modification show regions with modifications of the histones H3k27ac, H3K4me1-3 and H3K9ac. The tracks can also be loaded individually for each cell type or as histone-specific supertracks combined of all cell types. At the functional level the histone peaks are associated with transcriptional activity.

The tracks in the category TF binding site show ChIP-Seq peaks of more than 150 TFs. The data can be loaded in four different ways:

  1. single TF in a specific cell type
  2. all TFs in a specific tissue (supertrack)
  3. single TF in all cell-types/tissues (supertrack)
  4. all TFs in all cell types/tissues (supertrack)

The supertracks with more than one TF have a special feature. Right clicking on an alignment track gives you the option to either load a supertrack for a TF or a cell type-specific TF track. Right clicking on the peak of a TF, filters the options to the according TF. The next screenshot shows the option menu when right clicking on the peak region of GATA2:

TF option menu

The score of each region represents the signal strength of its peak. The score of the regions in the supertracks is the sum of the signal strength in the individual cell types.

The tracks in the categories Polymerase binding site and Other factor binding site show ChIP-Seq peaks of Polymerase and other factors binding to DNA that are not characterized as transcription factors. Loading and supertrack handling for these is the same as for the Transcription factor binding sites.

Evolutionary constrained elements

This BED file contains evolutionary constrained elements according to Kerstin Lindblad-Toh et al., Nature 478: 476-482 (2011). Kerstin Lindblad-Toh et al. compared the human genome sequence with 29 mammalian genome sequence and identified 4.2% of the human genome as evolutionary constrained.


Customizing Genome Browser

Adding and removing tracks

Adding tracks

To add tracks to either the annotation or user data area use the track selection button () as described in the sections above. You can either select each track separately by ticking the checkboxes or you can use the button () to the right of the track groups. Clicking the button shows several options to add tracks from the group. The different options are only shown if no more than 20 tracks would be added.


Please note that if you add more tracks than can be shown in your browser window, you will have to drag up, to view those tracks.

Removing tracks

To remove a track, either click on the track selection button () and deselect the track in the popup-menu or right click the track and select the 'Remove track' option from the pop-up menu. To remove several tracks from one group at once, click in the track selection popup-menu the button () and select "Remove all tracks".

The first number in parentheses behind the group name indicates the number of selected tracks and the second number in parentheses indicates the total number of tracks.

Changing the height of a track

You can change the height of any track in both the annotation and user data area. Right click on the track and select 'Set track height' from the menu. In the upcoming window you can set the track height in pixels. To change the height for several tracks at once, check the box 'Apply to multiple track' and select the corresponding tracks:

changing track heights

Changing the color of a track

You can also change the color of a track. Right click it, select 'Set color' and then choose a color using the color picker in the upcoming window. You can also change the transparency of a track if the track is stackable.To change the color and transparency for several tracks at once, check the box 'Apply to multiple tracks' and select the corresponding tracks:

setting track colors

Changing the scale of coverage tracks

Right-clicking a user data track will bring up an extended popup menu, where you can choose from different ways to scale the y-axis of a track. The default is 'Auto scale':

setting track scales

With 'Auto scale', the data are automatically scaled, which means that each track is scaled based on its maximum value in the currently displayed region. As a consequence, the maximum values displayed differ between tracks and will also differ if you move to a different location of the chromosome. The displayed data range for each track is shown on the right hand side of the track.

If you select 'Local common auto scale', all tracks for which this scaling method was selected will be scaled based on their common maximum in the displayed region. This can be helpful to assess the differences between data in the region you're looking at.

If you select 'Global common auto scale', each lane is in addition normalized based on the total length of regions and thus taking the numbers and lengths of all bed files into account where this scaling method was selected. This will give you more of a global view of your data.

Please note that if you select one of the common scales, data in your current view might be scaled down to being invisible if the maximal values in other tracks are high in comparison. The best way to check for the presence of any data at the current position is to use 'Auto scale', for comparing tracks either the 'local' or 'global common scales' are recommended.

'Manual scale' lets you enter a minimum and maximum value for the y-axis range. Depending on the values, some data might become invisible due to the scaling!

Selecting 'Log scale' will change the scale from linear to log 2 based. This can be applied to any of the 4 scaling choices.

The option 'Instant scale' is enabled by default. If you move the view by clicking in the window, holding the mouse button and moving the mouse, then the tracks are scaled instantly. If the option is disabled, then the view does not scale until you release the mouse button. If you added many data tracks and cannot move the view smoothly, it is recommended to disable this option.

The selected scaling options are displayed if you move your cursor over the bin range on the right side, and are also indicated by a small symbol beside the bin range:

Example

For the view below, tracks for DUX4 expression (control) were colored in blue, while the induced replicates were colored in red. The scale type was set to 'Local common auto scale' for all:

example

Changing the chart type of coverage tracks

Besides changing the height or color of a track you can also choose from four different chart types for user data tracks:

track types

'Bar chart' is the default, where each bar stands for a bin containing regions from the BED file. As in a histogram, the more regions from the BED file fall into a bin the higher the bar will be plotted.

'Line plot' uses the same binning approach as 'bar chart' but will display a line instead of bars. You could think of it as the outline of a bar chart.

'Heatmap' will display the values of the bin in a heatmap style, i.e. the value of each bin is indicated by the saturation of the drawn bar. Higher values will have 'darker' bars, lower values are 'lighter'. All bars are drawn in the same height.

'Region map' will just indicate that some data has been found for this position in the BED file without assigning any value to it.

The resolution of the user tracks, i.e. the number of bars or line segments drawn can be changed via the 'Settings' button (settings button) in the toolbar:

settings window

Please note that choosing a high number of bins might slow down the responsiveness of the Genome Browser when displaying many user tracks, as more calculations have to be performed!

Example

Here is an example combining several settings in one view: It shows the four different chart types:

all of them set to a resolution of 200 bins:

example

Changing the order of tracks

If you need a different order of tracks, e.g. to better correlate a user track with an annotation track, you can simply move tracks around by clicking on their name and then moving them to the new position and then releasing the mouse button. A line will appear in between other tracks to indicate the position where the track will be inserted:

changing the track order

You can move any track to any position in the display.

Locking tracks

The track view is divided into two containers, which can be independently dragged vertically. You can change the height/ratio of the containers, by dragging the horizontal divider. Double-clicking the divider, resizes the top container to show all tracks in the container at once, if possible. All tracks can be reordered between the two containers and the containers will auto scroll if you drag a track near to the container's top/bottom edges. With that you can "lock" one or multiple tracks in one of the containers, while other tracks can be dragged vertically in the other container without moving the "locked" tracks out of view. The default added tracks are automatically added to the top(1) container on startup, while all further added tracks will be added to the bottom(2) container.

locking tracks

Superimposing tracks

A powerful feature of Genome Browser is the capacity to superimpose tracks. This enables you to convey multiple levels of information in a single track. You can combine several annotation tracks or coverage tracks. If you prefer a superimposed display for the annotation track (as in the old 'Detailed graphics' in ElDorado), e.g. for primary transcripts, exons and promoters, you can overlay these tracks using simple drag and drop (select a track by clicking on the track label, hold and move the track over the other track and release the mouse button when the small plus sign appears):

superimposing tracks

The view below was generated by adding the 'Primary Transcript' and 'Exon tracks' and then combining those (using drag and drop) with the 'Promoter' track:

superimposed tracks

To separate any track from a combined track, simply click on its label and drag it to a different order in the display.

Please note that the 'Sequence', 'Position' and 'Transcript' and alignment tracks can't be superimposed on any other tracks.

You can also superimpose coverage tracks using the same drag & drop method as for the annotation tracks. In combination with different chart types and/or scalings this can be used to create highly informative combined tracks.

Example

For the view shown below two coverage tracks for DUX4 ChIP-Seq and ChIP-Seq-Peaks were added. The ChIP-Seq track's plot-type was set to 'Region Map' and the color was set to red with 0.5 transparency, while the ChIP-Seq-Peak track was set to blue with 0.5 transparency. The resolution was increased and subsequently the ChIP-Seq-Peak track was superimposed on top of the corresponding ChIP-Seq track. You can now see reads and peaks in a single track:

example

Track legend

The track legend is by default integrated in the tracks. You can change the display of the track legend in the settings panel by clicking on the button.

Integrated track legend

The integrated track legend shows a button and a label on the left side for each track. If you have superimposed several tracks, then several button and labels are displayed beneath each other on the superimposed tracks.

Track legend panel

The track legend panel is displayed on top of the tracks, but can be moved freely. Tracks can also be moved by dragging the label between the tracks and can be superimposed by dragging the label on another label. An advantage of the track legend panel is that the superimposed tracks have no minimum height.

Track legend on the left side

The track legend on the left side is displayed separately on the left side. Dragging the labels can also be used to move and superimpose the tracks. If you double-click on the border separating the track legend and the tracks, then the track legend fades out to the left. Double-clicking the separator again, fades in the track legend. You can also resize the left track legend by dragging the separator with the mouse.

Transcript track

Highlight transcripts by source

Right-clicking the transcript track and selecting 'Highlight transcripts by source' opens a small window where you can individually select and deselect transcript sources. Clicking the OK or Apply button, transcripts of deselected sources will be shown greyed out.

highlight sources popup

Remove genes and transcripts

Right-clicking the transcript track and selecting 'Remove genes and transcripts' opens a small window where you can individually select and deselect genes and transcripts shown in the current window. Deselected genes and transcripts will be removed from the track. Clicking on the arrow (Transcript info icon) beside the genes and transcripts will open a new browser tab with more detailed gene or transcript information.

transcript selection popup

Alignment track

Alignment settings

Right-clicking the alignment track and selecting 'Alignment settings' opens a window with several options to customize the alignment track for BED and BAM files.

Alignment settings

Max. coverage depth determines the maximal coverage in the alignment track and can be set to maximal 1000. A maximal coverage of 250 means that no more than 250 alignments are stacked at a single base position.

Enable max. alignment length determines if the alignments will be filtered according to the next option.

Max. alignment length filters out all alignments with a higher length.

Compact view draws the alignments more compact.

Show labels shows the read identifiers as labels..

Strand-specific coloring draws all alignments on the minus strand with an extra color. This color can be set under Minus strand color.

The upper track shows the compact view with strand-specific coloring and the lower tracks shows the default view:

Alignment settings example

SAM/BAM settings

Right-clicking the alignment track and selecting 'SAM/BAM settings' opens a window with several options to customize the alignment track for SAM/BAM files.

SAM/BAM settings

Selecting the option Draw mate pairs connects paired-end reads with a line and draws them side by side if both reads have been loaded. If the mate read is located on another chromosome, the read is highlighted with a colored border. The color for such an inter-chromosomal pair can also be customized. A useful feature for an inter-chromosomal pair is to right-click the mate and then select Jump to mate in order to jump to the location of the according mate. In this example the inter-chromosomal read pairs are highlighted with a red border:

Interchromosomal rearrangement

Selecting the option Highlight pair orientation highlights all read pairs deviating from the expected read pair orientation. The reads are highlighted with three different colors depending on whether the first, the second or both reads are incorrectly orientated. This option allows to detect structural variants more easily. In this example the first and the second reads are incorrectly orientated which indicates a duplication:

Paired end duplication example

In the next example only the second read is incorrectly orientated which indicates an inversion:

Paired end nversion example

Selecting the option Highlight insert size deviation colors all reads deviating from the expected insert size. The expected insert size range can be customized as well as the colors for inferred insert sizes smaller and larger than the expected insert size. The expected insert size range can also be estimated from the insert sizes of the currently loaded read pairs. A more precise range over all read pairs in the BAM file is calculated by the Genomatix Mapper. In this example the read pairs with a smaller insert size than expected are highlighted with a green border and indicate a potential insertion:

Insertion example

In this example the read pairs with a larger insert size than expected are highlighted with a blue border and indicate a potential deletion:

Deletion example

In order not to miss the highlighted reads in data with high coverage, the display option Show all sorted sorts the reads so that the highlighted reads appear at top. The option Show only highlighted displays only the highlighted reads. Please note that the sorting and filtering function will be only applied on the loaded reads. If your coverage depth exceeds the configured limit, please adjust the max. coverage depth (see Alignment settings).

The option Color mapping quality draws the color intensity of the reads depending on the mapping quality value. The higher the mapping quality the more intense the color. The mapping quality range is set by default from 0 to 60. Reads with mapping quality values smaller than the lowest range value are assigned the lowest color intensity and reads with mapping quality values higher than the highest range the highest color intensity. The color intensity for reads with values between the lowest and highest value is scaled linearly. In this example all reads on the left have a lower mapping quality than the reads on the right:

Mapping quality example


Extracting information

There are several ways to extract information about the region, annotation elements and sequences. E.g. you can select a region as mentioned above by clicking on the 'Select region button' () and selecting a region of interest using your mouse to query the ElDorado database, or to generate a BED file from the selected region. ElDorado gives you more information about the selected annotations and has several export functions to save the annotations as BED file or to extract their sequences.

Generate BED file from annotation

Besides generating a BED file from a specific region, BED files can also be generated from annotations. Open the context menu by right-clicking the annotation and choose the annotation via "Generate BED file".

Sequences

Sequences can be obtained from annotations and alignments. You can copy the sequence either to your clipboard or on the server and save it in the result management.

Annotation sequences

Right-clicking on annotation opens a context menu from where you can choose the annotation sequence to display. You can add the flanking sequences on both sides of the annotation which might be useful for SNP annotation.

Sequence window

Read sequences

Read sequences can also be displayed by right-clicking on the respective read in the alignment track. For BAM files the read sequence will be displayed with all mismatches, deletions and insertions in the reference (if the MD string is available). Mismatches are highlighted in lower-case. For BED files the reference sequence will be shown.

Open read sequenceRead sequence window

You can copy the sequence either as plain text or in FASTA format to your clipboard or to the server where you can save the sequence in the result management.

Links

More detailed information can be obtained for most annotation elements: exon, microRNA, module, primary transcript, promoter region, transcript, transcript start region, UTR and variation. Right-clicking on an annotation element opens the context menu with links to analyses or more details of the annotation element. E.g. you can start the comparative genomics task for a promoter region or analyze it for binding sites: