基本信息
源码名称:提取脚本
源码大小:0.83KB
文件格式:.pl
开发语言:Perl
更新时间:2021-07-15
   友情提示:(无需注册或充值,赞助后即可获取资源下载链接)

     嘿,亲!知识可是无价之宝呢,但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下,绝对物超所值哦!如有下载和支付问题,请联系我们QQ(微信同号):813200300

本次赞助数额为: 2 元 
   源码介绍


open (IN,"<$ARGV[0]") || die ("\nError: Couldn't open misa.pl results file (*.misa) !\n\n");


my $filename = $ARGV[0];

$filename =~ s/\.misa//;

open (SRC,"<$filename") || die ("\nError: Couldn't open source file containing original FASTA sequences !\n\n");
open (OUT,">$filename.p3in");


while (<SRC>){
	chomp;
	if (/^>(\S )/){
		$id=$1;
	}else{
		$seq{$id}.=$_;
	}
}


while(<IN>){
	next unless (/(\S )\t(\d )\t\S \t\S \t(\d )\t(\d )\t\d /g);
    my ($id,$ssr_nr,$size,$start) = ($1,$2,$3,$4);

	$count  ;
    print OUT "SEQUENCE_ID=$id"."_$ssr_nr\nSEQUENCE_TEMPLATE=$seq{$id}\n";
    print OUT "PRIMER_PRODUCT_SIZE_RANGE=100-280\n";
    print OUT "SEQUENCE_TARGET=",$start-3,",",$size 6,"\n";
    print OUT "PRIMER_MAX_END_STABILITY=250\n=\n"
}

print "\n$count records created.\n";
close IN;
close SRC;
close OUT;