1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
|
#!/usr/bin/perl -w
$filename = shift @ARGV;
open (IN, $filename) || die "failed to open '$filename'";
@outfields = ("BlastProDom", "Coil", "Gene3D", "HMMPanther",
"HMMPfam", "HMMSmart", "HMMTigr", "FPrintScan", "ScanRegExp",
"ProfileScan", "superfamily", "Seg", "SignalP", "TmHmm");
while (<IN>) {
@fields = split (/\t/);
$text{$fields[3]} .= $_;
}
close IN;
open (OUT, ">$filename") || die "failed to open '$filename'";
foreach $f (@outfields) {
if(defined($text{$f})) {
print OUT "$text{$f}";
}
$text{$f} = "";
}
foreach $f (sort(keys(%text))) {
if($text{$f} ne "") {
print OUT "$text{$f}";
}
}
close OUT;
|