本文提供了一个使用Perl语言处理FASTA格式生物序列数据的简单脚本示例。该脚本展示了如何读取FASTA文件，并获取序列ID、名称、序列本身及其长度等基本信息。
#! /usr/bin/perl
use Bio::SeqIO;
my $in=Bio::SeqIO->new(-file=>'/home/XXX/XXX/all.pep.fa',-format=>'fasta');
my $i=0;
while((my $s=$in->next_seq)&&($i <=10))
{ my $id=$s->id;
my $seq=$s->seq;
$i++;
my $name=$s->display_name;
my $len=$s->length;
print "$id\n$seq\n$name\n$len\n";
}
use Bio::SeqIO;
my $in=Bio::SeqIO->new(-file=>'/home/XXX/XXX/all.pep.fa',-format=>'fasta');
my $i=0;
while((my $s=$in->next_seq)&&($i <=10))
{ my $id=$s->id;
my $seq=$s->seq;
$i++;
my $name=$s->display_name;
my $len=$s->length;
print "$id\n$seq\n$name\n$len\n";
}
扫一扫
1308
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
