| 5 | 1/1 | 返回列表 |
| 查看: 2293 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
shuifeng1988金蟲 (小有名氣)
木蟲
|
[求助]
perl腳本求助,我想將一萬多個(gè)文件從fasta格序,改為phylip4格式
|
||
| 求助,如何利用bioperl批量修改序列格式!我在iny文件夾中有一萬多個(gè).fasta文件,我想將它們修改序列格式為phylip4,并存在out文件夾中!我電腦上已經(jīng)裝好bioperl,請問怎么寫腳本?? |
perl語言專欄 |

金蟲 (著名寫手)
|
google大神還提供了這個(gè):http://myhits.isb-sib.ch/cgi-bin/reformat 看起來是cgi寫的,就當(dāng)perl腳本咯~ |

金蟲 (小有名氣)
|
話說,我粗略的看了下上述兩種格式,簡單的腳本就能解決,用不著什么bioperl的大炮吧, google 大神給出的兩個(gè)貌似能解決問題的鏈接: python 版本 http://miss.ieph.net/archives/970 awk 版本 http://www.unix.com/shell-progra ... -phylip-format.html |
金蟲 (小有名氣)
木蟲
|
#/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 ] |

金蟲 (小有名氣)
木蟲
|
#/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 ] |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 一志愿211,0703化學(xué)310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 311求調(diào)劑 +13 | 冬十三 2026-03-15 | 14/700 |
|
|
[考研] 354求調(diào)劑 +6 | Tyoumou 2026-03-18 | 9/450 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|