#!/usr/bin/perl-w
use strict;
my $file="test.fa";
open(FH,$file) or die $!;
my $G_num=0;
my $C_num=0;
my $bases=0;
while(){
if($=~/^>/){next;}
if(/G/){$G_num+=1;}
if(/C/){$C_num+=1;}
$bases+=length($);
}
close FH;
my $GC=($G_num+$C_num)/$bases;
$GC=sprintf"%0.2f",$GC;
sprintf"GC content is:$GC\n";
~
嗯,有没有运行日志?有的话,你你你看一下 嗯,或者你发给我看一下也行