fasta序列宽度转换

#!/usr/bin/perl -w
use strict;
############################
#
#使用方法：windows的dos界面或者linux下：perl file_name
###########################
my ($input_file,$head,$line,$width,%seq);
print "输入文件名称，并按Enter\n";
$input_file=;
  chomp$input_file;
open IN,"<$input_file";
open OUT,">result";
   while (){
            chomp;
       if (/^>/){$head=$_;}else{$seq{$head}.=$_;}
       }
 print "程序负责单条序列转换成特定长度的序列\n";      
 print "输入期望的序列宽度，并按一下Enter键\n";     
              
           $width=;
            chomp$width;
foreach $line(sort keys %seq){
                     chomp$line;
             
              # $seq{$line}=~s/(?<=\w{30})(?                $seq{$line}=~s/((\w{$width})+?)(\w{0})/$1\n$3/gi;  #折叠  这一块需要优化，最好的方法是用环顾(?=|?<=)，比较一下就知道速度的差异了
               
               print OUT "$line\n$seq{$line}\n";
           }
           
       exit;