全外显子测序(WES)由入门到精通记录
背景
作者是一个没有生信背景知识的程序员,近期对全外显子测序有些感兴趣,因此记录一下学习过程,和大家共同进步。
预备知识
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
./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
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)由入门到精通记录相关推荐
- 史上最全的MonkeyRunner自动化测试从入门到精通(1)
史上最全的MonkeyRunner自动化测试从入门到精通(1) 史上最全的MonkeyRunner自动化测试从入门到精通(2) 史上最全的MonkeyRunner自动化测试从入门到精通(3) 史上最全 ...
- Axure RP原型设计|全实战项目演练从入门到精通-姜浩-专题视频课程
Axure RP原型设计|全实战项目演练从入门到精通-124人已学习 课程介绍 多项目实战,讲师手把手详细演示教学 从基础到精通的深入化进阶学习 做优秀的原型设计 课程收益 交 ...
- 最全Hadoop视频教程(从入门到精通 视频教程下载)Hadoop八天完全攻克Hadoop视频教程 Hadoop开发新版Hadoop视频教程
最新在学习hadoop .storm大数据相关技术,发现网上hadoop .storm 相关学习视频少之又少,更可恶的是不少机构使用段海涛老师Hadoop视频进行商业买卖,气愤至急,我这里整理了传智播 ...
- 全了!从Python入门到精通
总被读者问到,我看完了python入门的书,后面就不知道要学什么了.今天就给你们整理全套入门到进阶的教程. 这套教程非常全面而且详细,从Python入门到Python进阶.Django.Flask等W ...
- Unity3D游戏美术全攻略:从入门到精通
文/拉撒路 Unity现在已经用的很广泛啦,可是却一直没有什么美术向的教程. 程序用方面的内容在各个论坛都有讨论,但是美术似乎很弱势啊. 明明美术也很需要掌握引擎方面的内容嘛! 山谷里的野百合还有春天 ...
- 全网最全的Apollo教程从入门到精通
第1章 Apollo项目实战 能力目标 掌握Apollo安装 掌握Apollo配置操作 能够实现SpringBoot集成Apollo 能够实现Apollo适用功能[自动刷新.监听器.配置加密.1灰度发 ...
- 【Github教程】史上最全github用法:github入门到精通
原文 http://www.eoeandroid.com/thread-274556-1-1.html [初识Github] 首先让我们大家一起喊一句"Hello Github". ...
- 史上最全的MonkeyRunner自动化测试从入门到精通(10)
三.MonkeyRunner复杂的功能开始学习 (1)获取APK文件中ID的两种方式 Monkeyrunner的环境已经搭建完成,现在对Monkeyrunner做一个简介. Monkeyrunner工 ...
- python编程从入门到精通pdf-码哥—-Python从入门到精通(全60集),资源教程下载...
课程名称 码哥--python从入门到精通(全60集),资源教程下载 课程目录 1Python编程语言历史及特性 2Python编程语言初接触 3Python程序文件结构 4准备Python编程环境 ...
- python从入门到精通视频(全60集)-【网盘下载】Python从入门到精通视频(全60集)...
Python从入门到精通视频(全60集) JAVA? 今天 课程目录 python入门教程-1-Python编程语言历史及特性 python入门教程-2-Python编程语言初接触 python入门教 ...
最新文章
- 创建细分客户的无监督学习项目
- 拼多多程序员:因为看到同事被抬上救护车,我被开除了!
- hdu1074 状态压缩dp+记录方案
- 金庸小说人物关系表年表
- 使用StackBlitz和SAP Spartacus快速创建电商店铺页面
- PKU 学生反馈 3
- 运维前线:一线运维专家的运维方法、技巧与实践3.1 数据中心搬迁准备
- Ural_1586. Threeprime Numbers (DP)
- wxpython菜单的位置_wx.grid 怎么加菜单呀!
- Membership学习(三)Membership Providers介绍
- 【渝粤教育】国家开放大学2018年春季 0266-22T设计构成 参考试题
- [2018.10.17 T1] 斜率
- 操作系统ucore lab1
- JS+html+css实现的一个小小的贪吃蛇游戏
- 学计算机跨考航天航空,往届生跨考计算机经验
- 如何做好一次视频号直播
- 暗影骑士vn7安装linux,宏碁暗影骑士vn7 591g笔记本怎么样?
- RTMPLive多流媒体协议转无插件直播协议(RTMP)解决方案
- 【转】常用的隧道技术
- 无论LCD和OLED显示技术有多好,都无法替代这个古老的显示数码管
热门文章
- 8个高质量免抠素材网站
- (一)目标检测经典模型回顾
- 深度学习声纹识别_一种基于机器学习及深度学习的声纹降噪方法及系统与流程...
- Microsoft Word 教程:如何在 Word 中插入图片、图标?
- boost电路输出电流公式_BOOST电路参数计算公式
- 当今主流软件产品家族及其开发语言
- 广告数据定量分析:第一章——广告优化中的统计学
- 神州数码交换机配置基本命令
- 记一次https的免费ssl证书安装及配置过程--window系统
- android 时钟动态图标,神奇的Android动态时钟/时间UI界面