※ 引述《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