求助、这段代码怎么理解呢？？

bhbperl

#/usr/bin/perl -w
my $input=shift;
$/=">";
open(IN,$input) || die "cannot open the file, please check it!\n";
open(OUT,'>temp.fa');
while (<IN>) {
        chomp;
        if ($_=~/^(.+?)\n(.*)$/s) {
                my $head=$1;
                my $seq=$2;
                $seq=~s/\n//g;
                print OUT ">$head\n"."$seq\n";
        }
}
close IN;

felix0608

是对fasta文件做一些处理。把多行展示的序列弄成一行展示回复 1# bhbperl


   

bhbperl

是的，就是其中有的部分不理解，回复 2# felix0608


   

bhbperl

例如($_=~/^(.+?)\n(.*)$/s这一行。回复 2# felix0608


   

felix0608

这个代码写的有些繁琐。less *fa|perl -pe 'chomp;s/^|$/\n/g if /^>/;'|sed 1d > out.fa 这样也能实现你的目的
回复 2# felix0608


   

felix0608

这个代码写的有些繁琐。less *fa|perl -pe 'chomp;s/^|$/\n/g if /^>/;'|sed 1d > out.fa 这样也能实现你的目的
回复 4# bhbperl


   

felix0608

回复 4# bhbperl 正则匹配而已。