永发信息网

大神们,本人新手,求一个perl程序,谢谢

答案:2  悬赏:50  手机版
解决时间 2021-11-25 11:00
大神们,本人新手,求一个perl程序,谢谢
最佳答案
#!/usr/bin/perl -w
use strict;
die "perl $0 \n" unless(@ARGV == 1);
open IN,$ARGV[0];
$/=">";
;
$/="\n";
while(){
chomp $_;
my $seqName = $_;
$/=">";
my $seq = ;
chomp $seq;
$/ = "\n";
if($seqName =~/ENSG00000179088/){
$seq =~s/[\r\n]//g;
my $result = substr($seq,37,52);
print "提取的子序列为:\n$result\n";
}
}
close IN;

使用的时候,将假设将脚本存为 picseq.pl
然后序列所在fa文件为 test.fa
运行的时候再命令行输入 picseq.pl test.fa
回车即可。
如果有疑问,可以随时追问,希望能够帮到你
全部回答
$/=">";
<>;
while(<>){

chomp;
@tmp=split/\n/,$_,2;
if ($tmp[0] eq ENSG00000179088){
$tmp[1]=~s/\s+//g;
print ">$tmp[0]\n$tmp[1]\n";
$seq=substr($tmp[1],37,52);
print "$seq\n";
}
}
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
新疆天富热电股份公司怎么样呀?据说是上市企
太平鸟男装吊牌B2BB42408是是什么意思
当归长期吃有没有副作用
为什么我睡觉做梦,梦中很真实,可是醒来才知
七夕我要表白我的大学女同学,是说我爱你还是
圆的集合定义是什么?求解答
我想学黄辣丁鱼苗繁殖技术
南方航空,国内航班,能随身带液体化妆品上机
常州方程物流徐州专线
“林”在古代是什么意思?
初戴眼镜看东西是弯曲的吗
古时雍城是什么地方
新版新浪微博的"我的应用"在哪儿
我想买一辆江淮货车,不知道该买骏铃还是康铃
even, much as, even still都是尽管的意思?
推荐资讯
猎狗鼻子在行动读后感400字
糖斗广场舞茉莉老师简历
急!睡过头了,今天都不想上班。要怎么跟老板
公交北京市的运城街至劲松中路
汗宫养生地址在哪,我要去那里办事
为什么说上帝很伟大
词语欣漠安然什么意思?
一年级20个球每组4个可以分几组
宏兴灯饰五金超市地址在哪,我要去那里办事,
甘婷婷新水浒传床戏在第几集
微五币对换劲霸夹克是什么时间
天天乐棋牌地址在什么地方,我要处理点事!
正方形一边上任一点到这个正方形两条对角线的
阴历怎么看 ?