推荐一个SAM文件中flag含义解释工具--转载

SAM是Sequence Alignment/Map 的缩写。像bwa等软件序列比对结果都会输出这样的文件。samtools网站上有专门的文档介绍SAM文件。具体地址：http://samtools.sourceforge.net/SAM1.pdf

很多人困惑SAM文件中的第二列FLAG值是什么意思。根据文档介绍我们可以计算，但是为了方便大家，下面给大家提供一个脚本工具，大家直接输入flag值就可以知道它代表的含义了。

该脚本的使用方法如下截图所示：

脚本工具的使用方法：

将下面的代码保存在记事本里面，另存为一个html文件，如文件名：FlagExplain.html（拓展名一定要为.html）。双击既可以在浏览器里面打开了。

查看源代码

打印帮助

`01`	`<html><head><meta` `http-equiv="Content-Type"` `content="text/html; charset=GBK">`

`02`	`<title>Explain SAM Flags</title>`

`03`	`<script` `type="text/javascript">`

`04`	`lstFlags = [["read paired", 0x1],`

`05`	`["read mapped in proper pair", 0x2],`

`06`	`["read unmapped", 0x4],`

`07`	`["mate unmapped", 0x8],`

`08`	`["read reverse strand", 0x10],`

`09`	`["mate reverse strand", 0x20],`

`10`	`["first in pair", 0x40],`

`11`	`["second in pair", 0x80],`

`12`	`["not primary alignment", 0x100],`

`13`	`["read fails platform/vendor quality checks", 0x200],`

`14`	`["read is PCR or optical duplicate", 0x400]];`

15

`16`	`function explainFlags() {`

`17`	`var flagValue = parseInt(document.getElementById('tb').value); //returns 0 or NaN if can't parse`

`18`	`var summary = "";`

`19`	`for(var i = 0; i <` `lstFlags.length; i++) {`

`20`	`var` `checkbox` `=` `document.getElementById('cb' + i)`

`21`	`if(lstFlags[i][1] & flagValue) {`

`22`	`summary += "     " + lstFlags[i][0] + "<br>";`

`23`	`checkbox.checked = true;`

`24`	`} else {`

`25`	`checkbox.checked = false;`

26 }

27 }

28

`29`	`document.getElementById('summary').innerHTML = summary;`

30 }

31

`32`	`function checkboxClicked() {`

`33`	`//compute the new flag value`

`34`	`var newFlagValue = 0;`

`35`	`for(var i = 0; i <` `lstFlags.length; i++) {`

`36`	`var` `checkBox` `=` `document.getElementById('cb' + i);`

`37`	`if(checkBox.checked) {`

`38`	`newFlagValue \|= lstFlags[i][1];`

39 }

40 }

`41`	`var` `textbox` `=` `document.getElementById('tb');`

`42`	`textbox.value` `=` `newFlagValue;`

`43`	`explainFlags();`

44 }

`45`	`</script>`

46

`47`	`<noscript>This page requires JavaScript. Please enable it in your browser settings.</noscript>`

48 </head>

49 <body>

50

`51`	`This utility explains SAM flags in plain English. <br>`

52 <br>

53

`54`	`<form` `onsubmit="explainFlags(); return false;">`

55 Flag: 

`56`	`<input` `id="tb"` `type="text"` `size="10">      `

`57`	`<input` `type="submit"` `value="Explain"><br>`

58 <br>

`59`	`Explanation:<br>`

`60`	`<script` `type="text/javascript">`

`61`	`for(var i = 0; i <` `lstFlags.length; i++) {`

`62`	`document.write("<input` `type=checkbox` `name=cb" + i + "` `id='cb" + i + "'onclick='checkboxClicked();'>   " +lstFlags[i][0] + "</input><br>");`

63 }

64 </script><input type="checkbox" name="cb0" id="cb0"onclick="checkboxClicked();"> read paired<br><input type="checkbox"name="cb1" id="cb1" onclick="checkboxClicked();"> read mapped in proper pair<br><input type="checkbox" name="cb2" id="cb2"onclick="checkboxClicked();"> read unmapped<br><input type="checkbox"name="cb3" id="cb3" onclick="checkboxClicked();"> mate unmapped<br><input type="checkbox" name="cb4" id="cb4" onclick="checkboxClicked();"> read reverse strand<br><input type="checkbox" name="cb5" id="cb5"onclick="checkboxClicked();"> mate reverse strand<br><inputtype="checkbox" name="cb6" id="cb6" onclick="checkboxClicked();"> first in pair<br><input type="checkbox" name="cb7" id="cb7"onclick="checkboxClicked();"> second in pair<br><input type="checkbox"name="cb8" id="cb8" onclick="checkboxClicked();"> not primary alignment<br><input type="checkbox" name="cb9" id="cb9"onclick="checkboxClicked();"> read fails platform/vendor quality checks<br><input type="checkbox" name="cb10" id="cb10"onclick="checkboxClicked();"> read is PCR or optical duplicate<br>

65 <br>

`66`	`Summary:<br>`

`67`	`<div` `id="summary">`

`68`	`</div></form></body></html>`

随机推荐

实习培训——Java基础（4）
实习培训——Java基础(4) 1 多态多态是同一个行为具有多个不同表现形式或形态的能力. 多态就是同一个接口,使用不同的实例而执行不同操作,多态性是对象多种表现形式的体现. 现实中,比如我们按下 ...
四、Mosquitto 高级应用之用户配置
本文将讲解 Mosquitto 用户机制.如果还没有搭建 Mosquitto 服务的可以参考我的另外两篇文章<< 一.Mosquitto 介绍&安装>> << ...
python3中pymysql模块安装及连接数据库（同‘python安装HTMLTestRunner’）
https://pypi.org/project/PyMySQL/#files 安装完成之后就是连接数据库了机器上安装了mysql的数据库,并且已经创建了两张表用于测试 python3连接数据库及删 ...
通过脚本获取form表单的数值而不是submit
jQuery的serialize()方法通过序列化表单值,创建URL编码文本字符串,我们就可以选择一个或多个表单元素,也可以直接选择form将其序列化,如: <form action=" ...
dedecms开启报错
php.ini里面设置display_errors = On 开启错误提示,error_reporting = E_ALL & ~E_NOTICE 设置错误等级.也可以在php文件中ini_s ...
VS2010/MFC编程入门之三十一（常用控件：树形控件Tree Control 下）
前面一节讲了树形控件Tree Control的简介.通知消息以及相关数据结构,本节继续讲下半部分,包括树形控件的创建.CTreeCtrl类的主要成员函数和应用实例. 树形控件的创建 MFC为树形控件提 ...
小试---EF5.0简介
简介实体框架Entity Framework 是 ADO.NET 中的一组支持开发面向数据的软件应用程序的技术.是微软的一个ORM框架.简单的说就是把关系型数据库映射成面向对象模型. 一篇更加详细的 ...
复制控件出错--提示XX控件不存在
有时候比较喜欢复制粘贴,但是结果就是,经常出现XX控件找不到,或者不存在之类的话~~ 错误如下: 在相当郁闷的时候,尝试了一种很土土的方法,就是直接新建一个同名新页面,把前台和后台的代码复制进去~~ ...
彻底明白Flink系统学习5：window、Linux本地安装Flink
http://www.aboutyun.com/thread-26393-1-1.html 问题导读 1.如何在window下安装Flink? 2.Flink本地安装启动命令与原先版本有什么区别? 3 ...
P2414 [NOI2011]阿狸的打字机
P2414 [NOI2011]阿狸的打字机 AC自动机+树状数组优质题解 <------题目分析先AC自动机搞出Trie图然后根据fail指针建一只新树把树映射(拍扁)到一个序列上,用树 ...

推荐一个SAM文件中flag含义解释工具--转载

推荐一个SAM文件中flag含义解释工具--转载的更多相关文章

随机推荐

热门专题