Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
espinos:exome_analysis [2015/08/12 00:46] cespinos [1.3 Coverage Analysis] |
espinos:exome_analysis [2016/06/24 15:51] (current) fgarcia [1.1 Pipeline's description] |
||
|---|---|---|---|
| Line 1: | Line 1: | ||
| - | ====== Analysis of Whole-Exome Sequencing Data for the Diagnosis of Peripheral Neuropathies ====== | + | ====== Analysis of Whole-Exome Sequencing data for diagnosis of Inherited Peripheral Neuropathies ====== |
| + | |||
| + | |||
| ===== 1. Methods ===== | ===== 1. Methods ===== | ||
| Line 6: | Line 9: | ||
| ==== 1.1 Pipeline's description ==== | ==== 1.1 Pipeline's description ==== | ||
| - | {{:espinos:figura_pipeline.png?680|pipeline_scheme}} | + | {{:espinos:images_ana:figura_pipeline.png?680|pipeline_scheme}} |
| ==== 1.2 Primary Analysis ==== | ==== 1.2 Primary Analysis ==== | ||
| Line 12: | Line 15: | ||
| * Read counts from BAM files generated by both pipelines were obtained using [[http://www.bioinformatics.babraham.ac.uk/projects/fastqc/|FastQC]] in different quality control steps performed as well as using [[http://www.htslib.org/|SAMtools]]. | * Read counts from BAM files generated by both pipelines were obtained using [[http://www.bioinformatics.babraham.ac.uk/projects/fastqc/|FastQC]] in different quality control steps performed as well as using [[http://www.htslib.org/|SAMtools]]. | ||
| * For variants comparison, VCF files were manipulated using bash scripting to: | * For variants comparison, VCF files were manipulated using bash scripting to: | ||
| - | - Count total number of variants in whole exome and neuropathy-related genes. | + | - Count total number of variants in whole exome and in neuropathies-related genes. |
| - | - Obtain a list of matching genomic positions in which variants had been founded by both pipelines. | + | - Obtain a list of matching genomic positions in which variants had been found by both pipelines. |
| * Statistical metrics reports over variants were generated using [[http://wiki.opencb.org/projects/hpg/doku.php?id=variant:overview|VARIANT]]. | * Statistical metrics reports over variants were generated using [[http://wiki.opencb.org/projects/hpg/doku.php?id=variant:overview|VARIANT]]. | ||
| Line 20: | Line 23: | ||
| * Sequencing coverage was evaluated in final BAM files from which variants were called in both pipelines using [[https://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_gatk_tools_walkers_coverage_DepthOfCoverage.php|GATK DepthOfCoverage]]. | * Sequencing coverage was evaluated in final BAM files from which variants were called in both pipelines using [[https://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_gatk_tools_walkers_coverage_DepthOfCoverage.php|GATK DepthOfCoverage]]. | ||
| - | * For both pipelines, sequencing coverage was assessed in exome capture kit targets only. [[http://bedtools.readthedocs.org/en/latest/|BEDtools]] was used to intersect a list of examinated neuropathy-related genes genomic coordinates with a BED file containing exome capture kit targets provided by the manufacturer (SureSelect Human All Exon V5). | + | * For both pipelines, sequencing coverage was assessed in exome capture kit targets only. [[http://bedtools.readthedocs.org/en/latest/|BEDtools]] was used to intersect a list of significant neuropathies-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 [[https://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_gatk_tools_walkers_coverage_DepthOfCoverage.php|GATK DepthOfCoverage]] to calculate stats. Graphics were built with R's package [[http://docs.ggplot2.org/current/|ggplot2]] | + | * R scripts were used to manipulate coverage data generated by [[https://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_gatk_tools_walkers_coverage_DepthOfCoverage.php|GATK DepthOfCoverage]] to calculate stats and build box and bar plots with [[http://docs.ggplot2.org/current/|ggplot2]] package. |
| - | ===== 2. Results ===== | + | |
| + | ===== 2. Results ===== | ||
| ==== 2.1 Primary Analysis ==== | ==== 2.1 Primary Analysis ==== | ||
| - | |||
| * **[[Results.Primary.Analysis|Sequence processing]]** | * **[[Results.Primary.Analysis|Sequence processing]]** | ||
