#creation de la matrice des genotype selon les groupes genotypique foreach my $sample_geno_group(keys(%hash_geno_group)) { print $sample_geno_group. "\n"; foreach my $cds (keys %{$hash_geno_group{$sample_geno_group}}) { print "cds ".$cds. "\n"; next if ($cds eq "sample"); foreach my $RC(values %{$hash_geno_group{$sample_geno_group}{$cds}{"rc"}}) { print "rc ".$RC. "\n"; #verification du nombre minimun de reads pour chaque groupe genotypique if ($RC < $seuil_min_classe) { foreach my $ech_name(values @{$hash_geno_group{$sample_geno_group}{"sample"}}) { $hash_sample_cds_genotype{$cds}{$ech_name} = "NA\t" ; } delete($hash_geno_group{$sample_geno_group}); } else { foreach my $AB(values%{$hash_geno_group{$sample_geno_group}{$cds}{"ambigu"}}) { #verification du pourcentage maximun de reads ambigu pour chaque groupe genotypique if (($AB *100 / $RC) > $seuil_max_amb) { foreach my $ech_name(values @{$hash_geno_group{$sample_geno_group}{"sample"}}) { $hash_sample_cds_genotype{$cds}{$ech_name} = "NA\t" ; } delete($hash_geno_group{$sample_geno_group}); } } } }