| 9 | 1/1 | 返回列表 |
| 查看: 2291 | 回復(fù): 8 | |||
shuifeng1988金蟲(chóng) (小有名氣)
木蟲(chóng)
|
[求助]
perl腳本求助,我想將一萬(wàn)多個(gè)文件從fasta格序,改為phylip4格式
|
| 求助,如何利用bioperl批量修改序列格式!我在iny文件夾中有一萬(wàn)多個(gè).fasta文件,我想將它們修改序列格式為phylip4,并存在out文件夾中!我電腦上已經(jīng)裝好bioperl,請(qǐng)問(wèn)怎么寫(xiě)腳本?? |
perl語(yǔ)言專(zhuān)欄 |

金蟲(chóng) (小有名氣)
|
話說(shuō),我粗略的看了下上述兩種格式,簡(jiǎn)單的腳本就能解決,用不著什么bioperl的大炮吧, google 大神給出的兩個(gè)貌似能解決問(wèn)題的鏈接: python 版本 http://miss.ieph.net/archives/970 awk 版本 http://www.unix.com/shell-progra ... -phylip-format.html |
金蟲(chóng) (著名寫(xiě)手)
|
google大神還提供了這個(gè):http://myhits.isb-sib.ch/cgi-bin/reformat 看起來(lái)是cgi寫(xiě)的,就當(dāng)perl腳本咯~ |

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

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

金蟲(chóng) (小有名氣)
木蟲(chóng)

木蟲(chóng) (小有名氣)
|
#/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); } } |
金蟲(chóng) (小有名氣)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856 271求調(diào)劑 +3 | 展信悅_ 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +7 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 7/350 |
|
|
[考研] 302求調(diào)劑 +11 | 呼呼呼。。。。 2026-03-17 | 11/550 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +6 | 梨花珞晚風(fēng) 2026-03-17 | 6/300 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-17 | 8/400 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶(hù)573181 2026-03-20 | 5/250 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|