设为首页收藏本站
开启辅助访问
切换到窄版

 找回密码
 注册

QQ登录

只需一步,快速开始

搜索
查看: 2939|回复: 3

2012元旦悬赏求助,哪位大侠能帮一下!!!

[复制链接]
发表于 2012-1-2 01:27:35 | 显示全部楼层 |阅读模式
5金钱
我想去除一组序列(大约1000条)中指定区段序列,例如
file1中包含原始待处理序列,格式和内容如下:
>AT2G01210.1
MLASLIIFVALLCNVTVISGLN
>AT2G01820.1
MSNSHLGTLCFIISLLGLANFSLS
>AT2G01950.1
MTTSPIRVRIRTRIQISFIFLLTHL
file2中包含指定区段信息,格式和内容如下:
AT2G01210.1 1 12
AT2G01820.1 2 5,10 14,17 20
AT2G01950.1 3 7,12 15
要求得到file3,内容如下:
>AT2G01210.1
CNVTVISGLN
>AT2G01820.1
MLGTLLLFSLS
>AT2G01950.1
MTVRIRISFIFLLTHL
上面是个简单的例子,实际情况就是序列长一些,序列数多一些。
终端命令形式希望是cutseq.pl file1 file2 file3
哪位高手能够帮助小弟一下,写下perl命令,感激不尽,新的一年必将身体健康工作顺利!!!

 楼主| 发表于 2012-1-2 11:19:57 | 显示全部楼层
来个人看看吧nis30
回复

使用道具 举报

发表于 2012-1-5 11:11:43 | 显示全部楼层
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

QQ|申请友链|小黑屋|手机版|Archiver|生物信息学论坛 ( 蜀ICP备09031721号  

GMT+8, 2017-1-21 12:31 , Processed in 0.144939 second(s), 18 queries .

Powered by Discuz! X3

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表