if($readingframe2 =~ /MHGR/){
        print "$1\n";
        ($protein)= $readingframe2 =~ /MHGR/;