新手请求帮助！！！！！！

Patton007

最近在处理某些文件，如下：
>;ctg1
GAGCTCCACCGCGGTGGCGGCCGCTCTAGAACTAGTGGATCCCCCGGGCTGCAGGAATTC
GCGGCCGCGTCGACTTTATTCAAGCCTGTTTCCATCCAGACCATGTGGTGGGTTTTATTC
AAGCCCAGTTCCATACAGACCATGTGGTGGGTTTTATTCAAGCCTGTTTCCATCCAGACC
ATGTGGTGGGTTTTATTCAAGCCTGTTTCCATCCAGACCATGTGGTAGGTTTTATTCAAG
CCTGTTTCCATCCAGACTGAGTGGTGGGTTTTATTCAAGCCTGTTTCCATCCAGACCATG
TGGTAGGTTTTATTCAAGCCTGTTTCCATCCAGACCATGTGGTGGGTTTTATTCAAGCCT
GGTTCCATACAGACCATGTGGTGGGTTTTATTCAAGCCTGGTTCCATCCAGACCATGTGG
TAGGTTTTATTCAAGCCTGTTTCCATCCAGACCATGTGGTGGGTTTTATTCAAGCCTAGT
TCCATACAGACCATGTGGTGGGTTTTATTCAAGCCTGTTTCCATCCAGACCATGTGGTGG
GTTTTATTCAAGCCTTGTTCCATACAGACCATGTGGTGG
>;ctg2
GCTGGAGCTCCACCGCGGTGGCGGCCGCTCTAGAACTAGTGGATCCCCCGGGCTGCAGGA
ATTCGCGGCCGCGTCGACCTCGAGTTTTTTGGACAGAAATGTTGTTCAATGGCAACTTTG
AAGGATATAGCAGTAAGACTACGGTCCGTCAAAAACATACAGAAAATCACTGCGTCCATG
AAGATGGTTTCAGCTTCAAAATTTGCTAGGGCTGAACGTGACCTTCGTCCGGCCAGACCT
TATGGTCAAGGTGCTCGAGCATTTTATAAGACCGCAGACCTAGTACCCAGTGAGCTTAAA
CCTAGTGAGAATGATTTACTGATAGCAATCACATCAGACAGAGGTCTATGTGGTGCAGCT
AATAGTAGTATTGTGAAGGCGATCAGAGCTCAGCTGCGCAGTAATCCAGAGCTAGAAAAG
TCATGTAAATTACTCTTGGTTGGTGACAAGTCTCGTGCAATGATGCTTCGTCAGTACCGT
AACATGTTTTTGATGACTGTAAATGAAGTTGGCAAGAAAGCACCCACATTCGAGGATGCC
TCGACAATAGCCGAAGCAATCATCTCTACTGACTTCAAGTTTGACAAGGCTGTGATGTTT
TACAATACATTCAAGACAGT
.......
序列每行都是60个“ACGT”（由这四个字母组成的），每条完整的序列都是以">;"开始。我目的是想统计序列的长度，并把一定长度的序列挑出来，但因为很多条这样格式的序列，所以只能把序列读到一行中再进行序列长度的统计。
因为我是新手，写了好多次总是在读完一行60个字母的序列之后转入下一行是出现问题，所以请各位高手给予解答，谢谢高手指教！！！！！！

apile

open(FF,"<xxx.dat");
while(<FF>;){
   if(/^>;/){
      my $data=<FF>;;
      $data=~s/\012\015//g;
      my $len = length($data);
      print "length of ".$.." is ".$len."\n";
   }
}
close(FF);

Patton007

谢谢apile 兄！！！

可是我这句 $data=~s/\012\015//g;
不是很懂，“\012\015”代表什么意思？为什么要用空格代替？
刚开始看perl，不是很懂。
谢谢！

apile

用
=~s/[\012\015]//g;
可能比較好..
\012＝\n
\015=\r
(應該沒錯吧)
因為windows中的斷行...等於
\n\r有時候會等於\r\n..
所以避免問題...先將他們消掉 ...
chomp只會消掉\n...不會消掉 \r..
所以常常你在linux下open windows下建立的doc..
會看到一堆^M...

Patton007

非常感谢！

seneca

脱氧核糖核酸的碱基对序列？

Patton007

是的，呵呵。

bloves

希望搞生物的朋友都搞perl，