 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2295 | 回復: 8 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
shuifeng1988金蟲 (小有名氣)
木蟲
|
[求助]
perl腳本求助，我想將一萬多個文件從fasta格序，改為phylip4格式
|
||
| 求助，如何利用bioperl批量修改序列格式!我在iny文件夾中有一萬多個.fasta文件，我想將它們修改序列格式為phylip4，并存在out文件夾中！我電腦上已經(jīng)裝好bioperl，請問怎么寫腳本？？ |
回復此樓» 收錄本帖的淘帖專輯推薦
 perl語言專欄 |
» 猜你喜歡
山東省面上項目限額評審
已經(jīng)有4人回復
-
一志愿北京化工大學070300 學碩336求調(diào)劑
已經(jīng)有5人回復
-
生物學071000 329分求調(diào)劑
已經(jīng)有4人回復
-
一志愿華中科技大學071000，求調(diào)劑
已經(jīng)有4人回復
-
生物學一志愿985，分數(shù)349求調(diào)劑
已經(jīng)有4人回復
-
生物學調(diào)劑
已經(jīng)有4人回復
-
求調(diào)劑院校信息
已經(jīng)有4人回復
-
085600材料與化工306
已經(jīng)有4人回復
-
286求調(diào)劑
已經(jīng)有10人回復
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有9人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助有機牛人，有什么辦法可以強有力的促進還原消除啊~~
已經(jīng)有13人回復
- 求助：篩菌中有關(guān)試劑和培養(yǎng)基的問題
已經(jīng)有1人回復
- 同個任務同樣的節(jié)點數(shù)VASP4.6能算，而VASP5.2不能
已經(jīng)有9人回復
- rtl格式的文件怎么打開？
已經(jīng)有3人回復
- 參考文獻格式求助
已經(jīng)有6人回復
- 【分享】處理DOSCAR文件的腳本和程序
已經(jīng)有19人回復
- 【求助】namd、vmd推薦個學習修改和編寫腳本網(wǎng)站和地方吧
已經(jīng)有3人回復
- 【求助】用vasp4.6計算GaAs的能帶結(jié)構(gòu)
已經(jīng)有18人回復
shuifeng1988
金蟲 (小有名氣)
木蟲
- 應助: 0 (幼兒園)
- 金幣: 999
- 紅花: 1
- 帖子: 60
- 在線: 42.5小時
- 蟲號: 1076748
- 注冊: 2010-08-18
- 性別: GG
- 專業(yè): 基因組學
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
【答案】應助回帖
★ ★
jjdg(金幣+2): 感謝應助 2011-06-20 11:55:32
shuifeng1988(金幣+5): 2011-10-12 13:30:19
jjdg(金幣+2): 感謝應助 2011-06-20 11:55:32
shuifeng1988(金幣+5): 2011-10-12 13:30:19
|
話說，我粗略的看了下上述兩種格式，簡單的腳本就能解決，用不著什么bioperl的大炮吧， google 大神給出的兩個貌似能解決問題的鏈接： python 版本 http://miss.ieph.net/archives/970 awk 版本 http://www.unix.com/shell-progra ... -phylip-format.html |
2樓2011-06-19 22:45:55
【答案】應助回帖
★ ★
微塵、夢想(金幣+2): 謝謝參與應助！ 2011-06-20 19:38:40
shuifeng1988(金幣+5): 2011-10-12 13:30:27
微塵、夢想(金幣+2): 謝謝參與應助！ 2011-06-20 19:38:40
shuifeng1988(金幣+5): 2011-10-12 13:30:27
|
google大神還提供了這個：http://myhits.isb-sib.ch/cgi-bin/reformat 看起來是cgi寫的，就當perl腳本咯~ |
3樓2011-06-20 12:01:36
shuifeng1988
金蟲 (小有名氣)
木蟲
- 應助: 0 (幼兒園)
- 金幣: 999
- 紅花: 1
- 帖子: 60
- 在線: 42.5小時
- 蟲號: 1076748
- 注冊: 2010-08-18
- 性別: GG
- 專業(yè): 基因組學
★ ★
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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 材料學學碩080502 337求調(diào)劑-一志愿華中科技大學 +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 中南大學化學學碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] A區(qū)線材料學調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學666 2026-03-16 | 7/350 |
|
|
|
[考研] 0703化學調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 一志愿福大288有機化學，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿985，本科211，0817化學工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 0703化學調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|