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

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

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

金蟲 (小有名氣)
木蟲

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

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考研] 生物學(xué)調(diào)劑招人。! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|