html文件下的flag,推荐一个SAM文件中flag含义解释工具

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

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

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

bee4ab4cbee880da711a9dd9be99f385.png

脚本工具的使用方法:

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

Explain SAM Flags

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

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

["read unmapped", 0x4],

["mate unmapped", 0x8],

["read reverse strand", 0x10],

["mate reverse strand", 0x20],

["first in pair", 0x40],

["second in pair", 0x80],

["not primary alignment", 0x100],

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

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

function explainFlags() {

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

var summary = "";

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

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

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

summary += "     " + lstFlags[i][0] + "
";

checkbox.checked = true;

} else {

checkbox.checked = false;

}

}

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

}

function checkboxClicked() {

//compute the new flag value

var newFlagValue = 0;

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

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

if(checkBox.checked) {

newFlagValue |= lstFlags[i][1];

}

}

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

textbox.value = newFlagValue;

explainFlags();

}

This page requires JavaScript. Please enable it in your browser settings.

This utility explains SAM flags in plain English.

Flag:  

Explanation:

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

document.write("   " +lstFlags[i][0] + "
");

}

  read paired   read mapped in proper pair   read unmapped   mate unmapped   read reverse strand   mate reverse strand   first in pair   second in pair   not primary alignment   read fails platform/vendor quality checks   read is PCR or optical duplicate

Summary:

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值