關(guān)于blast結(jié)果的一點解讀

alignments 代表比對上的兩個序列

hits 表示兩個序列比對上的片段

Score 比對得分绞呈,如果序列匹配上得分所计,不一樣,減分渴邦,分值越高,兩個序列相似性越高

E Value 值越小拘哨,越可信谋梭,相對的一個統(tǒng)計值。這與你所使用的數(shù)據(jù)庫大小有關(guān)

Length 輸入序列的長度

Identities 一致性倦青,就是兩個序列有多少是一樣的

Query 代表輸入序列

Sbjct 代表數(shù)據(jù)庫中的序列

blast是區(qū)段比對瓮床,對于給定的兩個序列,blast會把具有相識性的片段(hit)找出來产镐,顯示的是hit的信息

E值的經(jīng)驗解釋如下隘庄。如果e<1e-50(或1×10-50),那么數(shù)據(jù)庫匹配應(yīng)該是同源關(guān)系的結(jié)果癣亚,這是一個非常高的置信度丑掺。如果e介于0.01和1e-50之間,則可以將匹配視為同源性的結(jié)果述雾。如果e介于0.01和10之間街州,則認(rèn)為匹配不重要,但可能暗示存在暫時的遠(yuǎn)程同源關(guān)系玻孟。需要額外的證據(jù)來確認(rèn)暫時的關(guān)系唆缴。如果e>10,考慮中的序列要么不相關(guān)黍翎,要么通過極遠(yuǎn)的關(guān)系相關(guān)面徽,這些關(guān)系低于當(dāng)前方法的檢測極限。

由于e值受數(shù)據(jù)庫大小的比例影響匣掸,一個明顯的問題是斗忌,隨著數(shù)據(jù)庫的增長质礼,給定序列匹配的e值也會增加。由于兩個序列之間真正的進(jìn)化關(guān)系保持不變织阳,隨著數(shù)據(jù)庫的增長眶蕉,序列匹配的可信度降低意味著隨著數(shù)據(jù)庫的擴大,人們可能會“失去”先前檢測到的同系物唧躲。因此造挽,需要一種替代E值計算的方法。

蛋白質(zhì)由20種不同的aas組成弄痹,如果你將兩個不相關(guān)的蛋白質(zhì)序列(或任何其他隨機a a序列)與任何長度對齊饭入,你將擁有5%的隨機特性(對于DNA和RNA序列,隨機特性為25%肛真,因為它們是由不同堿基a谐丢、t、c蚓让、g組成的)乾忱。

Identities?35%意味著序列中35%的aa與數(shù)據(jù)庫中的其他序列匹配。取決于你在尋找什么:——如果你有未知的蛋白質(zhì)序列历极,你想知道同源序列窄瘟,關(guān)于身份的信息(甚至35%)是有價值的,——如果你知道蛋白質(zhì)趟卸,你需要確認(rèn)序列蹄葱,身份35%是小的,可能表明在你的分析過程中出了問題锄列。

BLASTP中還有一個基于相似性的參數(shù)ppos图云。ppos是pident+(相似但不相同的aa匹配的百分比)。兩個pident高于20%和ppos高于30%的aa序列非常接近邻邮,可以稱為同源序列竣况。在NA序列中,PIDENT 40%及以上是可以的饶囚。

p值:取決于查詢和db長度帕翻,但我認(rèn)為小于10^-5的p值表示一個關(guān)系鸠补。

bitscore:很大程度上取決于查詢長度萝风。將bitscore與您的qlen進(jìn)行比較,我認(rèn)為如果一個命中的bitscore等于或大于qlen的0.7紫岩,那么查詢和主題就足夠接近了规惰。


https://www.biostars.org/p/187230/

http://boyun.sh.cn/bio/?p=1453

http://boyun.sh.cn/bio/?p=1632

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市泉蝌,隨后出現(xiàn)的幾起案子歇万,更是在濱河造成了極大的恐慌揩晴,老刑警劉巖,帶你破解...
    沈念sama閱讀 211,817評論 6 492
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件贪磺,死亡現(xiàn)場離奇詭異硫兰,居然都是意外死亡,警方通過查閱死者的電腦和手機寒锚,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 90,329評論 3 385
  • 文/潘曉璐 我一進(jìn)店門劫映,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人刹前,你說我怎么就攤上這事泳赋。” “怎么了喇喉?”我有些...
    開封第一講書人閱讀 157,354評論 0 348
  • 文/不壞的土叔 我叫張陵祖今,是天一觀的道長。 經(jīng)常有香客問我拣技,道長千诬,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 56,498評論 1 284
  • 正文 為了忘掉前任过咬,我火速辦了婚禮大渤,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘掸绞。我一直安慰自己泵三,他們只是感情好,可當(dāng)我...
    茶點故事閱讀 65,600評論 6 386
  • 文/花漫 我一把揭開白布衔掸。 她就那樣靜靜地躺著烫幕,像睡著了一般。 火紅的嫁衣襯著肌膚如雪敞映。 梳的紋絲不亂的頭發(fā)上较曼,一...
    開封第一講書人閱讀 49,829評論 1 290
  • 那天,我揣著相機與錄音振愿,去河邊找鬼捷犹。 笑死,一個胖子當(dāng)著我的面吹牛冕末,可吹牛的內(nèi)容都是我干的萍歉。 我是一名探鬼主播,決...
    沈念sama閱讀 38,979評論 3 408
  • 文/蒼蘭香墨 我猛地睜開眼档桃,長吁一口氣:“原來是場噩夢啊……” “哼枪孩!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 37,722評論 0 266
  • 序言:老撾萬榮一對情侶失蹤蔑舞,失蹤者是張志新(化名)和其女友劉穎拒担,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體攻询,經(jīng)...
    沈念sama閱讀 44,189評論 1 303
  • 正文 獨居荒郊野嶺守林人離奇死亡从撼,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,519評論 2 327
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了钧栖。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片谋逻。...
    茶點故事閱讀 38,654評論 1 340
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖桐经,靈堂內(nèi)的尸體忽然破棺而出毁兆,到底是詐尸還是另有隱情,我是刑警寧澤阴挣,帶...
    沈念sama閱讀 34,329評論 4 330
  • 正文 年R本政府宣布气堕,位于F島的核電站,受9級特大地震影響畔咧,放射性物質(zhì)發(fā)生泄漏茎芭。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 39,940評論 3 313
  • 文/蒙蒙 一誓沸、第九天 我趴在偏房一處隱蔽的房頂上張望梅桩。 院中可真熱鬧,春花似錦拜隧、人聲如沸宿百。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,762評論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽垦页。三九已至,卻和暖如春干奢,著一層夾襖步出監(jiān)牢的瞬間痊焊,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,993評論 1 266
  • 我被黑心中介騙來泰國打工忿峻, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留薄啥,地道東北人。 一個月前我還...
    沈念sama閱讀 46,382評論 2 360
  • 正文 我出身青樓逛尚,卻偏偏與公主長得像垄惧,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子黑低,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 43,543評論 2 349

推薦閱讀更多精彩內(nèi)容

  • Blast赘艳,全稱Basic Local Alignment Search Tool,即"基于局部比對算法的搜索工具...
    曉僉閱讀 14,896評論 1 26
  • ORA-00001: 違反唯一約束條件 (.) 錯誤說明:當(dāng)在唯一索引所對應(yīng)的列上鍵入重復(fù)值時克握,會觸發(fā)此異常蕾管。 O...
    我想起個好名字閱讀 5,256評論 0 9
  • 歡迎關(guān)注:oddxix 本章主要講序列分析與聯(lián)配 序列分析是生物信息學(xué)最主要的研究內(nèi)容之一,它可以分為兩個主要部分...
    oddxix閱讀 2,237評論 0 10
  • 國家電網(wǎng)公司企業(yè)標(biāo)準(zhǔn)(Q/GDW)- 面向?qū)ο蟮挠秒娦畔?shù)據(jù)交換協(xié)議 - 報批稿:20170802 前言: 排版 ...
    庭說閱讀 10,930評論 6 13
  • 官網(wǎng) 中文版本 好的網(wǎng)站 Content-type: text/htmlBASH Section: User ...
    不排版閱讀 4,370評論 0 5