| 查看: 431 | 回复: 2 | ||
wang4681060铁虫 (小有名气)
|
[求助]
perl 已有1人参与
|
|
我做了个blast，m8格式，我想从blast结果里提取query-start到query-end这一段query序列，求大神帮写个perl代码，不甚感激！！！ Contig203 gi|157325286|ref|YP_001468706.1|_TerL_Listeria_phage_B054 27.83 503 341 14 1255 2697 5 469 9e-58 201 就是从Contig203中提取1255到2697这一段序列，输出 >Contig203-gene_1 序列 第二个提取的序列>ContigXXX-gene_2 |
回复此楼» 猜你喜欢
筑牢营养安全线：以精准检测，护健康基石
已经有0人回复
-
推荐一些20种氨基酸检测的实际应用案例
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有133人回复
-
不合理蛙科研实验中的趣事：实验器材的 “乌龙”
已经有0人回复
-
不合理蛙科研实验中的趣事：和实验材料的 “斗智斗勇”
已经有0人回复
-
蛋白质检测：精准分析，解锁生物分子的密码
已经有0人回复
-
不合理蛙科研实验之小鼠实验：严谨设计，解析生命机制的重要载体
已经有0人回复
-
不合理蛙科研实验之重金属检测：精准筛查，守护健康与环境的防线
已经有0人回复
-
不合理蛙科研实验之“蛙测重金属我背锅三千”
已经有0人回复
-
不合理蛙科研实验之“鼠逃三次我发三篇SCI”
已经有0人回复
wizardfan
至尊木虫 (著名写手)
- MolEPI: 4
- 应助: 599 (博士)
- 贵宾: 1.818
- 金币: 24632.2
- 散金: 197
- 红花: 48
- 沙发: 2
- 帖子: 2254
- 在线: 400.7小时
- 虫号: 1879241
- 注册: 2012-07-05
- 性别: GG
- 专业: 生物信息学
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
kx444555: 金币+3, 鼓励交流 2014-09-30 12:36:29
感谢参与，应助指数 +1
kx444555: 金币+3, 鼓励交流 2014-09-30 12:36:29
|
基本思路 1. 读取fasta文件，存入一个hash，key为contig***, value为具体序列 2. 读取blast 结果文件 每行 my ($id, undef, undef, undef, undef, undef, undef, $start, $end) = split("\t", $line) 3. 获得序列 $seq = $hash{$id}; my $result = substr ($seq, $start, $end-$start+1); |
2楼2014-09-30 06:25:28
wang4681060
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 54.5
- 帖子: 154
- 在线: 50.2小时
- 虫号: 2045697
- 注册: 2012-10-06
- 性别: GG
- 专业: 肿瘤生物治疗
3楼2014-09-30 14:03:31