sam 文件的解释

高通两数据测序之后要map到基因组上，我通常用bwa bowtie和tophat 。其实我知道sam 文件有好几个field，每个field代表不同的意思，但是由于各种格式的文件确实太多了，我也记不住每个field具体的意思。每次要用的时候都是上网再去找，或者把相关的资料下载下来，再现学。

今天，我和一个同学讨论我们数据的情况。他们都是用bowtie 去map 的，我之前写了个脚本，是用bwa 做的map，我有点懒散，也懒得再去修改了，于是也没想着要再写一个用bowtie 做map的脚本。与同学的讨论中，我发现我之前对sam 文件的了解确实很肤浅，这是一个教训，以后每一步都要走踏实才好。

　　sam 文件中第二列的信息是FLAG 标签

　　0x0001　　　　reads is paired in sequence, no matter whether it is mapped in pair

　　0x0002　　　　the reads is mapped in a proper pair

　　0x0004　　　　the query sequence itself is unmapped

　　0x0008　　　　the mate sequence is unmapped

　　0x0010　　　　strand of the query is reverse complement

　　　　……

　　　　……

　　　　……

　　可以通过按位运算操作符筛选出特定的reads。

　　　　

　　0x0001　　　　　　1

　　0x0002　　　　　  10

　　0x0004　　　　　100

　　0x0008　　　　 1000

　　……

　　……

sam 文件的解释,布布扣,bubuko.com

时间： 2024-12-26 03:31:12

sam 文件的解释的相关文章

31、SAM文件中flag含义解释工具--转载

转载:http://www.cnblogs.com/nkwy2012/p/6362996.html SAM是Sequence Alignment/Map 的缩写.像bwa等软件序列比对结果都会输出这样的文件.samtools网站上有专门的文档介绍SAM文件.具体地址:http://samtools.sourceforge.net/SAM1.pdf 很多人困惑SAM文件中的第二列FLAG值是什么意思.根据文档介绍我们可以计算,但是为了方便大家,下面给大家提供一个脚本工具,大家直接输入flag值就可

linux 里 /etc/passwd 、/etc/shadow和/etc/group 文件内容解释

•/etc/passwd文件用于存放用户账户信息,每行代表一个账户,每个账户的各项信息用冒号分割,例如: root:x:0:0:root:/root:/bin/bash username:password:uid:gid:allname:homedir:shell 1.用户名. 2.密码,“x”代表密码存储在/etc/shadow中:当该值为其他任意非“x”值时,可以通过root用户切换(不需要密码),但是都无法通过非root用户切换到相应用户,因为无法获得起密码了. 3.uid,用户id,0-

T-SQL 操作文件 具体解释

/*******  导出到excel EXEC master..xp_cmdshell 'bcp SettleDB.dbo.shanghu out c:\temp1.xls -c -q -S"GNETDATA/GNETDATA" -U"sa" -P""' /***********  导入Excel SELECT * FROM OpenDataSource( 'Microsoft.Jet.OLEDB.4.0', 'Data Source="

ArcGIS学习记录—dbf shp shx sbn sbx mdb adf等类型的文件的解释

原文地址: ArcGIS问题:dbf shp shx sbn sbx mdb adf等类型的文件的解释 - Silent Dawn的日志 - 网易博客 http://gisman.blog.163.com/blog/static/34493388201022254341339/ 这些文件英文解释如下(看不懂可以看下面的中文解释): Shapefile file extensions shp-The main file that stores the feature geometry. Requi

数据库config.php文件内容解释及替换方法

在网上我们下载的程序,有很多并没有安装程序,而是需要我们手动修改数据库配置文件 有的还需要网站路径如域名等,下面由魅力网络为大家解释 <?php$pe['db_host'] = 'localhost';服务器名$pe['db_name'] = 'cps1';数据库名$pe['db_user'] = 'root';数据库用户名$pe['db_pw'] = '';数据库密码$pe['db_coding'] = 'utf8';编码$pe['url_model'] = 'pathinfo_safe';安

Android清单文件具体解释(二) ---- 应用程序权限声明

我们知道,Android系统的各个模块提供了很强大的功能(比方电话,电源和设置等),通过使用这些功能.应用程序能够表现的更强大.更灵活.只是,使用这些功能并非无条件的.而是须要拥有一些权限.接下来,我们就開始解说还有一个很重要的知识点--应用程序权限声明.当中主要包含应用程序的权限声明,自己定义应用程序的訪问权限和SDK版本号限定. 1.<uses-permission>--应用程序的权限申请 权限 描写叙述 android.permission.ACCESS_NETWORK_STATE 同意

.ssh/config 文件的解释算法及配置原则

前言 SSH 是连接远程主机最常用的方式,尽管连接到耽搁主机的基本操作非常直接,但当你开始使用大量的远程系统时,这就会成为笨重和复杂的任务. 幸运的是,OpenSSH 允许您提供自定义的客户端连接选项.这些选项可以被存储到一个配置文件中,这个配置文件可以用来定义每个主机的配置.这有助于保持每个主机的连接选项更好的独立和组织,也你让你在需要连接时避免在命令行中写繁琐的选项. 在本文中,我们将介绍 SSH 客户端配置文件的基本知识,了解 SSH 如何解释配置文件,但不会介绍 SSH 具体选项的含义及

Android清单文件具体解释(四) ---- backupAgent的使用方法

在<application>节点中有一个很重要的属性,那就是backupAgent.这里我们将它单独列出来,从基本含义,使用方法及其相关属性等方面来具体介绍一下. 1.backupAgent简单介绍 android:backupAgent用来设置备份代理.对于大部分应用程序来说,都或多或少保存着一些持久性的数据,比方数据库和共享文件.或者有自己的配置信息.为了保证这些数据和配置信息的安全性以及完整性.Android提供了这样一个机制. 我们能够通过这个备份机制来保存配置信息和数据以便为应用程序

git学习-git status, git diff, git rm, git mv等文件操作命令解释。

直接说每条命令了 git status是显示当前文件所处于的状态. ? /Users/alps/Sites/judianer/1 git:(master)>git status On branch master Your branch is up-to-date with 'origin/master'. nothing to commit, working directory clean 这个是干净的情况,就是全部push都提交了,没有改动过的文件. ? /Users/alps/Sites/j