Changeset 9


Ignore:
Timestamp:
Apr 22, 2010, 6:35:36 PM (10 years ago)
Author:
maarten
Message:

Fixed file names and parsing indels

Files:
4 edited

Legend:

Unmodified
Added
Removed
  • lib/Omim2dbsnp.pm

    r7 r9  
    4848                }
    4949                else {
    50                         $dbsnp = "";
     50                        $dbsnp = "\t";
    5151                        #carp( "omimid:" . $omim,
    5252                        #       " numberref:", $av, " value does not exist/n" );
  • lib/ProccesAllicVariants.pm

    r8 r9  
    1717sub trim($) {
    1818        my $string = shift;
    19         $string =~ s/^\s+//;
    20         $string =~ s/\s+$//;
     19        if ( length($string) > 0 ) {
     20                $string =~ s/^\s+//;
     21                $string =~ s/\s+$//;
     22        }
    2123        return $string;
     24}
     25
     26sub getfilename {
     27        my ( $title, $alt ) = @_;
     28        my $protein = "";
     29
     30        if ( $title =~ m/\;/ ) {
     31                if ( $title =~ /([\w\-]{2,10})\s{0,1}$/ ) {
     32                        $protein = $1;
     33
     34                }
     35        }
     36
     37        if ( $protein eq "" ) {
     38
     39                my $titlelong = $title . ";" . $alt . ";";
     40
     41                #remove breaks
     42                $titlelong =~ s/\r\n|\n|\r/;/g;
     43
     44                if ( $protein eq "") {
     45
     46                        if ( $alt =~ /^(\w{2,10})$/ ) {
     47                                $protein = $alt;
     48
     49                        }
     50
     51                        else {
     52                                if ( $titlelong =~ /[\;\:][ ]*?([\w\-]{3,10})[ ]*?\;/ ) {
     53                                        $protein = $1;
     54                                }
     55                        }
     56                }
     57        }
     58        return ( trim $protein);
    2259}
    2360
     
    2764        my $mutation = $av->aa_ori() . $av->position() . $av->aa_mut();
    2865        if ( length($mutation) == 0 ) {
    29                 $mutation = formatInDelMutation($av->additional_mutations());
    30                
     66                $mutation = formatInDelMutation( $av->additional_mutations() );
     67
    3168        }
    3269        return ($mutation);
     
    89126                        $result = $1;
    90127                }
     128
    91129                #n A-to-G transition at nucleotide 1730 in exon 13, resulting in
    92                 elsif ( $desc =~ /([ACTG]-to-[ACTG] trans[a-z]{4,9})[\w\s\-\.]{0,10}(nucleotide \d{1,5})[\w\s\-\.]{0,90}((exon|intron)\s\d{1,3})/i ) {
    93                         $result = $1." at ".$2." in ".$3;
    94                        
     130                elsif ( $desc =~
     131/([ACTG]-to-[ACTG] trans[a-z]{4,9})[\w\s\-\.]{0,10}(nucleotide \d{1,5})[\w\s\-\.]{0,90}((exon|intron)\s\d{1,3})/i
     132                  )
     133                {
     134                        $result = $1 . " at " . $2 . " in " . $3;
    95135
    96136                        #print($result);
    97                        
     137
    98138                        #T-to-C transition BLABLA in exon 30
    99139                }
     
    106146                elsif ( $desc =~ /([0-9]{1,5}[ACTG]-[ACTG] trans[a-z]{4,9})/i ) {
    107147                        $result = $1;
     148
     149                        #print($result);
     150
     151                        #T-to-C transition BLABLA in exon 30
     152                }
     153                elsif ( $desc =~ /([ACTG]{3}-to-[ACTG]{3})/i ) {
    108154                        $result = $1;
    109155
    110156                        #print($result);
    111                        
     157
    112158                        #T-to-C transition BLABLA in exon 30
    113                 }elsif($desc =~ /([ACTG]-to-[ACTG] trans[a-z]{4,9})[\w\s\-\.]{0,90}((exon|intron)\s\d{1,3})/i){
    114                         $result=$1." in ".$2;
     159                }
     160
     161                elsif ( $desc =~
     162/([ACTG]-to-[ACTG] trans[a-z]{4,9})[\w\s\-\.]{0,90}((exon|intron)\s\d{1,3})/i
     163                  )
     164                {
     165                        $result = $1 . " in " . $2;
    115166                }
    116167                elsif ( $desc =~ /([ACTG]{3} to [ACTG]{3})/ ) {
    117168                        $result = $1;
     169                }
     170                else {
     171
     172                        $result = "\t";
    118173                }
    119174
     
    136191        # add regexp for authors with  And between names
    137192
    138         push( @regexp, "\\b([A-Za-z\\-]*.[^\\)\\,]\\sand\\s[\\w-]*.)\\s*.\\(([\\d]{4})\\)" );
     193        push( @regexp,
     194                "\\b([A-Za-z\\-]*.[^\\)\\,]\\sand\\s[\\w-]*.)\\s*.\\(([\\d]{4})\\)" );
    139195
    140196        #TODO: add solution for mulitple authors (author 2009; auth3 2002 ,etc )
     
    152208        #de Vries and de Wet (1986, 1987)
    153209        push( @regexp,
    154                 "\\b([A-Za-z\\-]{3,11}\\sand\\s[\\w\\s-]{3,11})\\s*.\\(([\\d\\s,]{4,17})\\)" );
     210"\\b([A-Za-z\\-]{3,11}\\sand\\s[\\w\\s-]{3,11})\\s*.\\(([\\d\\s,]{4,17})\\)"
     211        );
    155212
    156213        #See 123456.1234
     
    163220
    164221                if ( @srcs = ( $desc =~ /$reg/g ) ) {
     222
    165223                        #print($reg);
    166224                        push( @returnarray, @srcs );
     
    192250                        my $auth = $ref->authors();
    193251                        if ( $auth =~ m/^$authAV/ ) {
    194                                 my $authY = $ref->location();
    195                                 my $tempresult="";
     252                                my $authY      = $ref->location();
     253                                my $tempresult = "";
    196254                                if ( $authY =~ m/$yearAV/ ) {
    197                                        
     255
    198256                                        $tempresult =
    199257                                          $self->{Omim2Pubmed}->getpubmed( $numb, $counter );
    200                                           if ($tempresult ne ""){
    201                                                 print($counter." ".$numb."\n");
    202                                                 $result.=$tempresult;
    203                                           }else{
    204                                                 $result.=$authAV."(".$yearAV.") ".$ref->location();
    205                                                
    206                                           }
     258                                        if ( $tempresult ne "" ) {
     259                                                $result .= "{PMID"
     260                                                  . $tempresult . ":"
     261                                                  . $authAV . " ("
     262                                                  . $yearAV . ")}";
     263                                        }
     264                                        else {
     265                                                $result .=
     266                                                  $authAV . "(" . $yearAV . ") " . $ref->location();
     267
     268                                        }
    207269                                }
    208270                        }
     
    238300
    239301sub formatAuthors {
    240         my ($authref) = shift;
     302        my ($authref)       = shift;
    241303        my @authorsAndYears = @$authref;
    242         my @formatedAuth = ();
     304        my @formatedAuth    = ();
    243305
    244306        for ( my $i = 0 ; $i < scalar(@authorsAndYears) ; $i = $i + 2 ) {
     
    260322        }
    261323        else {
    262                 $deceaseID = "decease ID not found";
     324                $deceaseID = "\t";
    263325        }
    264326        return ($deceaseID);
     
    275337#}
    276338
    277 sub formatInDelMutation{
    278                 my $mutation = shift;
    279 #               print("'".$mutation."'");
    280                 my $formatMutation="";
    281                 if  ( $mutation =~ /^1-BP\s*?(DEL)\,\s*?(\d{1,5})([A-Z])/i ) {
     339sub formatInDelMutation {
     340        my ($mutation,$desc) = @_;
     341        #create return vraiable
     342        my $formatMutation = "";
     343        #check a postion is availble
     344        if (! ($mutation =~ /,/)) {
     345                #carp("no postion\n $mutation \n");
     346                 if($desc =~ /loss of nucleotide(s | )(\d{1,5})/i ){
     347                        $mutation= $mutation.",".$2;
     348                 }
     349        }
     350        #codon check
     351        if ( $mutation =~ /,[\w ]*?(\d+)/){
     352                my $codonnumber =$1;
     353#               carp("test 1 \n");
     354                if ( $desc =~ /codon $codonnumber/){
     355#                       carp("test 2 \n");
     356                        my $nuclitdenumber=($codonnumber*3)-2;
     357                        $mutation=~s/$codonnumber/$nuclitdenumber/;
     358       
     359                }
     360        }
     361        #check indels
     362        if ( $mutation =~ /(DEL|INS).*(DEL|INS)/){
     363#               carp("INDEL 2 \n");
     364                my $indelregex="(\\d{1,5})-BP\\s*?(DEL|INS)\\/(\\d{1,5})-BP\\s*?(DEL|INS)";
     365                if($mutation=~/$indelregex/){
     366#                               carp("INDEL 3 \n");
    282367                       
    283                         $formatMutation=$2.lc($1).$3;
    284                 }elsif($mutation =~ /^(\d{1,5})-BP\s*?(DEL|INS)\,[\sNT]*?(\d{1,5})([A-Z]*.?)/i){
    285                         #correction for calculating deletions
    286                         my $cor_del=1;
    287                         if(lc($2)eq"ins"){
    288                                 $cor_del=0;
     368                        my $del=0;
     369                        my $ins=0;
     370                        #derm length of del en ins and asign to right variable
     371                        if ($2 eq "DEL"){
     372                         $del=$1;
     373                                $ins=$3;
     374                        }else{
     375                                 $del=$3;
     376                                 $ins=$1;
    289377                        }
    290                         $formatMutation=$3."_".($3+$1-$cor_del).lc($2).$4;
    291                 }elsif($mutation =~ /^(\d{1,5})([A-Z]*?)\s*?(DEL)/i){
    292                         $formatMutation=$1.lc($3).$2;
    293                 }elsif($mutation =~/(INT(\d{0,3})-(\d{0,3}) (DEL))/i){
    294                
    295                         $formatMutation="IVS".$2."+?_IVS".$3."+?".lc($4);
    296                         $mutation =~ s/INT\d{0,3}-\d{0,3} DEL//;
    297                         if(length($mutation)>5){
    298                                 $formatMutation.=",".formatInDelMutation($mutation);
     378                        my $mutationclean=$mutation;
     379                        $mutationclean=~s/$indelregex/ /;
     380#                       carp("\n\n $mutationclean \n\n");
     381                        if ( my @results=$mutationclean =~ /([ACTGN]*)(\d{1,5})([ACTGN]*)/i ) {
     382#                                               carp("INDEL 4 \n");
     383                                $formatMutation=$results[1]."_".(($results[1])+($del-1))."del".$results[0]."ins".$results[2];
    299384                        }
    300385                }
    301                 else{
    302                         $formatMutation=$mutation;
    303                 }
    304                
    305                
    306                 return ($formatMutation);
     386        }
     387        #               print("'".$mutation."'");
     388       
     389        if ( $mutation =~ /^1-BP\s*?(DEL)\,\s*?(\d{1,5})([A-Z])/i ) {
     390
     391                $formatMutation = $2 . lc($1) . $3;
     392        }
     393        elsif ( $mutation =~
     394                /^(\d{1,5})-BP\s*?(DEL|INS)\,[\sNT]*?(\d{1,5})([A-Z]*.?)/i )
     395        {
     396
     397                #correction for calculating deletions
     398                my $cor_del = 1;
     399                if ( lc($2) eq "ins" ) {
     400                        $cor_del = 0;
     401                }
     402                $formatMutation = $3 . "_" . ( $3 + $1 - $cor_del ) . lc($2) . $4;
     403        }
     404        elsif ( $mutation =~ /^(\d{1,5})([A-Z]*?)\s*?(DEL)/i ) {
     405                $formatMutation = $1 . lc($3) . $2;
     406        }
     407        elsif ( $mutation =~ /(INT(\d{0,3})-(\d{0,3}) (DEL))/i ) {
     408
     409                $formatMutation = "IVS" . $2 . "+?_IVS" . $3 . "+?" . lc($4);
     410                $mutation =~ s/INT\d{0,3}-\d{0,3} DEL//;
     411                if ( length($mutation) > 5 ) {
     412                        $formatMutation .= "," . formatInDelMutation($mutation);
     413                }
     414        }
     415        elsif($formatMutation eq "") {
     416                $formatMutation = $mutation;
     417        }
     418
     419        return ($formatMutation);
    307420}
    308421
     
    311424        my $result      = "";
    312425        my @authAndYear = @$authref;
    313         my ( $author, $year, $pubmedid ) = "";
     426        my ( $author, $year, $pubmedidref ) = "";
    314427        for ( my $i = 0 ; $i < scalar(@authAndYear) ; $i = $i + 2 ) {
    315                 $author   = $authAndYear[$i];
    316                 $year     = $authAndYear[ $i + 1 ];
    317                 $pubmedid = $self->getPubmedID( $refsref, $author, $year, $numb );
    318                 $result .= "{PMID" . $pubmedid . ":" . $author . " (" . $year . ")}";
     428                $author = $authAndYear[$i];
     429                $year   = $authAndYear[ $i + 1 ];
     430                $result .= $self->getPubmedID( $refsref, $author, $year, $numb );
     431
     432        }
     433        if ( $result eq "" ) {
     434                $result = "\t";
    319435        }
    320436
     
    327443        my $result = "";
    328444        foreach my $av (@avs) {
    329 #print($av->number() . "\t");
    330 if ($av->title() ne "REMOVED FROM DATABASE"){
    331                 my $mutation        = getMutation( $av );
    332                
    333                 my $authorsAndYears = getAuthor( $av->description() );
    334                 $authorsAndYears = formatAuthors($authorsAndYears);
    335                 my $articles =
    336                   $self->createAuthorsOutput( $refsref, $authorsAndYears, $numb );
    337        
    338                 #lool up allelic variant in dbsnp and return accesionnumber 
    339                 my $dbsnp = $self->{Omim2DBsnp}->getdbsnp( $numb, $av->number() );
    340                                
    341                 #print(Dumper($authorsAndYears));
    342                
    343                 my $desc = $av->description();
    344                
    345                 $desc =~ tr/\n/ /;
    346                 if ( length($desc) < 60 ) {
    347 
    348                         my $pubmedid = "test";
    349 
    350                         # $self->getPubmedID( $refsref, $authAV, $yearAV, $numb );
    351                         $result .=
    352                             $numb . "\t"
    353                           . $av->number() . "\t"
    354                           . $av->title() . "\t"
    355                           . $mutation . "\t"
    356                           . $desc
    357                           . $articles . "\n";
    358                 }
    359                 else {
     445
     446                #print($av->number() . "\t");
     447                if ( $av->title() ne "REMOVED FROM DATABASE" ) {
     448                        my $mutation = getMutation($av);
     449
     450                        my $authorsAndYears = getAuthor( $av->description() );
     451                        $authorsAndYears = formatAuthors($authorsAndYears);
     452                        my $articles =
     453                          $self->createAuthorsOutput( $refsref, $authorsAndYears, $numb );
     454
     455                        #lool up allelic variant in dbsnp and return accesionnumber
     456                        my $dbsnp = $self->{Omim2DBsnp}->getdbsnp( $numb, $av->number() );
     457
     458                        #print(Dumper($authorsAndYears));
     459
     460                        my $desc = $av->description();
     461
     462                        $desc =~ tr/\n/ /;
     463
     464                        #                       if ( length($desc) < 60 ) {
     465                        #
     466                        #                               my $pubmedid = "test";
     467                        #
     468                        #                               # $self->getPubmedID( $refsref, $authAV, $yearAV, $numb );
     469                        #                               $result .=
     470                        #                                   $numb . "\t"
     471                        #                                 . $av->number() . "\t"
     472                        #                                 . $av->title() . "\t"
     473                        #                                 . $mutation . "\t"
     474                        #                                 . $desc
     475                        #                                 . $articles . "\n";
     476                        #                       }
     477                        #                       else {
    360478
    361479                        my $deceaseID = getDeceaseID($desc);
     
    370488                          . $deceaseID . "\t"
    371489                          . $mutation . "\t"
    372                           . $dbsnp. "\t"
     490                          . $dbsnp . "\t"
    373491                          . $articles . "\t"
    374492
    375 
    376493                          . $nucleotidemutation . "\n";
    377                 }
    378 }
     494
     495                        #                       }
     496                }
    379497        }
    380498
  • omimparser.pl

    r7 r9  
    1313
    1414        use Carp;
    15 sub trim($)
    16 {
    17         my $string = shift;
    18         $string =~ s/^\s+//;
    19         $string =~ s/\s+$//;
    20         return $string;
    21 }
    2215
    2316
    2417
    25 sub getfilename{
    26   my (  $title,$alt)=@_;
    27   my $protein="";
    28 
    29   if($title =~ m/\;/){
    30     $title=~/(\w{2,10})\s{0,1}$/;
    31     $protein=$1;
    32   }else{
    33  
    34   my $titlelong=$title.$alt;
    35   #remove breaks
    36   $titlelong=~tr/\n/ /;
    37   #split on ;
    38   my @splitstring=split(/\;/, $titlelong);
    39   # print Dumper(@splitstring);
    40  
    41   #get the protein on second position of array
    42    $protein=($splitstring[1]);
    43   }
    44  
    45   return(trim$protein);
    46 
    47 
    48 }
    4918
    5019
     
    7443$PAV-> addOmim2pubmed($omim2pubmed);
    7544$PAV->  addOmim2dbsnp($dbsnp);
     45open (FILE, ">output/filename_errors.txt") or die $!;
     46
    7647  while ( my $omim_entry = $omim_parser->next_phenotype() ) {
    7748    # This prints everything.
     
    133104  #print Dumper(@avs);
    134105        if (scalar(@avs)>0){
    135               my $geneid=getfilename($title,$alt);
    136               open (FILE, ">output/".$geneid.".txt") or die $!;
     106                     
     107          my $geneid=ProccesAllicVariants::getfilename($title,$alt);
    137108         
    138               print FILE ($PAV->proccesAllicVariants(\@avs,\@refs,$numb));
    139               close (FILE);
     109             # print FILE ($numb."\t".$geneid."||".$title."<_".$alt."\n");
    140110
    141 }
     111               
     112                if( -e "output/".$geneid.".txt"){
     113                        print FILE ("\nWARNING FILENAME ALREADY EXITS  with name: $geneid and $numb\n $title \n $alt");
     114               
     115                }
     116                if($geneid eq ""){
     117                        if(! -e "output/".$numb.".txt"){
     118                                $geneid=$numb;
     119                        }
     120                }
     121                             
     122      open (CSV, ">output/".$geneid.".txt") or die $!;
     123       
     124      print FILE ($PAV->proccesAllicVariants(\@avs,\@refs,$numb));
     125             close (CSV);
     126
     127
    142128      # Array of Bio::Phenotype::OMIM::OMIMentryAllelicVariant objects.
    143129      # class Bio::Phenotype::OMIM::OMIMentryAllelicVariant
     
    160146    # do something ...
    161147  }
     148                }
     149close (FILE);
  • test/test.pl

    r8 r9  
    1010use Omim2dbsnp;
    1111use Bio::Phenotype::OMIM::OMIMparser;
    12 use Test::More tests => 62;
     12use Test::More tests =>71;
    1313
    1414#lib "../lib";
     
    197197
    198198is(
     199        ProccesAllicVariants::extractnuclitdemutation( " bedu drg r g r syhrt eyh5tr r wefess AAG-to-AAC change bla", "75" ),
     200        "AAG-to-AAC",
     201        "nucleotide as tripple AAG-to-AAC"
     202);
     203
     204
     205is(
    199206        ProccesAllicVariants::extractnuclitdemutation("on on 1 chromosome in the proband was a T-to-C transition predicted to result in a cys1417-to-arg substitution.The mutation was located in exon 30","10"),"T-to-C transition in exon 30","T-to-C transition BLABLA in exon 30");
    200207#acta .00010
     
    225232print("\ntesting format of deletions inserts and other awkward mutations\n");
    226233is(
    227         ProccesAllicVariants::formatInDelMutation("1-BP DEL, 1764G") ,
     234        ProccesAllicVariants::formatInDelMutation("1-BP DEL, 1764G","") ,
    228235        "1764delG",
    229236        "1-BP DEL, 1764G"
    230237);
    231 
    232 is(
    233         ProccesAllicVariants::formatInDelMutation( "2-BP DEL, 3283TC") ,
     238#600046.0004
     239is(
     240        ProccesAllicVariants::formatInDelMutation("3-BP DEL","In a French-Canadian family with familial high density lipoprotein deficiency (604091) previously reported by Marcil et al. (1995), Brooks-Wilson et al. (1999) found a 3-bp deletion that resulted in loss of nucleotides 2017-2019 and deletion of a leucine at position 633, which is conserved in mouse and C. elegans") ,
     241        "2017_2019del",
     242        "3-BP DEL"
     243);
     244
     245#tected a change in codon 747 in
     246#601199.0008
     247
     248is(
     249        ProccesAllicVariants::formatInDelMutation( "2-BP DEL/1-BP INS, CCC747TC","athyroidism (239200), Pearce et al. (1995) detected a change in codon 747 in ") ,
     250        "2239_2240delCCCinsTC",
     251        "2239_2240delCCCinsTC"
     252);
     253
     254is(
     255        ProccesAllicVariants::formatInDelMutation( "2-BP DEL/1-BP INS,747TC","athyroidism (239200), Pearce et al. (1995) detected a change in codon 747 in ") ,
     256        "2239_2240delinsTC",
     257        "2239_2240delinsTC"
     258);
     259
     260is(
     261        ProccesAllicVariants::formatInDelMutation( "2-BP DEL, 3283TC","") ,
    234262        "3283_3284delTC",
    235263        "3283_3284delTC"
     
    237265
    238266is(
    239         ProccesAllicVariants::formatInDelMutation( "40-BP DEL, NT1294") ,
     267        ProccesAllicVariants::formatInDelMutation( "40-BP DEL, NT1294","") ,
    240268        "1294_1333del",
    241269        "1294_1333del"
     
    243271
    244272is(
    245         ProccesAllicVariants::formatInDelMutation( "1675A DEL") ,
     273        ProccesAllicVariants::formatInDelMutation( "1675A DEL","") ,
    246274        "1675delA",
    247275        "1675delA"
    248276);
    249277is(
    250         ProccesAllicVariants::formatInDelMutation( "1-BP INS, 1135A") ,
     278        ProccesAllicVariants::formatInDelMutation( "1-BP INS, 1135A","") ,
    251279        "1135_1136insA",
    252280        "1135_1136insA"
     
    254282
    255283is(
    256         ProccesAllicVariants::formatInDelMutation( "110-BP INS/14-BP DEL") ,
     284        ProccesAllicVariants::formatInDelMutation( "110-BP INS/14-BP DEL","") ,
    257285        "110-BP INS/14-BP DEL",
    258286        "110-BP INS/14-BP DEL"
    259287);
    260288is(
    261         ProccesAllicVariants::formatInDelMutation( "INT12-14 DEL") ,
     289        ProccesAllicVariants::formatInDelMutation( "INT12-14 DEL","") ,
    262290        "IVS12+?_IVS14+?del",
    263291        "IVS12+?_IVS14+?del"
     
    265293
    266294is(
    267         ProccesAllicVariants::formatInDelMutation( "INT12-14 DEL, INT16-31 DEL") ,
     295        ProccesAllicVariants::formatInDelMutation( "INT12-14 DEL, INT16-31 DEL"," nothing") ,
    268296        "IVS12+?_IVS14+?del,IVS16+?_IVS31+?del",
    269297        "IVS12+?_IVS14+?del,IVS16+?_IVS31+?del"
     
    271299#BRCA1.0033: 5-BP INS, NT3171. Als de nucleotiden in de tekst gespecificeerd worden, zouden we die graag terug zien, bijv. 3171_3172insAGTCG.
    272300is(
    273         ProccesAllicVariants::formatInDelMutation( "5-BP INS, NT3171") ,
     301        ProccesAllicVariants::formatInDelMutation( "5-BP INS, NT3171","") ,
    274302        "",
    275303        ""
    276304);
    277305is(
    278         ProccesAllicVariants::formatInDelMutation( "") ,
     306        ProccesAllicVariants::formatInDelMutation( "","") ,
    279307        "",
    280308        ""
     
    282310
    283311is(
    284         ProccesAllicVariants::formatInDelMutation( "") ,
     312        ProccesAllicVariants::formatInDelMutation( "","") ,
    285313        "",
    286314        ""
     
    288316
    289317is(
    290         ProccesAllicVariants::formatInDelMutation( "") ,
     318        ProccesAllicVariants::formatInDelMutation( "","") ,
    291319        "",
    292320        ""
     
    303331
    304332is( ( $PAV->getPubmedID( \@refs, "Ilkovski", "2001", $numb ) ),
    305         "11333380", "get omim record " );
     333        "{PMID11333380:Ilkovski (2001)}", "get omim record " );
    306334        #               $self->createAuthorsOutput( $refsref,$authorsAndYears, $numb  );
    307335         $authandyear=["Ilkovski", "2001"];
     
    322350@avs = $omim_entry->each_AllelicVariant();
    323351
    324 print(Dumper(@refs));
    325 is( ( $PAV->getPubmedID( \@refs, "Pitteloud", "2004", $numb ) ),
    326         "Pitteloud(2004) J. Clin. Endocr. Metab. 89 1053-1058 (2004)", "get omim record without omim number " );
    327        
     352
    328353
    329354is( ( $PAV->getPubmedID( \@refs, "Brown", "1989", $numb ) ),
    330         "2563196", "get omim record Brown " );
     355        "{PMID2563196:Brown (1989)}", "get omim record Brown " );
    331356               
    332 is( ( $PAV->getPubmedID( \@refs, "Lumbroso", "1996", $numb ) ),
    333         "8626869", "get omim record Lumbroso 100" );
    334        
    335         is( ( $PAV->getPubmedID( \@refs, "McPhaul", "1991", $numb ) ),
    336         "2010552", "get omim record  McPhaul 112" );
    337        
    338                 is( ( $PAV->getPubmedID( \@refs, "Migeon", "1981", $numb ) ),
    339         "6947233", "get omim record Migeon 120" );
    340        
    341                 is( ( $PAV->getPubmedID( \@refs, "Mononen", "2000", $numb ) ),
    342         "11103816", "get omim record  Mononen 123" );
    343        
    344                         is( ( $PAV->getPubmedID( \@refs, "Montie", "2009", $numb ) ),
    345         "19279159", "get omim record  Montie 124" );
    346        
    347 #                               is( ( $PAV->getPubmedID( \@refs, "", "", $numb ) ),
    348 #       "", "get omim record  " );
     357
    349358       
    350359is( ( $PAV->getPubmedID( \@refs, "Mooney", "2003", $numb ) ),
    351         "12682377", "get omim record Mooney " );
     360        "{PMID12682377:Mooney (2003)}", "get omim record Mooney " );
    352361       
    353362is( ( $PAV->getPubmedID( \@refs, "Trifiro", "1989", $numb ) ),
    354         "Trifiro (1989) (Abstract) Am. J. Hum. Genet. 45 (suppl.): A225, 1989", "get omim record whitout PMID" );
     363        "Trifiro(1989) (Abstract) Am. J. Hum. Genet. 45 (suppl.) A225 (1989)", "get omim record whitout PMID" );
    355364                       
    356365       
     
    370379
    371380is( ( $PAV->getPubmedID( \@refs, "Allikmets", "1997", $numb ) ),
    372         "Pitteloud(2004) J. Clin. Endocr. Metab. 89 1053-1058 (2004)", "get omim record without omim number " );
     381        "{PMID9295268:Allikmets (1997)}{PMID9054934:Allikmets (1997)}", "get omim record with 2 entries " );
     382#THYMIDINE PHOSPHORYLASE: TYMP<_;;TP;;
     383is(ProccesAllicVariants::getfilename("THYMIDINE PHOSPHORYLASE: TYMP",";;TP;;"),"TYMP","get name");
     384        is(ProccesAllicVariants::getfilename("LMBR1 DOMAIN-CONTAINING PROTEIN 1: LMBRD1",";;LMBD1;;"),"LMBRD1","get name");
     385                        is(ProccesAllicVariants::getfilename("V-MYC AVIAN MYELOCYTOMATOSIS VIRAL-RELATED ONCOGENE, NEUROBLASTOMA-DERIVED;","MYCN\n;; afas sfsdas "),"MYCN","get name");
     386        is(ProccesAllicVariants::getfilename("LMBR1 DOMAIN-CONTAINING PROTEIN 1: LMBRD1",";;LMBD1;;"),"LMBRD1","get name");
     387        is(ProccesAllicVariants::getfilename("TRANSIENT RECEPTOR POTENTIAL CATION CHANNEL, SUBFAMILY M, MEMBER 6;","TRPM6"),"TRPM6","get name");
Note: See TracChangeset for help on using the changeset viewer.