脚本程序

原始序列

运行程序：perl calculateDNA_length_GCcontent.pl

第一步合并结果

第二步计算结果

——————————————————————————————————————————————————

#!/usr/bin/perl -w

use strict;

my $content="";

open In,"sampleDNA.fasta";

open Out,">sampleDNA.fasta.out";

while(<In>){

chomp;

if(/^>/){

if($content ne ""){

print Out "$content\n";

$content="";}

print Out "$_\n";

}

else{$content.=$_;}

}

print Out "$content\n";

close In;close Out;

open In,"sampleDNA.fasta.out";

open Out,">sampleDNAresult.out";

while(<In>){

chomp;

my $GC=0;my $length=0;my $GC_count=0;

if(/^>(\S+)/){

print Out "$1\t";}

else {$_=~s/\s//g;

$length=length($_);

while($_=~/c|g/ig){$GC++;};

$GC_count=$length?($GC/$length):0;

$GC_count*=100;

printf Out "%d\t%.2f",$length,$GC_count;

printf Out "%%\n";}

}

close In;close Out;