提取psite位点附近序列

清华大佬耗费三个月吐血整理的几百G的资源,免费分享!....>>>

use strict;
use warnings;

my $seq = "MTVLLLVVLQMWKATAGHSIAVSQDDGADDWETDPDFVNDVSEKEQRWGAKTVKGSGHQE";

my $peptide = getPeptide($seq,2,7,7);
print $peptide;

sub getPeptide {
	my ($seq, $site, $left, $right) = @_;
	my $len = length($seq);
	die " > seq length: $len\n" if $len < $site;
	die "'$site' < 1\n" if $site < 1;	
	$seq = ("*" x $left).$seq.("*" x $right);
	my $peptide = substr($seq, $site - 1, $right + $left + 1);
	return $peptide;
}