From: <she...@us...> - 2003-12-19 15:03:11
|
Update of /cvsroot/gmod/schema/chado/load/bin In directory sc8-pr-cvs1:/tmp/cvs-serv28420 Modified Files: genbank2gff3.PLS Log Message: Added a GFF 3 header to the output so it will be parsed correctly by Bio::DB::GFF Index: genbank2gff3.PLS =================================================================== RCS file: /cvsroot/gmod/schema/chado/load/bin/genbank2gff3.PLS,v retrieving revision 1.2 retrieving revision 1.3 diff -C2 -d -r1.2 -r1.3 *** genbank2gff3.PLS 18 Dec 2003 16:41:58 -0000 1.2 --- genbank2gff3.PLS 19 Dec 2003 15:03:08 -0000 1.3 *************** *** 1,3 **** ! #!perl -w use Config; --- 1,3 ---- ! #!/usr/bin/perl use Config; *************** *** 59,62 **** --- 59,72 ---- #my @top_sfs = $unflattener->get_SeqFeatures; my @skipped_lines; + + + # print a GFF header + my $date = localtime; + my $sname = $seq->accession || $seq->display_name; + my $end = length $seq->seq; + print OUT "##gff-version 3\n", + "##date $date\n", + "##sequence-region ",join(' ', $sname, 1, $end),"\n"; + foreach my $sf (@sfs) { *************** *** 239,243 **** --- 249,256 ---- } } + !NO!SUBS! + + close OUT or die "Can't close $file: $!"; chmod 0755, $file or die "Can't reset permissions for $file: $!\n"; |