open(FILE,"C:/Users/hp/Desktop/UNIQ.fasta");

$i = 0;


while($line = <FILE>)

{ chomp $line;

if($line =~/(\>gi\|\d{1,})/)

{

if($i >= 1)

{

$hash{$def[$i]} = $list[$i];

}

$i++;

$def[$i] = $line;

}

else

{

$list[$i] .= $line;

}



} 


open(FILE,">C:/Users/hp/Desktop/Amino_acid_count_data.txt");

print 
FILE "DEFINITION\tLENGTH\tG_RATIO\tE_RATIO\tD_RATIO\tA_RATIO\tV_RATIO\tR_RATIO\

tS_RATIO\tK_RATIO\tN_RATIO\tT_RATIO\tM_RATIO\tI_RATIO\tQ_RATIO\tH_RATIO\tP_RATI

O\tL_RATIO\tW_RATIO\tC_RATIO\tY_RATIO\tF_RATIO\tMiss_RATIO\n";

$o = 0;

foreach $k (keys %hash) {

$o++;

    #print "$k => $hash{$k}\n";

Passref($k,$hash{$k},$o);

		}


		
		
		
		


		
		###########################sub boundary################
		sub Passref

{

my($define,$seq,$count) = @_;

chomp ($define);

print $define,"\n";

print  $seq,"\n";

print  $count,"\n";

$len = length ($seq);

#print "\n$len\n";

$count_G = ( $seq =~ tr/G//);

#print "$count_G\n";

$count_E = ( $seq =~ tr/E//);

$count_D = ( $seq =~ tr/D//);

$count_A = ( $seq =~ tr/A//);

$count_V = ( $seq =~ tr/V//);

$count_R = ( $seq =~ tr/R//);

$count_S = ( $seq =~ tr/S//);

$count_K = ( $seq =~ tr/K//);

$count_N = ( $seq =~ tr/N//);

$count_T = ( $seq =~ tr/T//);

$count_M = ( $seq =~ tr/M//);

$count_I = ( $seq =~ tr/I//);

$count_Q = ( $seq =~ tr/Q//);

$count_H = ( $seq =~ tr/H//);

$count_P = ( $seq =~ tr/P//);

$count_L = ( $seq =~ tr/L//);

$count_W = ( $seq =~ tr/W//);

$count_C = ( $seq =~ tr/C//);

$count_Y = ( $seq =~ tr/Y//);

$count_F = ( $seq =~ tr/F//);

$count_Miss = ( $seq =~s/[^YCWLPHQIMTNKSRVADEGF]//ig);

#print "$count_Miss\n";



$ratio_G = ($count_G/$len); 

#print "\n$ratio_G\n"; 

$ratio_E = ($count_E/$len);

$ratio_D = ($count_D/$len);

$ratio_A = ($count_A/$len);

$ratio_V = ($count_V/$len);

$ratio_R = ($count_R/$len);

$ratio_S = ($count_S/$len);

$ratio_K = ($count_K/$len);

$ratio_N = ($count_N/$len);

$ratio_T = ($count_T/$len);

$ratio_M = ($count_M/$len);

$ratio_I = ($count_I/$len);

$ratio_Q = ($count_Q/$len);

$ratio_H = ($count_H/$len);

$ratio_P = ($count_P/$len);

$ratio_L = ($count_L/$len);

$ratio_W = ($count_W/$len);

$ratio_C = ($count_C/$len); 

$ratio_Y = ($count_Y/$len);

$ratio_F = ($count_F/$len);

$ratio_Miss = ($count_Miss/$len);

open(FH,">>C:/Users/hp/Desktop/Amino_acid_count_data.txt");

print 

FH "$define\t$len\t$ratio_G\t$ratio_E\t$ratio_D\t$ratio_A\t$ratio_V\t$ratio_R\t

$ratio_S\t$ratio_K\t$ratio_N\t$ratio_T\t$ratio_M\t$ratio_I\t$ratio_Q\t$ratio_H\

t$ratio_P\t$ratio_L\t$ratio_W\t$ratio_C\t$ratio_Y\t$ratio_F\t$ratio_Miss\n";

#print "$define\t$len\t$ratio_G\t$ratio_E\t$ratio_D\t$ratio_A\t$ratio_V\t$ratio

_R\t$ratio_S\t$ratio_K\t$ratio_N\t$ratio_T\t$ratio_M\t$ratio_I\t$ratio_Q\t$rati

o_H\t$ratio_P\t$ratio_L\t$ratio_W\t$ratio_C\t$ratio_Y\t$ratio_F\t$ratio_Miss\n"
;


}