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
|
########### PARTIE II : Analyse des resultats BLAST##########
## question (1) : calcul du nombre du resultat BLAST####
open(INFILE ,'<','resultat_Blast1000.tab');
my $line = <INFILE>; my @lines= <INFILE>; my @keyz;my $compte=0;
foreach (@lines){ if($_=~m/clone/g) { push (@keyz,$_); } }
$compte=$#keyz+1;
###### Question 1 parie II (le nombre de resultts blst)####
print"le nombre totale des resultats BLAST est :$compte\n";
#print "@keyz\n";
#########print "@lines\n"; #### attention de confondre lines pour le blast en pluriel###
my @enfin; my $score; my $especeKEGG ; my $evalue ; my $fragment;
my @ressign;my $num=0;
my @dax_score;
my @dax_evalue;
foreach my $dax(@lines)
{
my @enfin = split(m/\t/s,$dax);
#print "clone : $enfin[0] espece KEGG:$enfin[1] evalue : $enfin[10] score: $enfin[11] \n";
#print "$enfin[0]\n";
push(@dax_evalue,$enfin[10]);
push(@dax_score,$enfin[11]);
}
#print "@dax_score\n";
#print "@dax_evalue\n\n\n";##### probleme d'affichage de resultat~~~~
my $comteur=0;
my $prov;
my $papa=0;
foreach $prov(0..$#dax_score){
if (($dax_score[$prov]>50)&&($dax_evalue[$prov]<10**-5)){
$comteur++;
my @blastsign;
my @evalsign;
push (@blastsign,$dax_score[$prov]);
push (@evalsign,$dax_evalue[$prov]);
print "scores des resultats blast significtifs:@blastsign:@evalsign";
}
}
#######question 2 partie II (le nombre de resultats blast significatifs ) #####
#print "le nombre de resultat BLAST significatifs est :$comteur\n"; |
