Subsections

Identify Variants (TAS)

The Identify Variants (TAS) tool takes sequencing reads as input and returns identified variants as part of a Track List.

The tool runs an internal workflow, which starts with mapping the sequencing reads to the human reference sequence, followed by removal of duplicate mapped reads (to reduce biases introduced by target enrichment). The resulting read mapping is analyzed by the Structural Variant Caller to infer indels and other structural variants from unaligned end read patterns. Subsequently, the mapping is realigned, guided by the indels detected by the Structural Variant Caller. The locally realigned read mapping is analyzed by the Low Frequency Variant Detection tool. The Low Frequency Variant Detection tool produces a track of unfiltered variants; these are subjected to a number of post filters to remove variants that are likely due to artifacts or noise. The variants called by the Low Frequency Variant Detection tool that pass the post filtering criteria can be found in the Identified variants track. Variants inferred by the Structural Variant Caller, and not detected by the Low Frequency Variant Detection tool, are also subjected to a number of post filters; those that pass the post filter criteria can be found in the Indels indirect evidence track.

In addition, a targeted region report is created to inspect the overall coverage and mapping specificity in the targeted regions.

Before starting the workflow, you will need to import in the CLC Workbench a file with the genomic regions targeted by the amplicon or hybridization kit. Such a file (a BED or GFF file) is usually available from the vendor of the enrichment kit and sequencing machine. Use the Import | Tracks tool to import it in your Navigation Area.

Run the Identify Variants (TAS) workflow

To run the Identify Variants (TAS) workflow, go to:

        Template Workflows | Biomedical Workflows (Image biomedical_twf_folder_open_16_n_p) | Targeted Amplicon Sequencing (Image targeted_sequencing_closed_16_n_p) | Somatic Cancer (Image somatic_folder_closed_16_n_p) | Identify Variants (TAS) (Image identify_variants_tas_16_n_p)

  1. Select the sequencing reads from the sample that should be analyzed (figure 21.27).

    Image identify_variants_wizardstep1_tas
    Figure 21.27: Please select all sequencing reads from the sample to be analyzed.

    If several samples should be analyzed, the tool has to be run in batch mode. This is done by checking "Batch" and selecting the folder that holds the data you wish to analyze.

  2. In the Target regions dialog you specify the target regions for your application (figure 21.28). The variant calling will be restricted to these regions.

    Image identify_variants_wizardstep2_tas
    Figure 21.28: Select the track with the targeted regions from your experiment.

  3. In the next dialog, you have to select which reference data set should be used to identify variants (figure 21.29).

    Image identify_variants_tas
    Figure 21.29: Choose the relevant reference Data Set to identify variants in your sample.

  4. In the next wizard step (figure 21.30) you can specify the parameters for variant detection.

    Image identify_variants_wizardstep3_tas
    Figure 21.30: Specify the parameters for variant detection.

  5. In the next wizard step (figure 21.31) you specify the parameters for the QC reporting on the targeted regions.

    Image identify_variants_wizardstep4_tas
    Figure 21.31: Specify minimum coverage for the QC reporting on the targeted regions.

  6. In the last wizard step you can check the selected settings by clicking on the button labeled Preview All Parameters. In the Preview All Parameters wizard you can only check the settings, and if you wish to make changes you have to use the Previous button from the wizard to edit parameters in the relevant windows.

  7. Choose to Save your results and click Finish.

Output from the Identify Variants (TAS) workflow

The Identify Variants (TAS) tool produces five different types of output:

It is important that you do not delete any of the produced files individually as some of the outputs are linked to other outputs. If you would like to delete the outputs, please always delete all of them at the same time.

First have a look at the mapping report to see if the coverage is sufficient in regions of interest (e.g. > 30 ). Furthermore, check that at least 90% of reads are mapped to the human reference sequence. In case of a targeted experiment, also check that the majority of reads are mapping to the targeted region.

Afterwards please open the Genome Browser View file (see  21.32).

The Genome Browser View includes the track of identified variants in context to the human reference sequence, genes, transcripts, coding regions, targeted regions and mapped sequencing reads.

Image identify_variants_result1_tas
Figure 21.32: The Genome Browser View allows you to inspect the identified variants in the context of the human genome.

Open the variant track as a table to see information about all identified variants (see 21.33).

Image identify_variants_result2_tas
Figure 21.33: Genome Browser View with an open track table to inspect identified variants more closely in the context of the human genome.