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 beaver
authorelserj <elserj@localhost>
Fri, 23 Dec 2016 01:09:34 +0000 (01:09 +0000)
committerelserj <elserj@localhost>
Fri, 23 Dec 2016 01:09:34 +0000 (01:09 +0000)
svn path=/; revision=653

interactome_scripts/find_species.pl

index ea74fc30567f801e55a9a868beacb7dded4006dc..327d902f521f238d3a69d07abba4743ba196dd0b 100755 (executable)
@@ -54,6 +54,10 @@ sub find_species {
                $species = "Castor_canadensis";
        }elsif ($temp =~ /Castor_canadensis.no.isoforms/) {
                $species = "Castor_canadensis.no.isoforms";
+       }elsif ($temp =~ /Castor_canadensis.binpacker/) {
+               $species = "Castor_canadensis.binpacker";
+       }elsif ($temp =~ /Castor_canadensis.velvet/) {
+               $species = "Castor_canadensis.velvet";
        }elsif ($temp =~ /Cavia_porcellus/) {
                $species = "Cavia_porcellus";
        }elsif ($temp =~ /Chlamy/) {
@@ -364,6 +368,14 @@ sub find_gene {
        }elsif ($species eq "Castor_canadensis.no.isoforms") {
                my ($gene_id,$m_number, $g_number) = split(/\|/, $gene_header);
                $gene = $gene_id;
+       }elsif ($species eq "Castor_canadensis.binpacker") {
+               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.velvet") {
+               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 "Cavia_porcellus") {
                my ($isoform, $type, $scaffold, $temp_gene, $transcript, $gene_biotype, $transcript_biotype, $gene_symbol, $description) = split(/\s/, $gene_header);
                $gene = $isoform;
@@ -769,6 +781,17 @@ sub find_gene_synonym {
                my ($gene_number,$gene_id) = split(/::/, $gene_info);
                my ($trinity_id, $m_number) = split(/\|/, $gene_id);
                $synonym = "$gene_number,$m_number";
+       }elsif ($species eq "Castor_canadensis.no.isoforms") {
+               my ($gene_id,$m_number, $g_number) = split(/\|/, $gene_header);
+               $gene = $gene_id;
+       }elsif ($species eq "Castor_canadensis.binpacker") {
+               my ($m_id, $g_id, $temp) = split(/\s/, $gene_header);
+               my ($root_id, $g_number) = split(/\|/, $m_id);
+               $synonym = "$root_id" . "_" . "$g_number";
+       }elsif ($species eq "Castor_canadensis.velvet") {
+               my ($m_id, $g_id, $temp) = split(/\s/, $gene_header);
+               my ($root_id, $g_number) = split(/\|/, $m_id);
+               $synonym = "$root_id" . "_" . "$g_number";
        }elsif ($species eq "Cavia_porcellus") {
                my ($isoform, $type, $scaffold, $temp_gene, $transcript, $gene_biotype, $transcript_biotype, $gene_symbol, $description) = split(/\s/, $gene_header);
                $temp_gene =~ s/gene\://;