Hello!

To see the file structure, click on "tree".

Note that updates take place every 10 minutes, commits may not be seen immediately.
Add more beavers
authorelserj <elserj@localhost>
Wed, 28 Dec 2016 20:50:08 +0000 (20:50 +0000)
committerelserj <elserj@localhost>
Wed, 28 Dec 2016 20:50:08 +0000 (20:50 +0000)
svn path=/; revision=654

interactome_scripts/find_species.pl

index 327d902f521f238d3a69d07abba4743ba196dd0b..5dcb80d63f1770a3ec5129b9df0a67f205885e5c 100755 (executable)
@@ -58,6 +58,10 @@ sub find_species {
                $species = "Castor_canadensis.binpacker";
        }elsif ($temp =~ /Castor_canadensis.velvet/) {
                $species = "Castor_canadensis.velvet";
+       }elsif ($temp =~ /Castor_canadensis.brent.proteins/) {
+               $species = "Castor_canadensis.brent.proteins";
+       }elsif ($temp =~ /Castor_canadensis.brent.transcripts/) {
+               $species = "Castor_canadensis.brent.transcripts";
        }elsif ($temp =~ /Cavia_porcellus/) {
                $species = "Cavia_porcellus";
        }elsif ($temp =~ /Chlamy/) {
@@ -376,6 +380,12 @@ sub find_gene {
                my ($m_id, $g_id, $temp) = split(/\s/, $gene_header);
                my ($root_id, $m_number) = split(/\|/, $m_id);
                $gene = "$root_id" . "_" . "$m_number";
+       }elsif ($species eq "Castor_canadensis.brent.proteins") {
+               my ($scaffold, $type, $coords) = split(/\s/, $gene_header);
+               $gene = $scaffold;
+       }elsif ($species eq "Castor_canadensis.brent.transcripts") {
+               my ($scaffold, $type, $coords) = split(/\s/, $gene_header);
+               $gene = $scaffold;
        }elsif ($species eq "Cavia_porcellus") {
                my ($isoform, $type, $scaffold, $temp_gene, $transcript, $gene_biotype, $transcript_biotype, $gene_symbol, $description) = split(/\s/, $gene_header);
                $gene = $isoform;