From 708ed6b63e65ae5ad1d12a38803fe89fa2d6d44d Mon Sep 17 00:00:00 2001
From: Ewan Birney <birney@sanger.ac.uk>
Date: Mon, 17 Apr 2000 14:08:15 +0000
Subject: [PATCH] added clone stats to get out coding potential for clones
---
scripts/clonestats.pl | 47 +++++++++++++++++++++++++++++++++++++++++++
1 file changed, 47 insertions(+)
create mode 100755 scripts/clonestats.pl
diff --git a/scripts/clonestats.pl b/scripts/clonestats.pl
new file mode 100755
index 0000000000..653914fe31
--- /dev/null
+++ b/scripts/clonestats.pl
@@ -0,0 +1,47 @@
+#!/usr/local/lib/perl
+
+
+my $locator = "Bio::EnsEMBL::DBSQL::Obj/host=obi-wan;dbname=ensembl;user=ensro;";
+$db = Bio::EnsEMBL::DBLoader->new($locator);
+@clones = $db->get_all_Clone_id();
+
+
+foreach my $clone_id ( @clones ) {
+
+ eval {
+ my $clone = $db->get_Clone($clone_id);
+ @genes = $clone->get_all_Genes();
+ $genenumber = scalar @genes;
+ my @exons;
+
+ foreach $g ( @genes ) {
+ push(@exons,$g->each_unique_Exon);
+ }
+
+ $exonnumber= scalar @exons;
+ $length = 0;
+ foreach $exon ( @exons ) {
+ $length += $exon->lenght;
+ }
+
+ $gene{$clone_id} = $genenumber;
+ $exon{$clone_id} = $exonnumber;
+ $exonl{$clone_id} = $length;
+ };
+
+ if( $i > 20 ) {
+ last;
+ }
+
+ $i++;
+}
+
+@clones = sort { $exon{$a} <=> $exon{$b} } keys %exon;
+
+foreach $clone_id (@clones) {
+ print "$clone_id\t$exon{$clone_id}\t$gene{$clone_id}\t$exonl{$clone_id}\n";
+}
+
+
+
+
--
GitLab