作者anysomething (努力一点比在实在！)
看板Perl
标题Re: [问题]有关资料撷取部分的一个小问题
时间Tue May 13 14:33:21 2008
※ 引述《LinusChen (骗自己)》之铭言：
: 我的做法是直接判断该行是否是 ">"开头, 不是就直接输出...
谢谢LinusChen提供的建议，判断">"的部分我不太会写
现在有一个问题是，我是要
>ENSP00000332127
MGLLAFRDVALEFSPEEWECLDPAQRSLYRDVMLENYRNLISLGEDSFNMQFLFHSLAMS
KPELIICLEARKEPWNVNTEKTARHSVLSSYLTEDILPEQGLQVSFQKVMLRRYERCCLE
KLRLRNDWEIPCEDVLASPLPSAMILSFLRPPQKQKHVKPTEPIQSKTLQWL
单独输出一个档案
我将code改写了
while(<IN>)
{
chomp;
$_ =~ /^(>.*)/
$ENSP = $1;
print "$ENSP", "\n";
if ($_ =~ /^>(.*)/ )
{
$openENSP = $1;
}
else
{
$sequence = $_;
open(FILE, "+>/home/philippe/CutPeptide/$openENSP.fas");
print FILE ">$ENSP\n$sequence\n";
close FILE;
}
}
结果输出变成
>ENSP00000332127
KLRLRNDWEIPCEDVLASPLPSAMILSFLRPPQKQKHVKPTEPIQSKTLQWL
中间的部分被忽略了!
请问是我回圈部分出了问题吗??
想请高手帮我解答一下!
谢谢!!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.109.55.234