 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2303 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
shuifeng1988金蟲 (小有名氣)
木蟲
|
[求助]
perl腳本求助，我想將一萬多個(gè)文件從fasta格序，改為phylip4格式
|
||
| 求助，如何利用bioperl批量修改序列格式!我在iny文件夾中有一萬多個(gè).fasta文件，我想將它們修改序列格式為phylip4，并存在out文件夾中！我電腦上已經(jīng)裝好bioperl，請問怎么寫腳本？？ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 perl語言專欄 |
» 猜你喜歡
一志愿北化求調(diào)劑
已經(jīng)有4人回復(fù)
-
085601材料工程找調(diào)劑
已經(jīng)有6人回復(fù)
-
318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿雙一流機(jī)械285分求調(diào)劑
已經(jīng)有4人回復(fù)
-
332求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
已經(jīng)有8人回復(fù)
-
375求調(diào)劑
已經(jīng)有3人回復(fù)
-
南京大學(xué)化學(xué)調(diào)劑
已經(jīng)有7人回復(fù)
-
070300化學(xué)354求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 求助有機(jī)牛人，有什么辦法可以強(qiáng)有力的促進(jìn)還原消除啊~~
已經(jīng)有13人回復(fù)
- 求助：篩菌中有關(guān)試劑和培養(yǎng)基的問題
已經(jīng)有1人回復(fù)
- 同個(gè)任務(wù)同樣的節(jié)點(diǎn)數(shù)VASP4.6能算，而VASP5.2不能
已經(jīng)有9人回復(fù)
- rtl格式的文件怎么打開？
已經(jīng)有3人回復(fù)
- 參考文獻(xiàn)格式求助
已經(jīng)有6人回復(fù)
- 【分享】處理DOSCAR文件的腳本和程序
已經(jīng)有19人回復(fù)
- 【求助】namd、vmd推薦個(gè)學(xué)習(xí)修改和編寫腳本網(wǎng)站和地方吧
已經(jīng)有3人回復(fù)
- 【求助】用vasp4.6計(jì)算GaAs的能帶結(jié)構(gòu)
已經(jīng)有18人回復(fù)
shuifeng1988
金蟲 (小有名氣)
木蟲
- 應(yīng)助: 0 (幼兒園)
- 金幣: 999
- 紅花: 1
- 帖子: 60
- 在線: 42.5小時(shí)
- 蟲號: 1076748
- 注冊: 2010-08-18
- 性別: GG
- 專業(yè): 基因組學(xué)
★ ★
jjdg:編輯內(nèi)容 2011-10-12 04:45
jjdg:編輯內(nèi)容 2011-10-12 04:45
jjdg(金幣+2): 感謝分享 2011-10-12 04:46:43
shuifeng1988: 回帖置頂 2011-12-27 16:51:38
jjdg:編輯內(nèi)容 2011-10-12 04:45
jjdg:編輯內(nèi)容 2011-10-12 04:45
jjdg(金幣+2): 感謝分享 2011-10-12 04:46:43
shuifeng1988: 回帖置頂 2011-12-27 16:51:38
|
#/usr/bin/perl use Bio::AlignIO; my @name = glob("*.fasta" ); foreach $name (@name) { $in = Bio::AlignIO -> new(-file => "$name", '-format' => 'fasta'); $out = Bio::AlignIO -> new(-file => ">..\\sequence03\\$name", '-format' => 'phylip'); while (my $aln = $in -> next_aln()) {$out ->write_aln($aln); } } [ Last edited by jjdg on 2011-10-12 at 04:45 ] |
5樓2011-10-11 21:34:12
xioooli
金蟲 (小有名氣)
- 程序強(qiáng)帖: 1
- 應(yīng)助: 2 (幼兒園)
- 金幣: 744
- 帖子: 82
- 在線: 33小時(shí)
- 蟲號: 940951
- 注冊: 2010-01-11
- 性別: GG
- 專業(yè): 藥化
【答案】應(yīng)助回帖
★ ★
jjdg(金幣+2): 感謝應(yīng)助 2011-06-20 11:55:32
shuifeng1988(金幣+5): 2011-10-12 13:30:19
jjdg(金幣+2): 感謝應(yīng)助 2011-06-20 11:55:32
shuifeng1988(金幣+5): 2011-10-12 13:30:19
|
話說，我粗略的看了下上述兩種格式，簡單的腳本就能解決，用不著什么bioperl的大炮吧， google 大神給出的兩個(gè)貌似能解決問題的鏈接： python 版本 http://miss.ieph.net/archives/970 awk 版本 http://www.unix.com/shell-progra ... -phylip-format.html |
2樓2011-06-19 22:45:55
huycwork
金蟲 (著名寫手)
- 程序強(qiáng)帖: 22
- 應(yīng)助: 0 (幼兒園)
- 金幣: 953
- 散金: 663
- 紅花: 8
- 沙發(fā): 13
- 帖子: 1080
- 在線: 264.1小時(shí)
- 蟲號: 1257243
- 注冊: 2011-04-06
- 專業(yè): 金融學(xué)
【答案】應(yīng)助回帖
★ ★
微塵、夢想(金幣+2): 謝謝參與應(yīng)助！ 2011-06-20 19:38:40
shuifeng1988(金幣+5): 2011-10-12 13:30:27
微塵、夢想(金幣+2): 謝謝參與應(yīng)助！ 2011-06-20 19:38:40
shuifeng1988(金幣+5): 2011-10-12 13:30:27
|
google大神還提供了這個(gè)：http://myhits.isb-sib.ch/cgi-bin/reformat 看起來是cgi寫的，就當(dāng)perl腳本咯~ |
3樓2011-06-20 12:01:36
shuifeng1988
金蟲 (小有名氣)
木蟲
- 應(yīng)助: 0 (幼兒園)
- 金幣: 999
- 紅花: 1
- 帖子: 60
- 在線: 42.5小時(shí)
- 蟲號: 1076748
- 注冊: 2010-08-18
- 性別: GG
- 專業(yè): 基因組學(xué)
jjdg:編輯內(nèi)容 2011-10-12 04:46
|
#/usr/bin/perl use Bio::AlignIO; my @name = glob("*.fasta" ); foreach $name (@name) { $in = Bio::AlignIO -> new(-file => "$name", '-format' => 'fasta'); $out = Bio::AlignIO -> new(-file => ">..\\sequence03\\$name", '-format' => 'phylip'); while (my $aln = $in -> next_aln()) {$out ->write_aln($aln); } } [ Last edited by jjdg on 2011-10-12 at 04:46 ] |
4樓2011-10-11 21:33:56
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)354求調(diào)劑 +5 | 101次希望 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 291求調(diào)劑 +4 | Y-cap 2026-03-29 | 4/200 |
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-29 | 4/200 |
|
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +4 | 26要上岸 2026-03-28 | 4/200 |
|
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 307求調(diào)劑 +8 | 超級伊昂大王 2026-03-24 | 9/450 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|