基本信息
源码名称:提取脚本
源码大小:0.83KB
文件格式:.pl
开发语言:Perl
更新时间:2021-07-15
友情提示:(无需注册或充值,赞助后即可获取资源下载链接)
嘿,亲!知识可是无价之宝呢,但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下,绝对物超所值哦!如有下载和支付问题,请联系我们QQ(微信同号):813200300
本次赞助数额为: 2 元×
微信扫码支付:2 元
×
请留下您的邮箱,我们将在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;