24小时热门版块排行榜

>论坛更新日志 (3024)
>虫友互识 (324)
>文献求助 (273)
>导师招生 (173)
>论文投稿 (106)
>休闲灌水 (100)
>考博 (87)
>仿真模拟 (63)
>基金申请 (52)
>硕博家园 (43)
>博后之家 (39)
>找工作 (31)
>招聘信息布告栏 (30)
>公派出国 (27)
>考研 (26)
>SciFinder/Reaxys (25)

返回列表

zhanggxxyy

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 57
帖子: 14
在线: 9.3小时
虫号: 5805175
注册: 2017-03-02
专业: 动物遗传学

[求助] 基因组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文件。）可是实现不了，求指教