| 查看: 456 | 回复: 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 |
回复此楼» 猜你喜欢
留学--博士招生
已经有16人回复
-
化学工程，本211，求调剂，ab区不限
已经有4人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有51人回复
-
湖北师范大学招调剂啦
已经有7人回复
-
邵阳学院食品与化学工程学院硕士调剂
已经有0人回复
-
天津商业大学 生物与医药课题组招生
已经有0人回复
-
生物化工学硕招收调剂
已经有0人回复
-
广东石油化工学院2026年硕士研究生需要多名生物，制药工程，食品专业的调剂生
已经有0人回复
-
317分 一志愿江南大学 化学工程学硕 求调剂
已经有6人回复
-
化工申博
已经有1人回复
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 (幼儿园)
- 金币: 82
- 帖子: 156
- 在线: 51.3小时
- 虫号: 2045697
- 注册: 2012-10-06
- 性别: GG
- 专业: 肿瘤生物治疗
3楼2014-09-30 14:03:31
