 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 2533 | 回復: 7 | |||
cyanobacter金蟲 (小有名氣)
|
[求助]
高通量數(shù)據(jù)做本地blast時，如何將97%以上的比對結(jié)果全部顯示出來？ 已有3人參與
|
|
各位蟲友，我有一批16S的高通量數(shù)據(jù)，然后我現(xiàn)在想做一個本地blast，把我手上的一部分序列（isolate.fasta）與高通量數(shù)據(jù)進行blast，然后把所有高于97%的結(jié)果全部顯示出來，請大家?guī)兔Α?br />
我是這樣做的，首先用格式化本地數(shù)據(jù)庫（即16S高通量數(shù)據(jù)515F-907R.format.fq.fasta），用的是formatdb命令 formatdb -i 515F-907R.format.fq.fasta -p F -t My database -o T 然后生成三個擴展名分別是：.pin、.psq、.phr的文件，然后就開始進行本地blast，用blastall命令 blastall -p blastn -i isolate.fasta -d 515F-907R.format.fq.fasta -o myresult.blastout -T T -e 1e-10 -m 9 這樣運行出來會有一個所有isolate.fasta的序列與庫序列比對的結(jié)果，但是每個比對結(jié)果只顯示前250個，有些可能比到第250個相似性也是99%以上，就說說還有許多庫序列與query的相似性高于97%，但是這部分就不顯示出來了。 請問一下大家，可以用什么程序解決這個問題，可不可以設(shè)置一個97%的閾值，比對高于這個值得全部顯示出來，低于97%的不顯示。剛剛接觸生物信息方面的工作，很多不太清楚，謝謝大家?guī)兔Γ? |
回復此樓» 收錄本帖的淘帖專輯推薦
 生物信息 |
» 猜你喜歡
307求調(diào)劑
已經(jīng)有9人回復
-
288求調(diào)劑
已經(jīng)有15人回復
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有5人回復
-
294求調(diào)劑材料與化工專碩
已經(jīng)有14人回復
-
0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有18人回復
-
265求調(diào)劑
已經(jīng)有3人回復
-
296求調(diào)劑
已經(jīng)有6人回復
-
0703化學調(diào)劑
已經(jīng)有13人回復
-
085600材料與化工 求調(diào)劑
已經(jīng)有15人回復
-
生物學調(diào)劑招人！��！
已經(jīng)有4人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請教，轉(zhuǎn)錄組測序與基因組測序結(jié)果分析
已經(jīng)有4人回復
- 有木有一個基因組疾病的數(shù)據(jù)庫？
已經(jīng)有5人回復
- 測序完之后，怎么看測序結(jié)果？請指教
已經(jīng)有6人回復
- RACE拼接后擴增全長的測序結(jié)果不一樣
已經(jīng)有5人回復
- 測轉(zhuǎn)錄組遇到的一個問題 求高手解答
已經(jīng)有4人回復
- 基因組測序結(jié)果問題
已經(jīng)有7人回復
- 各位，我800bp的基因測序后有3個點突變，后續(xù)能做表達用嗎？
已經(jīng)有7人回復
- 我擴了某個物種的某個基因，測序結(jié)果比對后，我驚呆了
已經(jīng)有9人回復
- Illumina高通量RNA測序數(shù)據(jù)分析中不同標準化方法的綜合比較
已經(jīng)有6人回復
- 鑒定一個基因的特定功能需要從哪些方面開展研究才能確定？
已經(jīng)有10人回復
- 基因測序?qū)Ρ戎�，該怎樣判斷下一步怎么做�?/a>
已經(jīng)有3人回復
- 植物轉(zhuǎn)錄組測序結(jié)果
已經(jīng)有9人回復
- 急！急！急！對轉(zhuǎn)錄組測序結(jié)果進行驗證，如何選定基因
已經(jīng)有9人回復
- 轉(zhuǎn)錄組數(shù)據(jù)分析免費軟件
已經(jīng)有11人回復
- 測序結(jié)果的疑問
已經(jīng)有7人回復
- 【求助/交流】玉米基因組庫
已經(jīng)有7人回復
- 【求助/交流】primer-blast
已經(jīng)有3人回復
- 生物信息學專題－生物版，醫(yī)學版和信息科學版共同創(chuàng)建
已經(jīng)有87人回復
cyanobacter
金蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 639.8
- 散金: 50
- 紅花: 4
- 帖子: 153
- 在線: 220.6小時
- 蟲號: 1083494
- 注冊: 2010-08-27
- 專業(yè): 環(huán)境微生物學
2樓2014-08-29 22:15:34
3樓2014-08-29 22:19:20
liuderong
鐵桿木蟲 (正式寫手)
- BioEPI: 1
- 應(yīng)助: 310 (大學生)
- 金幣: 8202.1
- 紅花: 31
- 帖子: 582
- 在線: 142.1小時
- 蟲號: 1882145
- 注冊: 2012-07-07
- 性別: GG
- 專業(yè): 生物化學
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
wizardfan: 金幣+2, 感謝應(yīng)助 2014-08-30 06:32:26
cyanobacter: 金幣+90, ★★★★★最佳答案, 太感謝了，問題解決了！ 2014-09-01 10:51:14
感謝參與，應(yīng)助指數(shù) +1
wizardfan: 金幣+2, 感謝應(yīng)助 2014-08-30 06:32:26
cyanobacter: 金幣+90, ★★★★★最佳答案, 太感謝了，問題解決了！ 2014-09-01 10:51:14
|
樓主試一下-v參數(shù)吧！ -v參數(shù)和-b參數(shù)： 這兩個參數(shù)都是限制輸出結(jié)果的數(shù)量的。 -v (integer)：規(guī)定輸出中每一個query的比對列表最多顯示subject個數(shù)（即"Sequences producing significant alignments:"后面列出的subjects數(shù)目），缺省為500條。 -b (integer)：規(guī)定輸出中每個query最多顯示與多少條subject的比對條形圖（即每條query的結(jié)果中">"的個數(shù)），缺省為250條。 如果同時使用"-m 8"參數(shù)，則輸出結(jié)果中的subjects數(shù)量和"-b"參數(shù)規(guī)定的數(shù)量一致。 在database數(shù)據(jù)中能和query比上的subjects過多的時候，這兩個參數(shù)就能夠幫助我們把其中比對結(jié)果最好的一部分挑出來，屏蔽掉相對差的結(jié)果。當然有些時候我們是不希望屏蔽掉這些結(jié)果的，比如在某個大基因組的Contig數(shù)據(jù)集中統(tǒng)計一條轉(zhuǎn)座子的重復次數(shù)，就需要把"-v"和"-b"參數(shù)定的足夠大以顯示所有結(jié)果。 |
4樓2014-08-29 23:52:41
5樓2014-08-30 08:26:49
cyanobacter
金蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 639.8
- 散金: 50
- 紅花: 4
- 帖子: 153
- 在線: 220.6小時
- 蟲號: 1083494
- 注冊: 2010-08-27
- 專業(yè): 環(huán)境微生物學
6樓2014-08-30 18:31:40
7樓2014-08-31 13:05:57
零下一度lj
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 40
- 帖子: 8
- 在線: 2.2小時
- 蟲號: 6116406
- 注冊: 2017-03-23
- 專業(yè): 遺傳學與生物信息學
送紅花一朵 |
你解決這個問題了嗎？我做的是找出不同的片段，也是要找出相似度很低的片段，但是比對的結(jié)果太多，而且因為是同源，所以給出的結(jié)果都是相似度很好的片段。跪求解答 發(fā)自小木蟲Android客戶端 |
8樓2017-04-09 16:04:43
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cyanobacter 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|---|---|---|---|---|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +14 | enenenhui 2026-03-13 | 15/750 |
|
|
|
[考研] 0703化學調(diào)劑 ，六級已過，有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
|
[考研] 321求調(diào)劑 +8 | 何潤采123 2026-03-18 | 10/500 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 0703化學調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|