SAM是Sequence Alignment/Map 的縮寫(xiě)。像bwa等軟件序列比對(duì)結(jié)果都會(huì)輸出這樣的文件。samtools網(wǎng)站上有專門(mén)的文檔介紹SAM文件。具體地址：http://samtools.sourceforge.net/SAM1.pdf

很多人困惑SAM文件中的第二列FLAG值是什么意思。根據(jù)文檔介紹我們可以計(jì)算，但是為了方便大家，下面給大家提供一個(gè)腳本工具，大家直接輸入flag值就可以知道它代表的含義了。

該腳本的使用方法如下截圖所示：

腳本工具的使用方法：

將下面的代碼保存在記事本里面，另存為一個(gè)html文件，如文件名：FlagExplain.html(拓展名一定要為.html)。雙擊既可以在瀏覽器里面打開(kāi)了。

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:

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來(lái)咯，堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)

總結(jié)

以上是生活随笔為你收集整理的html文件下的flag,推荐一个SAM文件中flag含义解释工具的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。