1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59
|
#! /usr/bin/perl
use strict;
open (G, ">resu.txt");
open (E, "./holoe.data") or die " impossible d'ouvrir le fichier holoe\n";
open (F, "./M.data/M.data5");
my $compteur =0;
my $gpe ="";
my @hml =<E>;
my $gpe = "";
my $i=0;
my $test="";
my $long =@hml;
open (DATA, "Mdata/M.data5") or die "Can't open file M.data5 $!";
my @mdata = <DATA>;
chomp(@mdata);
my $mdata = join("\t", @mdata);
my @ge=split(/\/\//, $mdata);
my $ge ="";
foreach $ge(@ge)
{
if($ge ne "")
{
chomp($ge);
$ge =~ s/'/\\'/g;
$_ = $ge;
my ($unigId) = /ID\s+([^\t]*)/;
my ($geNam) = /TITLE\s+([^\t]*)/;
my ($geSym) = /GENE\s+([^\t]*)/;
my ($Var) = /LOCUSLINK\s+([^\t]*)/;
my ($sequences) = /ID\s+.*?(SEQUENCE.*)$/;
print "$unigId, $geSym, $Var\n";
if ($Var ne "")
{print G "variable1 \t$Var\n"; #mon prog tourne jusque la!
for ($i=0; $i<=$long ; $i++)
{ $test=~/^([^\t]*)\t[^\t]*\t([^\t]*)\t[^\t]*\t[^\t]*\t[^\t]*\t.*/;
print $test;
if ($Var eq $2)
{
$gpe = $1;
print G "Gpe\t$gpe\n";
}
}
}
}
} |
Partager