提供了参考序列,需要根据起始位置和终止位置提取参考序列上固定的碱基序列,碱基位置及参考序列如下表

参考基因组 起始位置 终止位置
NC_002645.1 293 12550
NC_006577.2 206 21753
NC_006213.1 210 21496
NC_005831.2 287 20475
NC_004718.3 265 21485
NC_019843.3 279 21514
NC_045512.2 266 21555

参考序列文件如下:

参考基因组 参考序列
NC_002645.1 ATCG…
NC_006577.2 ATCG…
NC_004718.3 ATCG…


提取脚本如下:

use strict;
use warnings;my $in=shift;
open IN,"<$in" or die "pleaes give a file.\n";
print "参考基因组\t序列\n";
while(<IN>)
{chomp;my @a=split("\t",$_);open IN2,"<$a[0]";while(<IN2>){chomp;my $length=$a[2]-$a[1]+1;my $seq=substr($_,$a[1],$length);print "$a[0]\t$seq\n";}  close IN2;
}
close IN;

perl函数substr应用举例相关推荐

  1. perl Carp模块使用举例(转)

    perl Carp模块使用举例 perl编程 2009-08-20 19:16:42 阅读262 评论0   字号:大中小 订阅 近日分析extmail的源代码,看到Carp的perl模块,有些搞不懂 ...

  2. php error_reporting(0);,php error_reporting()函数的用法举例(错误捕捉)

    php error_reporting()函数的用法举例,供大家学习参考. 学习CI框架过程中遇到个问题: A PHP Error was encountered Severity: Notice M ...

  3. R语言missing()函数用法及举例

    R语言missing()函数用法及举例 函数作用就是如果y没有具体的值,那么输出为真,如果y有具体的值,那么输出为假,这么说可能有点抽象,举个例子: 对于这样一段代码 myplot <- fun ...

  4. Perl 字符串截取函数substr

    在编写perl程序的时候,我们有时需要截取一个字符串中的部分内容. 这个时候通常会使用substr函数实现这个功能. 函数参数: substr(待截取字符串,截取起始位置,截取终了位置); view ...

  5. js字符串截取函数substr substring slice使用对比

    2019独角兽企业重金招聘Python工程师标准>>> JS中常用三个的字符串截取函数:substr substring slice,调用方式如下: stringObject.sli ...

  6. 关于string的成员函数substr

    在LKSNAPSHOT.cpp中解析接收串口的字符串时有一句: string cmd = msg.substr(cmdstr.length()+1,cmdstr.length()-cmdstr.len ...

  7. Perl函数pack/unpack(二进制读写)

    函数名 pack 调用语法 formatstr = pack(packformat, list); 解说 把一个列表或数组以在实际机器存贮格式或C等编程语言使用的格式转化(包装)到一个简单变量中.参数 ...

  8. C++ 字符串复制函数substr

    关于substr函数,回过头来看自己最初写的技术记录才发现有一个错误 正确的格式是substr((int型)下标X,(int型)长度Len); 解释为截取从下标X开始的Len个字符串到一个string ...

  9. Perl函数及参数介绍

    from: http://www.cnblogs.com/polestar/archive/2012/05/03/2480118.html 子例程,函数 子例程,函数的作用域是指在程序中能够看到它的位 ...

最新文章

  1. 2014年7月17日学习笔记--PHP的循环结构学习
  2. FaaS — Serverless Computing(无服务器计算)
  3. python selenium定位元素方法,python + selenium 练习篇 - 定位元素的方法
  4. 第三次学JAVA再学不好就吃翔(part60)--Arrays类
  5. jquery checkbox 选中、改变状态、change 和 click 事件
  6. PCIe的事务传输层的处理(TLP)
  7. 躺枪实惨!与新冠病毒“重名”,这款墨西哥国民啤酒宣布暂停生产...
  8. 【Java】数据库编程
  9. pyspark指定schema
  10. 毕设 数据库ER图绘制过程中的参考资料
  11. Crucible安装使用
  12. 三重积分--------球坐标系
  13. 那些年逃离北京的软件测试员,现在都怎么样了?
  14. Tiled结合Unity实现瓦片地图
  15. 怎么压缩图片 ? 掌握这几种免费压缩图片的方法就够了
  16. 四种内容付费盈利模式介绍
  17. echarts实现甘特图
  18. SCI与EI检索是什么意思,包括哪些专业?
  19. 信息无障碍研究机构---教育、公益等
  20. Visual Studio 2022安装与使用教程

热门文章

  1. altium 原理图reset unique id
  2. mysql连接我的sql主机_linux,centos安装mysql,windows主机远程连接mysql教程
  3. 虚拟机字节码执行引擎-方法调用
  4. [BZOJ3032]七夕祭
  5. MMDetection(五)
  6. 【运行报错】Redis:Operation against a key holding the wrong kind of value
  7. 【加速乐SCDN】正式上线,“安全+加速”一步到位!
  8. javascript练习8:判断用户名是否合法,正确应由英文字母或数字组成,且长度是6-12位
  9. php qqwry.dat_php读取qqwry.dat ip地址数据库文件程序
  10. 百度地图——点标记及点击弹窗显示