view snpEff_2_1a/scripts/file2GenomeName.pl @ 0:f8eaa3f8194b default tip

Uploaded snpEff_v2_1a_core.tgz from Pablo Cingolani
author greg
date Fri, 20 Apr 2012 14:47:09 -0400
parents
children
line wrap: on
line source

#!/usr/bin/perl

#-------------------------------------------------------------------------------
#
# Transform file name to genome short name
#
#-------------------------------------------------------------------------------

$file = $ARGV[0];
print "$file";

$base = `basename $file`;
chomp $base;
print "\t$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); }
$full = "$gen.$short";
print "\t$full\t$gen\t$short";


print "\n";