读取Unique reads

从贴好的bam文件中读取unique mapped reads，来进行下一步的分析。
tophat贴的结果对于非unique的reads会反复出现在bam/sam中，因此不挑unique reads后续结果会不准确。
Tophat结果中Unique reads的Quality一栏为255，而BWA贴出的结果此栏会根据reads本身的Phred值变化，因此确定unique的唯一标准是optional fields中的NH:i:1。
直接
samtools view accepted.bam | grep "\bNH:i:1\b" > accepted_unique.sam
即可。其中”\bNH:i:1\b"如果忘记\b，会挑出NH:i:12之类的非unique reads。

读取Unique reads相关推荐

RNA 1. SCI 文章中读取 GEO 数据
依稀记得10年前一个样本的RNA-SEQ费用还是蛮高的,而现在就是洒洒水啦,所以这种转录组的数据也已经成为文章的主流,基本上就是WES结合RNA,但是个人感觉,真的能关联到突变基因和表达的水平上去的, ...
读取STM32单片机-ID操作
每一片 STM32 芯片内部拥有一个独一无二的 Unique Device ID, 96 Bit. 这个 ID 号可以提供给开发者很多优越的功能, 例如: 1. 可以把 ID 做为用户最终产品的序列号 ...
Microbiome：芝麻菜中肠杆菌科主导核心微生物组并贡献抗生素抗性组(简单套路16S+meta+培养组发高分文章)
文章目录日报文章思路总结摘要主要结果图1. 三类样本的细菌组成图2. 宏基因组中肠杆菌群体结构和丰度图3. 叶际和根际中肠杆菌科的核心微生物组图4. 芝麻菜抗性组评估图5. 可食用植 ...
驱动专题：第四章MTD及Flash驱动 3. nandflash 详细硬件特性
[详解]如何编写Linux下Nand Flash驱动版本:v2.2.1 Crifan Li 摘要本文先解释了Nand Flash相关的一些名词,再从Flash硬件机制开始,介绍到Nand Flas ...
SQL Server 2005 术语词汇表
术语定义 ActiveX 数据对象 (ActiveX Data Objects) 一种易于使用的应用程序编程接口 (API),用于封装 OLE DB 以在诸如 Visual Basic.Visual ...
Mothur1_Mothur的简介及运行准备
本文主要介绍生物信息学软件 Mothur及其分析准备,接下来我们将持续更新Mothur中文教程. Mothur是由密歇根大学微生物学和免疫学系的Patrick Schloss教授和其研究小组开发的一款 ...
nand flash 个人觉得写得比较好的文章
[详解]如何编写Linux下Nand Flash驱动版本:v2.2.1 Crifan Li 摘要本文先解释了Nand Flash相关的一些名词,再从Flash硬件机制开始,介绍到Nand Flas ...
nand flash 经典全面 ------如何编写Linux下Nand Flash驱动
Crifan Li 摘要本文先解释了Nand Flash相关的一些名词,再从Flash硬件机制开始,介绍到Nand Flash的常见的物理特性,且深入介绍了Nand Flash的一些高级功能,然后开 ...
2020转录组RNA-SEQ上游分析
文章目录安装配置conda 使用清华源下载sh脚本并安装设置镜像源 conda环境创建激活和退出环境 conda安装软件质量评估 @ fastQC fastq格式 FsatQC软件 multi ...
NGS数据分析实践：05. 测序数据的基本质控 [2] - MultiQC
NGS数据分析实践:05. 测序数据的基本质控 [2] - MultiQC 2. MultiQC 2.1 帮助信息及运行代码 2.2 报告解读 2.3 小结文接上篇:NGS数据分析实践:05. 测序 ...

读取Unique reads

读取Unique reads相关推荐

最新文章

热门文章