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
|
open my $fh, '<', $fichier or die "Impossible de lire le fichier $fichier\n";
my %data;
my $i=0;
my $j=0;
while(my $ligne = <$fh>){
$i++;
chomp $ligne;
if($ligne =~ /gene.*/){
$j++;
my ($gene,$champ1,$champ2,$champ3) = split "\t", $ligne;
my $element =$champ1."\t".$champ2."\t".$champ3;
push @{$data{$gene}}, [$element];
}
}
for my $id (keys %data) {
my @array = @{$data{$id}};
$data{$id} = \@array;
}
print Dumper \%data;
for my $id (keys %data) {
print "--$id\n";
my @array = @{$data{$id}};
for my $pair (@array) {
print ("\tstart : $pair->[0]\n");
#faire la somme de chaque colonne et quand la somme de la colonne est égale à 2 mettre 1
my @hit = split("\t",$pair->[0]);
}
}
close($fh); |
Partager