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