This is an old revision of the document!


Analysis of Whole-Exome Sequencing Data for the Diagnosis of Peripheral Neuropathies

1. Methods

1.1 Pipeline's description

pipeline_scheme

1.2 Primary Analysis

  • Read counts from BAM files generated by both pipelines were obtained using FastQC in different quality control steps performed as well as using SAMtools.
  • For variants comparison, VCF files were manipulated using bash scripting to:
    1. Count total number of variants in whole exome and neuropathy-related genes.
    2. Obtain a list of matching genomic positions in which variants had been founded by both pipelines.
  • Statistical metrics reports over variants were generated using VARIANT.

1.3 Coverage Analysis

  • Sequencing coverage was evaluated in final BAM files from which variants were called in both pipelines using GATK DepthOfCoverage.
  • For both pipelines, sequencing coverage was assessed in exome capture kit targets only. BEDtools was used to intersect a list of examined neuropathy-related genes genomic coordinates with a BED file containing exome capture kit targets (SureSelect Human All Exon V5) provided by the manufacturer.
  • R scripts were used to manipulate coverage data generated by GATK DepthOfCoverage and calculate stats. Plots were built with R's package ggplot2

2. Results

2.1 Primary Analysis

2.2 Coverage Analysis

espinos/exome_analysis.1439333563.txt.gz ยท Last modified: 2015/08/12 00:52 by cespinos
www.chimeric.de Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0