免费注册 查看新帖 |

Chinaunix

  平台 论坛 博客 文库
最近访问板块 发新帖
查看: 1825 | 回复: 5
打印 上一主题 下一主题

求程序_生物信息_PERL [复制链接]

论坛徽章:
0
跳转到指定楼层
1 [收藏(0)] [报告]
发表于 2013-01-25 22:32 |只看该作者 |倒序浏览
各位高手,小弟最近在做系统进化分析,需要用PERL来处理一些文件,基础薄弱,还请各位高手帮小弟编写一个小程序!!
程序的目的:将interested_domain.txt文件(小弟已经贴上了)中,每个基因所对应的序列长度小于等于80的,将这样的序列取出来(还需要基因名字,文件格式还是fasta格式,即>后面是基因的名字,下面是基因的序列)。
序列文件片段如下:

红色框是每个基因对应的序列(有长有短),蓝色框里的数字,后面的较大数字减去前面的较小数字再加1就是序列的长度。
文件:
interested_domain.rar (6.84 KB, 下载次数: 27)
麻烦各位大神了!

论坛徽章:
0
2 [报告]
发表于 2013-01-25 23:43 |只看该作者
为啥有两个数字?

论坛徽章:
0
3 [报告]
发表于 2013-01-25 23:47 |只看该作者
回复 2# zhaoshuaijiang
这些序列是我从另外一个文件中截取的,前面那个数字是结构域的起始位点,后面的数字是结构域的结束位点。

论坛徽章:
0
4 [报告]
发表于 2013-01-25 23:51 |只看该作者
就是用这两个数字截取的
jingwu110 发表于 2013-01-25 23:47
回复 2# zhaoshuaijiang
这些序列是我从另外一个文件中截取的,前面那个数字是结构域的起始位点,后面的数 ...

论坛徽章:
0
5 [报告]
发表于 2013-01-26 23:15 |只看该作者
需要找出的是 这两个数字的差值小于等于80?

论坛徽章:
0
6 [报告]
发表于 2013-01-27 13:57 |只看该作者
经过请教 问题已经解决了,谢谢各位!
您需要登录后才可以回帖 登录 | 注册

本版积分规则 发表回复

  

北京盛拓优讯信息技术有限公司. 版权所有 京ICP备16024965号-6 北京市公安局海淀分局网监中心备案编号:11010802020122 niuxiaotong@pcpop.com 17352615567
未成年举报专区
中国互联网协会会员  联系我们:huangweiwei@itpub.net
感谢所有关心和支持过ChinaUnix的朋友们 转载本站内容请注明原作者名及出处

清除 Cookies - ChinaUnix - Archiver - WAP - TOP