| 查看: 2116 | 回復(fù): 12 | ||
bettyshan木蟲 (正式寫手)
夢(mèng)幻biobrick
|
[求助]
批量序列比對(duì)分析 已有2人參與
|
| 想同時(shí)比對(duì)幾百個(gè)基因在某一個(gè)基因組中是否存在同源基因,如何實(shí)現(xiàn)? |

木蟲 (小有名氣)
木蟲 (正式寫手)
夢(mèng)幻biobrick

鐵桿木蟲 (正式寫手)
木蟲 (正式寫手)
夢(mèng)幻biobrick

木蟲 (正式寫手)
夢(mèng)幻biobrick

木蟲 (正式寫手)
夢(mèng)幻biobrick

鐵桿木蟲 (正式寫手)
|
用本地blast是要下載所有鏈霉菌的蛋白序列的,NCBI里好像沒有單獨(dú)將鏈霉菌的蛋白整理出來,你可以下載所有細(xì)菌的蛋白,ftp://ftp.ncbi.nlm.nih.gov/refseq/release/bacteria/這個(gè)FTP里頭,進(jìn)去之后網(wǎng)頁(yè)拉到下面一點(diǎn),下載那些帶“nonredundant” 和“faa”的,所有231個(gè)文件,蛋白估計(jì)也就10多個(gè)G,總能下載完的。下載完之后可以從這些序列中挑出鏈霉菌的序列,然后構(gòu)建blast數(shù)據(jù)庫(kù),這一步估計(jì)要寫一些計(jì)算機(jī)代碼,不知道你會(huì)不會(huì)。 上面的方法最準(zhǔn)確但是有些麻煩,你也可以到uniprot網(wǎng)站http://www.uniprot.org/里頭,搜索streptomyces,出來的就是鏈霉菌的蛋白,大約是60萬(wàn)條,搜索完后這個(gè)網(wǎng)站會(huì)提供序列下載,點(diǎn)擊download就行了,下載fasta格式的序列可以直接用來構(gòu)建blast數(shù)據(jù)庫(kù)。 有blast數(shù)據(jù)庫(kù)之后比對(duì)就比較簡(jiǎn)單了,難的估計(jì)是分析比對(duì)結(jié)果。你說的e-value比較好辦, identity和 matching我不知道該怎么控制,不過用blast的-m8參數(shù)可以輸出Tab類型的結(jié)果,要篩選 identity和 matching應(yīng)該也不難。 我以前做過類似你說的東西,不過我找的是真菌蛋白,當(dāng)時(shí)找到的真菌序列大約是200萬(wàn),我估計(jì)鏈霉菌的也是百萬(wàn)級(jí)別。百萬(wàn)級(jí)別的序列不算多,想要完成你說的這類任務(wù),你得習(xí)慣同這種大數(shù)據(jù)打交道。 |
木蟲 (正式寫手)
夢(mèng)幻biobrick

木蟲 (正式寫手)
夢(mèng)幻biobrick

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +3 | gdlf9999 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠(chéng)心求調(diào)劑 +12 | 吃吃吃才有意義 2026-03-19 | 13/650 |
|
|
[考研] 085600材料與化工調(diào)劑 +7 | A-哆啦Z夢(mèng) 2026-03-23 | 12/600 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +5 | niko- 2026-03-22 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +3 | WW.' 2026-03-23 | 5/250 |
|
|
[考研] 316求調(diào)劑 +7 | 梁茜雯 2026-03-19 | 7/350 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 9/450 |
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|