使用Perl从UniProt中下载蛋白质序列

以下代码由PHP站长网 52php.cn收集自互联网

现在PHP站长网小编把它分享给大家，仅供参考

use strict;
use warnings;
use LWP::Simple qw(get);

my $seq = downloadUniprot("P31749");
print $seq;

sub downloadUniprot {
	my $ac = shift;
	my $url = "http://www.uniprot.org/uniprot/".$ac.".fasta";
	my $content = get($url);
	die "Couldn't get '$ac'!" unless defined $content;
	my @line = split/R/,$content;
	shift @line;
	my $seq = "";
	map {$seq .= $_} @line;
	return $seq;
}

以上内容由PHP站长网【52php.cn】收集整理供大家参考研究

如果以上内容对您有帮助，欢迎收藏、点赞、推荐、分享。

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!