my $sample_cds_genotype = $prefixe ."_cds_genotype";
my $sample_cds_RT = $prefixe ."_cds_RT";
my $sample_cds_R1 = $prefixe ."_cds_R1";
my $sample_cds_R2 = $prefixe ."_cds_R2";

my %hash_sample_cds_genotype;
my %hash_sample_cds_RT;
my %hash_sample_cds_R1;
my %hash_sample_cds_R2;
my %hash_geno_group;
my %hash_read_p1_count;
my %hash_read_p2_count;

#alimentation du hash des groupe genotypique connu  pour les echatillons
open(GENOTYPED,"<$list_geno") or die ("Cannot open $list_geno\n");
while (defined(my $_ = <GENOTYPED>))
	{
		my @geno_group_infos = split(m/\t/,$_);
		next if ($_ =~ m/^ECHANTILLON/);
		chomp($_);
		my $sample_nm = $geno_group_infos[0];
		my $sample_geno_group = $geno_group_infos[1];
		$hash_geno_group{$sample_geno_group}{"sample"}=$sample_nm;		
	}
close(GENOTYPED);

print Dumper(\%hash_geno_group);

#String des fichiers inputs
foreach my $fichier (@fichiers)
{
	my $filename = $fichier;
	my $sample_name;


	#alimentation des hashs de resultat avec cet echatillon

	open(RESULTAT,"<$filename") or die ("Cannot open $filename\n");

	while (defined(my $_ = <RESULTAT>))
	{
		$sample_name = $filename."\t";
		#my @attribut_name = split(m/out_/,$sample_name);
		#$sample_name = $attribut_name[1]."\t";
		my @resutat_infos = split(m/\t/,$_);
		next if ($_ =~ m/^CDS/);		
		chomp($_);
		my $cds = $resutat_infos[0] ."\t";
		my $reads_p1 =$resutat_infos[1]."\t";
		my $reads_p2 =$resutat_infos[2]."\t";
		my $count_reads_p1 =$resutat_infos[3];
		my $count_reads_p2 =$resutat_infos[4];
		my $count_reads_ambigu =$resutat_infos[5];
		my $count_reads_classe =$resutat_infos[6];
		my $reads_totaux =$resutat_infos[8]."\t";
		$hash_sample_cds_RT{$cds}{$sample_name} = $reads_totaux ;
		$hash_sample_cds_R1{$cds}{$sample_name} = $reads_p1 ;
		$hash_sample_cds_R2{$cds}{$sample_name} = $reads_p2 ;
		$hash_read_p1_count{$cds}{$sample_name} = $count_reads_p1;
		$hash_read_p2_count{$cds}{$sample_name} = $count_reads_p2;
		foreach my $group(keys(%hash_geno_group))
		{
			print $group."\n";
			if (exists($hash_geno_group{$group}{"sample"}{$sample_name}))
			{
				$hash_geno_group{$group}{$cds}{"ambigu"} += $count_reads_ambigu;
				$hash_geno_group{$group}{$cds}{"rc"} += $count_reads_classe;
				$hash_geno_group{$group}{"sample"}{$sample_name}=($count_reads_p1+$count_reads_p2);
			}
		}
	}

	close(RESULTAT);
	print "Finished\n";
}