| 查看: 480 | 回复: 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 |
回复此楼» 猜你喜欢
宠物寄生虫防治：氟雷拉纳阻断GABA受体，持效12周驱杀跳蚤蜱虫
已经有0人回复
-
为呼吸打开通道：依伐卡托的科学之旅
已经有1人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有138人回复
-
Dordaviprone（ONC201）：小分子药物在中线胶质瘤中的应用
已经有0人回复
-
依喜替康：新型喜树碱衍生物的研究进展
已经有1人回复
-
膀胱癌靶向治疗新选择：厄达替尼作用机制与耐药研究进展
已经有0人回复
-
从水母到实验室：腔肠素的发现历程与生物发光奥秘
已经有1人回复
-
线粒体氧化磷酸化的新靶点：S-Gboxin的发现与研究进展
已经有0人回复
-
PROTAC药物开发选择VHL配体：MDK-7526以87%出口向量占有率成首选
已经有0人回复
-
MCC950：NLRP3炎症小体特异性抑制剂的科研应用与前景
已经有0人回复
-
康替唑胺研发17年：如何解决多重耐药菌感染难题
已经有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 (幼儿园)
- 金币: 106.5
- 帖子: 156
- 在线: 51.3小时
- 虫号: 2045697
- 注册: 2012-10-06
- 性别: GG
- 专业: 肿瘤生物治疗
3楼2014-09-30 14:03:31