作者rom1025 (περ'ιπεπια)
看板Perl
标题Re: [问题]有关资料撷取部分的一个小问题
时间Tue May 13 20:24:58 2008
※ 引述《anysomething (努力一点比在实在！)》之铭言：
: ※ 引述《LinusChen (骗自己)》之铭言：
: : 我的做法是直接判断该行是否是 ">"开头, 不是就直接输出...
: 谢谢LinusChen提供的建议，判断">"的部分我不太会写
: 现在有一个问题是，我是要
: >ENSP00000332127
: MGLLAFRDVALEFSPEEWECLDPAQRSLYRDVMLENYRNLISLGEDSFNMQFLFHSLAMS
: KPELIICLEARKEPWNVNTEKTARHSVLSSYLTEDILPEQGLQVSFQKVMLRRYERCCLE
: KLRLRNDWEIPCEDVLASPLPSAMILSFLRPPQKQKHVKPTEPIQSKTLQWL
试试看这个可不可以用
open (IN, 'X:\filename.txt');
while ($line = <IN>){
chomp($line);
if ($line=~ /^>(ENSP.*?)\s/) {
$ENSP = $1;
open (OUT, ">>$ENSP.fasta");
print OUT ">$ENSP\n";
}else{
print OUT "$line\n";
}
}
close IN;
close OUT;
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 59.104.147.233