 |
 |
 |
| 24小時熱門版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 2534 | 回復(fù): 7 | |||
cyanobacter金蟲 (小有名氣)
|
[求助]
高通量數(shù)據(jù)做本地blast時，如何將97%以上的比對結(jié)果全部顯示出來？ 已有3人參與
|
|
各位蟲友，我有一批16S的高通量數(shù)據(jù)，然后我現(xiàn)在想做一個本地blast，把我手上的一部分序列（isolate.fasta）與高通量數(shù)據(jù)進(jìn)行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的文件，然后就開始進(jìn)行本地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ī)兔Γ? |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 生物信息 |
» 猜你喜歡
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
081700化工學(xué)碩調(diào)劑
已經(jīng)有3人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有4人回復(fù)
-
329求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一）
已經(jīng)有3人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有14人回復(fù)
-
320求調(diào)劑0856
已經(jīng)有3人回復(fù)
-
307求調(diào)劑
已經(jīng)有9人回復(fù)
-
288求調(diào)劑
已經(jīng)有15人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請教，轉(zhuǎn)錄組測序與基因組測序結(jié)果分析
已經(jīng)有4人回復(fù)
- 有木有一個基因組疾病的數(shù)據(jù)庫？
已經(jīng)有5人回復(fù)
- 測序完之后，怎么看測序結(jié)果？請指教
已經(jīng)有6人回復(fù)
- RACE拼接后擴增全長的測序結(jié)果不一樣
已經(jīng)有5人回復(fù)
- 測轉(zhuǎn)錄組遇到的一個問題 求高手解答
已經(jīng)有4人回復(fù)
- 基因組測序結(jié)果問題
已經(jīng)有7人回復(fù)
- 各位，我800bp的基因測序后有3個點突變，后續(xù)能做表達(dá)用嗎？
已經(jīng)有7人回復(fù)
- 我擴了某個物種的某個基因，測序結(jié)果比對后，我驚呆了
已經(jīng)有9人回復(fù)
- Illumina高通量RNA測序數(shù)據(jù)分析中不同標(biāo)準(zhǔn)化方法的綜合比較
已經(jīng)有6人回復(fù)
- 鑒定一個基因的特定功能需要從哪些方面開展研究才能確定？
已經(jīng)有10人回復(fù)
- 基因測序?qū)Ρ戎�，該怎樣判斷下一步怎么做�?/a>
已經(jīng)有3人回復(fù)
- 植物轉(zhuǎn)錄組測序結(jié)果
已經(jīng)有9人回復(fù)
- 急！急！急！對轉(zhuǎn)錄組測序結(jié)果進(jìn)行驗證，如何選定基因
已經(jīng)有9人回復(fù)
- 轉(zhuǎn)錄組數(shù)據(jù)分析免費軟件
已經(jīng)有11人回復(fù)
- 測序結(jié)果的疑問
已經(jīng)有7人回復(fù)
- 【求助/交流】玉米基因組庫
已經(jīng)有7人回復(fù)
- 【求助/交流】primer-blast
已經(jīng)有3人回復(fù)
- 生物信息學(xué)專題－生物版，醫(yī)學(xué)版和信息科學(xué)版共同創(chuàng)建
已經(jīng)有87人回復(fù)
cyanobacter
金蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 639.8
- 散金: 50
- 紅花: 4
- 帖子: 153
- 在線: 220.6小時
- 蟲號: 1083494
- 注冊: 2010-08-27
- 專業(yè): 環(huán)境微生物學(xué)
2樓2014-08-29 22:15:34
3樓2014-08-29 22:19:20
liuderong
鐵桿木蟲 (正式寫手)
- BioEPI: 1
- 應(yīng)助: 310 (大學(xué)生)
- 金幣: 8202.1
- 紅花: 31
- 帖子: 582
- 在線: 142.1小時
- 蟲號: 1882145
- 注冊: 2012-07-07
- 性別: GG
- 專業(yè): 生物化學(xué)
【答案】應(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é)果。當(dāng)然有些時候我們是不希望屏蔽掉這些結(jié)果的，比如在某個大基因組的Contig數(shù)據(jù)集中統(tǒng)計一條轉(zhuǎn)座子的重復(fù)次數(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)境微生物學(xué)
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è): 遺傳學(xué)與生物信息學(xué)
送紅花一朵 |
你解決這個問題了嗎？我做的是找出不同的片段，也是要找出相似度很低的片段，但是比對的結(jié)果太多，而且因為是同源，所以給出的結(jié)果都是相似度很好的片段。跪求解答 發(fā)自小木蟲Android客戶端 |
8樓2017-04-09 16:04:43
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cyanobacter 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 12/600 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
[考研] 287求調(diào)劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|