統計fasta文件每條序列的長度【Perl】

#!/usr/bin/perl -w 
use strict;

my $file="myseq.fa";
open(FH,$file) or die $!;
my %hash;
my $id=0;
while(<FH>){
	chomp;
	my $temp;
	if($_=~/^>/){
		$_=~ s/^>//g;
		my $temp=$_;
		next;
	}
	$hash{$_}=length($_);
}
close FH;
foreach my $k(sort keys %hash){
	print $k,"=>",$hash{$k},"\n";
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章