计算基因上外显子碱基覆盖度（exon coverage depth）：Samtool工具使用

假设想要计算ATP1A4基因上的外显子碱基覆盖度

首先查询这个基因所有exon的起始和终止位置，查询链接：http://grch37.ensembl.org/Homo_sapiens/Transcript/Exons?db=core;g=ENSG00000132681;r=1:160121360-160156767;t=ENST00000368081

如下图所示，可以看到该基因有6个转录本，假定选择第二个转录本ENST00000470705.1，点击进去，显示该转录本有5个外显子

具体外显子的详细信息如下所示：

获得该基因的所有外显子的起始终止位置后，就可以用Samtool提取该位置，并且计算碱基覆盖度了，用到的命令如下：

1、先提取外显子的起始终止位置：

/software/samtools view -h file.bam 1:160147201-16014454 > ATP1A4_exon1.sam #提取基因ATP1A4上转录本ENST00000470705.1的第一个外显子片段

2、计算每个exon的碱基覆盖度：

/software/samtools depth ATP1A4_exon1.sam > ATP1A4_exon1_coverage.txt #计算基因ATP1A4上转录本ENST00000470705.1的第一个外显子碱基覆盖度

转载于:https://www.cnblogs.com/chenwenyan/p/8417973.html

计算基因上外显子碱基覆盖度（exon coverage depth）：Samtool工具使用相关推荐

【perl计算基因在基因组上的距离】
perl计算基因在基因组上的距离介绍根据输入注释文件及目的基因文件,找出目的基因在基因组上的距离并输出到文件内一.文件格式基因注释文件:为细菌全基因组利用prokka注释文件,格式如下 gen ...
html计算平均分,Calculate phastCon Score for a gene —- 计算基因的phastCon平均分，判断基因保守型...
Calculate phastCon Score for a gene -- 计算基因的phastCon平均分,判断基因保守型 PhastCon socre is the score from 0 t ...
文献速递 | 新研究揭示一系列能够提高Cas9基因编辑和碱基编辑特异性的选择性核输出抑制剂
CRISPR基因编辑 CRISPR基因编辑是一项"年轻"的技术,但在疾病治疗中已经展示出广泛的应用潜力.在临床方案中,CRISPR-Cas9往往为组成型表达,容易导致Cas9核酸酶 ...
Michael Jordan在“21世纪的计算大会”上的主题演讲：大数据的分治和统计推断
演讲人:Michael Jordan博士美国加州大学伯克利分校电气工程与计算机系.统计学系教授美国科学院(NAS).美国工程院(NAE).美国文理科学院(AAAS)院士演讲标题:Divide-a ...
AWS计算云上下载代码命令 git clone url
AWS计算云上下载代码命令 git clone url https://blog.csdn.net/qq_42672770/article/details/81317778
php 计算两点时间距离,PHP计算地球上两点之间的距离（示例详解）
给定经度和纬度,求地球上两点之间的距离.首先我们需要了解该问题的解决思路,然后再用PHP代码来实现计算. 此问题可以用半正矢(haversine)公式求解: 大圆距离或正交距离是球面(或地球表面)上两 ...
谷歌量子计算登上Nature封面，首次实现量子优越性，里程碑式突破
这是量子计算领域的里程碑事件:谷歌在被称为「量子优越性」方向上的重大突破研究,登上了<自然>杂志 150 周年版的封面. 它或许标志着量子计算正在走向实用化:谷歌已经利用一台 54 量子比 ...
ITK:计算网格上的测地距离
ITK:计算网格上的测地距离内容提要输出结果输入输出 C++实现代码内容提要从网格上提供的seed 顶点计算测地距离. 输出结果输入输出 C++实现代码 #include " ...
php 地图两点距离计算,计算地图上两点间的距离PHP类
计算地图上两点间的距离,使用的是谷歌地图 class GeoHelper { /** * @param int $lat1 * @param int $lon1 * @param int $lat2 ...
数据库：计算地球上两个坐标点之间里程
今天给大家分享计算地球上两个坐标点之间里程不同数据库版本的脚本. 1.SQLServer脚本 –-计算地球上两个坐标点(经度,纬度)之间距离sql函数 CREATE FUNCTION [dbo].[f ...

计算基因上外显子碱基覆盖度（exon coverage depth）：Samtool工具使用

计算基因上外显子碱基覆盖度（exon coverage depth）：Samtool工具使用相关推荐

最新文章

热门文章