- or download this
#!/usr/bin/perl -w
use strict;
die 'Usage: ctpA [DIR]' if @ARGV > 1;
my $dir = shift || '.';
- or download this
use File::Basename;
open local(*OUT), "> $dir/pA.cat" or die $!;
- or download this
for (glob "${dir}/*.seq" ) { - or download this
local $/;
open local(*IN), "< $_" or die $!;
my $seq = <IN>;
close IN or die $!;
- or download this
$seq =~ s/\s//g; - or download this
$seq =~ /([AN]*)$/;
printf OUT "%d %s.seq\n", length $1, basename $_;
}
- or download this
close OUT or die $!;