24小时热门版块排行榜    

CyRhmU.jpeg
查看: 677  |  回复: 1

zhanggxxyy

新虫 (初入文坛)

[求助] 基因组DNA序列按照位置提取

我有1个基因组DNA序列文件(记为文件A, fasta格式), 现在我要从该文件里面大批量提取A文件的不同指定位置的序列,
比如 我要提取以下2个不同位置的序列(记为文件B,fasta格式)如下   
编号:开始位置-结束位置
1:107-1618
2:1657-2760

把这些提取信息(编号:开始位置-结束位置)都放在文件B,然后根据这些提取信息去提取目标文件A里相对应的序列,
提取出来的序列写入一个新文件里面(记为文件C,fasta格式),每一条提取出来的序列的名称要为所使用的相对应提取信息以及序列具体是什么,比如提取信息为1:107-1618,
提取出来的序列的文件C如下
>1:107-1618
ATCGGGGCTTTTTTTTTACTGCTGAGTCTCAGCGTAGTC
请问如何用perl  代码来实现以上过程.下面是我的程序(见文件夹里.pl文件。) 可是实现不了,求指教
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

生信菜鸟

新虫 (小有名气)

2楼2018-07-05 10:44:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhanggxxyy 的主题更新
信息提示
请填处理意见