Hello!

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

Note that updates take place every 10 minutes, commits may not be seen immediately.
Fix a couple of typos
authorelserj <elserj@localhost>
Thu, 1 Dec 2016 03:53:32 +0000 (03:53 +0000)
committerelserj <elserj@localhost>
Thu, 1 Dec 2016 03:53:32 +0000 (03:53 +0000)
svn path=/; revision=649

interactome_scripts/find_species.pl

index 097079e729d679e61c0d7631eaff7df6700961fd..98c0e30665d9e52a304cdfcf8ba0471173db2ffc 100755 (executable)
@@ -88,7 +88,7 @@ sub find_species {
                $species = "Elaphocordyceps_paradoxa";
        }elsif ($temp =~ /Epichloe_festucae/) {
                $species = "Epichloe\_festucae";
-       }elsif ($temp = /Erinaceus_europaeus/) {
+       }elsif ($temp =~ /Erinaceus_europaeus/) {
                $species = "Erinaceus_europaeus";
        }elsif ($temp =~ /Escherichia\_coli/) {
                $species = "Escherichia_coli";
@@ -144,8 +144,8 @@ sub find_species {
                $species = "Mus_musculus";
        }elsif ($temp =~ /Musa/) {
                $species = "Musa_acuminata";
-       }elsif ($temp =~ /Mustela_putorious.furo/) {
-               $species = "Mustela_putorious.furo";
+       }elsif ($temp =~ /Mustela_putorius\.furo/) {
+               $species = "Mustela_putorius.furo";
        }elsif ($temp =~ /Nectria\_haematococca/) {
                $species = "Nectria_haematococca";
        }elsif ($temp =~ /Nelumbo/) {
@@ -358,7 +358,7 @@ sub find_gene {
                my ($gene_info,$type,$gc,$gene_length) = split(/\s/, $gene_header);
                my ($gene_number,$gene_id) = split(/::/, $gene_info);
                my ($trinity_id, $m_number) = split(/\|/, $gene_id);
-               $gene = $trinity_id; #???
+               $gene = $gene_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;