我来一个菜鸟水准的
[user132@mu01 perl]$ cat test_3
>salind_6404
CCTGCGGAAGATCGGCACTAGAATAGCCAGAACCGTTTCTCTGAGGCTTCCGGCCTTCCC
TCCCACTAATAATTCTGAGG
>Rosalind_5959
CCATCGGTAGCGCATCCTTAGTCCAATTAAGTCCCTATCCAGGCGCTCCGCCGAAGGTCT
ATATCCATTTGTCAGCAGACACGC
>Rosalind_0808
CCACCCTCGTGGTATGGCTAGGCATTCAGGAACCGGAGAACGCTTCAGACCAGCCCGGAC
TGGGAACCTGCGGGCAGTAGGTGGAAT
[user132@mu01 perl]$ cat test_3.pl
[Perl] 纯文本查看 复制代码 #!/usr/bin/perl
open IN,"<$ARGV[0]" or die $!;
$/ = ">";<IN>;
my $max_gc;
my $max_id;
while(<IN>){
s/\r?\n>//;
my($id,$seq) =split /\n/,$_,2;
my $length = length $seq;
my $gc =($seq =~ s/C/C/g) + ($seq =~ s/G/G/g);
my $GC =$gc/$length*100;
if($max_gc < $GC){
$max_gc = $GC;
$max_id = $id;
}
# print "id$id\tlength$length\tgc$gc\tmax_gc$max_gc\n";
}
print "$max_id\t$max_gc\n";
close IN;
[user132@mu01 perl]$ perl test_3.pl test_3
Rosalind_0808 59.5505617977528
|