背景

作者是一个没有生信背景知识的程序员,近期对全外显子测序有些感兴趣,因此记录一下学习过程,和大家共同进步。

预备知识

WES基础知识
什么是基因
什么是DNA
二代测序基础知识
等位基因

软件安装

bwa

git clone https://github.com/lh3/bwa.git
cd bwa
make

由于bwa的makefile没有写intall的脚本,所以我们直接将编译的目录加入到环境变量中,

vim ~/.bashrc
export PATH=${PATH}:/usr/local/cuda/bin/:/home/fangl5/Downloads/bwa-0.7.17/

samtools

下载samtools-1.14.tar.bz2
解压后

    ./configuremakesudo make install

java 11

sudo apt-get install openjdk-11-jdk

picard

wget https://github.com/broadinstitute/picard/releases/download/2.26.10/picard.jar

gatk

sudo apt-get install git-lfs
git clone https://github.com/broadinstitute/gatk.git
./gradlew bundle

将gatk加入到环境变量PATH中,和bwa中操作一致。

SOAPnuke

git clone https://github.com/BGI-flexlab/SOAPnuke.git
cd SOAPnuke
make

加入环境变量中

下载数据

hg38.ga

wget http://hgdownload.cse.ucsc.edu/goldenPath/hg38/bigZips/hg38.fa.gz

gatk参考文件

我直接将google cloud中文件夹下载下来了

gsutil -m cp -r \"gs://genomics-public-data/resources/broad/hg38/v0/" \.

gatk-tutotial 相关数据 google cloud地址

genomics-public-data

annovar 数据下载

http://hgdownload.cse.ucsc.edu/goldenpath/hg38/database/

实践fastq to vcf

# 大于2G的参考基因组要用参数-a bwtsw
# 约40min
bwa index -a bwtsw hg38.fa
java -jar ~/Downloads/picard.jar CreateSequenceDictionary R=./hg38.fa O=./hg38.dict
samtools faidx hg38.fa
gatk IndexFeatureFile -I ./Mills_and_1000G_gold_standard.indels.hg38.vcf

https://app.terra.bio/#workspaces/help-gatk/GATK4-Germline-Preprocessing-VariantCalling-JointCalling/notebooks/launch/1-gatk-germline-variant-discovery-tutorial.ipynb

Q&A

核苷酸序列或氨基酸序列是什么?

WGS\WES\TS区别

http://www.genomesop.com/wp-content/uploads/2018/11/genomic_comparison3.png

参考资料

GATK4 流程分析- 从fastq到vcf
GATK_Discovery_Tutorial-Worksheet-AUS2016.pdf
基因组测序、外显子测序和靶向测序有什么样的区别,如何选择?
华大外显子测序
Exome sequencing data analysis for diagnosing a genetic disease
全外显子测序(WES)分析流程
使用fastp进行数据质控
GATK4 多个样本GenotypeGVCFs前用 CombineGVCFs还是GenomicsDBImport
intel 2015 面向 GATK 最佳实践流水线部署的基础设施
家系分析软件汇总
WGS数据分析入门篇 碱基矿工
全外显子测序(WES)分析流程 这一篇,用到了bed文件。
gatk spark加速

Whole-Genome-Sequencing
全外显子测序(wes)数据分析详细流程(小样本)

全外显子测序(WES)由入门到精通记录相关推荐

  1. 史上最全的MonkeyRunner自动化测试从入门到精通(1)

    史上最全的MonkeyRunner自动化测试从入门到精通(1) 史上最全的MonkeyRunner自动化测试从入门到精通(2) 史上最全的MonkeyRunner自动化测试从入门到精通(3) 史上最全 ...

  2. Axure RP原型设计|全实战项目演练从入门到精通-姜浩-专题视频课程

    Axure RP原型设计|全实战项目演练从入门到精通-124人已学习 课程介绍         多项目实战,讲师手把手详细演示教学 从基础到精通的深入化进阶学习 做优秀的原型设计 课程收益     交 ...

  3. 最全Hadoop视频教程(从入门到精通 视频教程下载)Hadoop八天完全攻克Hadoop视频教程 Hadoop开发新版Hadoop视频教程

    最新在学习hadoop .storm大数据相关技术,发现网上hadoop .storm 相关学习视频少之又少,更可恶的是不少机构使用段海涛老师Hadoop视频进行商业买卖,气愤至急,我这里整理了传智播 ...

  4. 全了!从Python入门到精通

    总被读者问到,我看完了python入门的书,后面就不知道要学什么了.今天就给你们整理全套入门到进阶的教程. 这套教程非常全面而且详细,从Python入门到Python进阶.Django.Flask等W ...

  5. Unity3D游戏美术全攻略:从入门到精通

    文/拉撒路 Unity现在已经用的很广泛啦,可是却一直没有什么美术向的教程. 程序用方面的内容在各个论坛都有讨论,但是美术似乎很弱势啊. 明明美术也很需要掌握引擎方面的内容嘛! 山谷里的野百合还有春天 ...

  6. 全网最全的Apollo教程从入门到精通

    第1章 Apollo项目实战 能力目标 掌握Apollo安装 掌握Apollo配置操作 能够实现SpringBoot集成Apollo 能够实现Apollo适用功能[自动刷新.监听器.配置加密.1灰度发 ...

  7. 【Github教程】史上最全github用法:github入门到精通

    原文 http://www.eoeandroid.com/thread-274556-1-1.html [初识Github] 首先让我们大家一起喊一句"Hello Github". ...

  8. 史上最全的MonkeyRunner自动化测试从入门到精通(10)

    三.MonkeyRunner复杂的功能开始学习 (1)获取APK文件中ID的两种方式 Monkeyrunner的环境已经搭建完成,现在对Monkeyrunner做一个简介. Monkeyrunner工 ...

  9. python编程从入门到精通pdf-码哥—-Python从入门到精通(全60集),资源教程下载...

    课程名称 码哥--python从入门到精通(全60集),资源教程下载 课程目录 1Python编程语言历史及特性 2Python编程语言初接触 3Python程序文件结构 4准备Python编程环境 ...

  10. python从入门到精通视频(全60集)-【网盘下载】Python从入门到精通视频(全60集)...

    Python从入门到精通视频(全60集) JAVA? 今天 课程目录 python入门教程-1-Python编程语言历史及特性 python入门教程-2-Python编程语言初接触 python入门教 ...

最新文章

  1. 创建细分客户的无监督学习项目
  2. 拼多多程序员:因为看到同事被抬上救护车,我被开除了!
  3. hdu1074 状态压缩dp+记录方案
  4. 金庸小说人物关系表年表
  5. 使用StackBlitz和SAP Spartacus快速创建电商店铺页面
  6. PKU 学生反馈 3
  7. 运维前线:一线运维专家的运维方法、技巧与实践3.1 数据中心搬迁准备
  8. Ural_1586. Threeprime Numbers (DP)
  9. wxpython菜单的位置_wx.grid 怎么加菜单呀!
  10. Membership学习(三)Membership Providers介绍
  11. 【渝粤教育】国家开放大学2018年春季 0266-22T设计构成 参考试题
  12. [2018.10.17 T1] 斜率
  13. 操作系统ucore lab1
  14. JS+html+css实现的一个小小的贪吃蛇游戏
  15. 学计算机跨考航天航空,往届生跨考计算机经验
  16. 如何做好一次视频号直播
  17. 暗影骑士vn7安装linux,宏碁暗影骑士vn7 591g笔记本怎么样?
  18. RTMPLive多流媒体协议转无插件直播协议(RTMP)解决方案
  19. 【转】常用的隧道技术
  20. 无论LCD和OLED显示技术有多好,都无法替代这个古老的显示数码管

热门文章

  1. 8个高质量免抠素材网站
  2. (一)目标检测经典模型回顾
  3. 深度学习声纹识别_一种基于机器学习及深度学习的声纹降噪方法及系统与流程...
  4. Microsoft Word 教程:如何在 Word 中插入图片、图标?
  5. boost电路输出电流公式_BOOST电路参数计算公式
  6. 当今主流软件产品家族及其开发语言
  7. 广告数据定量分析:第一章——广告优化中的统计学
  8. 神州数码交换机配置基本命令
  9. 记一次https的免费ssl证书安装及配置过程--window系统
  10. android 时钟动态图标,神奇的Android动态时钟/时间UI界面