| 查看: 734 | 回复: 1 | ||
[求助]
基因组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文件。） 可是实现不了，求指教  |
回复此楼» 猜你喜欢
时间戳又变了8-15
已经有20人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
生信菜鸟
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 184.6
- 散金: 10
- 帖子: 73
- 在线: 10.6小时
- 虫号: 6037865
- 注册: 2017-03-18
- 性别: GG
- 专业: 分子与进化生态学
2楼2018-07-05 10:44:34
