1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
#!/usr/bin/perl
#-------------------------------------------------------------------------------
#
# Transform file name to genome short name
#
#-------------------------------------------------------------------------------
$file = $ARGV[0];
$base = `basename $file`;
chomp $base;
if( $base =~ /(.*?)\.(.*)\..?dna\.(.*)\.fa\.gz/ ) { ($gen, $short) = ($1, $2); }
elsif( $base =~ /(.*?)\.(.*)\.gtf\.gz/ ) { ($gen, $short) = ($1, $2); }
elsif( $base =~ /(.*?)\.(.*)\.pep\.all\.fa\.gz/ ) { ($gen, $short) = ($1, $2); }
if( $short =~ /(.*)\.\d+/ ) { $id = $1; }
$full = "$gen.$short";
print "$file\t$base\t$full\t$gen\t$short\t$id\t$gen.$id\n";
|
