| Server IP : 79.136.114.73 / Your IP : 3.14.251.87 Web Server : Apache/2.4.7 (Ubuntu) PHP/5.5.9-1ubuntu4.29 OpenSSL/1.0.1f System : Linux b8009 3.13.0-170-generic #220-Ubuntu SMP Thu May 9 12:40:49 UTC 2019 x86_64 User : www-data ( 33) PHP Version : 5.5.9-1ubuntu4.29 Disable Function : pcntl_alarm,pcntl_fork,pcntl_waitpid,pcntl_wait,pcntl_wifexited,pcntl_wifstopped,pcntl_wifsignaled,pcntl_wexitstatus,pcntl_wtermsig,pcntl_wstopsig,pcntl_signal,pcntl_signal_dispatch,pcntl_get_last_error,pcntl_strerror,pcntl_sigprocmask,pcntl_sigwaitinfo,pcntl_sigtimedwait,pcntl_exec,pcntl_getpriority,pcntl_setpriority, MySQL : ON | cURL : ON | WGET : ON | Perl : ON | Python : ON | Sudo : ON | Pkexec : ON Directory : /usr/bin/ |
Upload File : |
#!/usr/bin/perl
use strict; use warnings;
my %Feature = (
'Exon' => 'CDS',
'Einit' => 'CDS',
'Eterm' => 'CDS',
'Esngl' => 'CDS',
);
print "##gff-version 3\n";
my $seq;
my %H;
while (<>) {
if (/^>(\S+)/) {
print "#region $1\n";
$seq = $1;
} else {
my @f = split;
if (@f == 4) {
my $strand = $f[1] < $f[2] ? '+' : '-';
if ($strand eq '-') {($f[1], $f[2]) = ($f[2], $f[1])}
print join("\t", $seq, 'snap', $Feature{$f[0]}, $f[1], $f[2], '.',
$strand, "Name=$f[3]"), "\n";
} elsif (@f == 9) {
print join("\t", $seq, 'snap', $Feature{$f[0]}, $f[1], $f[2], '.',
$f[3], "Name=$f[8]"), "\n";
} else {die "input does not appear to be ZFF"}
}
}
__END__