#!/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";
}
統計fasta文件每條序列的長度【Perl】
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章
[轉載] Perl的內置變量(2)
$- 當前頁可打印的行數,屬於Perl格式系統的一部分$! 根據上下文內容返回錯誤號或者錯誤串$” 列表分隔符$# 打印數字時默認的數字輸出格式$$ Perl解釋器的進