| 查看: 521 | 回复: 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 |
回复此楼» 猜你喜欢
科研人必读：新一代 KRAS G12C 抑制剂 Olomorasib 的研究进展
已经有0人回复
-
BT-11（Omilancor）研究化合物全解析：从分子结构到实验应用要点
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有248人回复
-
富马酸贝达喹啉：二芳基喹啉类工具化合物的化学特征
已经有0人回复
-
BT投稿with editor状态快两周了，正常吗
已经有25人回复
-
Silevertinib（BDTX-1535）结构拆解：炔基双环片段为何成为合成瓶颈
已经有0人回复
-
RO4929097 解析：双苯并氮杂䓬骨架与 Notch 通路实验设计要点
已经有0人回复
-
分子量 1157 的 PROTAC 怎么配不析出：SGK3-PROTAC1 的溶解、分装与图谱复核
已经有0人回复
-
logP 6.1 的羧酸怎么配母液：GSK 650394 的溶解窗口与稀释顺序
已经有0人回复
-
从"水不溶"到 71 mg/mL：雷特格韦的 pH 依赖溶解与螯合色谱边界
已经有0人回复
-
133 Ų 的极性表面积，如何决定 TNG462 的母液配制与色谱条件
已经有0人回复
wizardfan
至尊木虫 (著名写手)
- MolEPI: 4
- 应助: 599 (博士)
- 贵宾: 1.818
- 金币: 24660.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 (幼儿园)
- 金币: 147
- 帖子: 157
- 在线: 51.7小时
- 虫号: 2045697
- 注册: 2012-10-06
- 性别: GG
- 专业: 肿瘤生物治疗
3楼2014-09-30 14:03:31
