24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1569  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

湖人湖人

铁杆木虫 (小有名气)

[求助] 请帮忙编写一个perl程序处理文本(提取GI号) 已有1人参与

现在我碰到这样一个问题,我有一个文本文件1.txt,是Blastx的结果,内容如下:
gi|293567070|gb|FS329526.1|FS329526,gi|219980050|emb|CAX06392.1|,33.96,106,70,0,406,723,256,361,5e-014,72.0
gi|293567070|gb|FS329526.1|FS329526,gi|219980050|emb|CAX06392.1|,30.43,161,107,3,244,711,28,188,7e-013,68.2
gi|293567070|gb|FS329526.1|FS329526,gi|219980050|emb|CAX06392.1|,32.31,130,84,2,334,711,180,309,1e-012,67.4
gi|293567065|gb|FS329521.1|FS329521,gi|388506924|gb|AFK41528.1|,42.27,220,115,1,107,766,67,274,9e-043, 167
gi|293567064|gb|FS329520.1|FS329520,gi|388510998|gb|AFK43565.1|,80.84,167,32,0,239,739,16,182,1e-079, 289
。。。。。。
后面还有很多序列,总共大概3M左右。
现在想提取前面核酸序列的GI号:gi|293567070,gi|293567065。。。输出为txt文本,其中重复的只输出一个,
另外,再输出一个文本,提取后面蛋白序列的GI号:gi|219980050,gi|388506924。。。重复的也是只输出一个。
各位大虾,请问怎样编写一个perl程序处理上面的问题啊?谢谢了!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

christina1213

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
湖人湖人: 金币+10, 有帮助, 用Excel删除重复项 2012-06-25 21:57:16
xzhdty: 金币+1, 谢谢 2012-06-26 06:19:02
湖人湖人: 金币+10, ★★★很有帮助 2012-06-26 21:02:39
#!/usr/bin/perl
use strict;
use warnings;

open (IN, "<1.txt";
open (OUT, ">>out.txt";

while (my $line= {
        my @cut=split("|gb|",$line);
        print OUT $cut[0]."\n";
}
close IN;
close OUT;
去重的话用EXCEL吧,一下就搞定了,不用写程序。
3楼2012-06-25 20:04:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

longwen36

铁杆木虫 (正式写手)

用正则表达式吧,规律一样的话。用matlab我会。
2楼2012-06-25 12:49:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

感谢参与,应助指数 +1
python的解法:正则+map数据结构应该是最简单的.
不知道perl有没有map类似的数据结构.
正则提取数据,map去重
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
4楼2012-06-26 00:04:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mycaas

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
cat 1.txt |perl -ne 'chomp;while(/(gi\|\d+)/g){print "$1\n";}' |sort |uniq
5楼2012-06-27 00:59:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601材料工程找调剂 +13 oatmealR 2026-03-29 14/700 2026-03-29 16:41 by 学员8dgXkO
[考研] 一志愿211,335分,0856,求调剂院校和导师 +5 倾____萧 2026-03-27 6/300 2026-03-29 16:35 by 唐沐儿
[考研] 085600 286分 材料求调剂 +8 麻辣鱿鱼 2026-03-27 9/450 2026-03-29 15:36 by Serene1974
[考研] 一志愿郑州大学,080500学硕,总分317分求调剂 +8 举个栗子oi 2026-03-24 9/450 2026-03-29 13:08 by peike
[考研] 317分 一志愿南理工材料工程 本科湖工大 求调剂 +9 芋泥小铃铛 2026-03-28 9/450 2026-03-29 11:46 by longlotian
[考研] 348求调剂 +5 小懒虫不懒了 2026-03-28 5/250 2026-03-29 10:34 by 唐沐儿
[考研] 305求调剂 +8 RuiFairyrui 2026-03-28 8/400 2026-03-29 08:22 by fmesaito
[考研] 求调剂 +3 QiMing7 2026-03-25 4/200 2026-03-28 14:30 by QiMing7
[考研] 286求调剂 +12 PolarBear11 2026-03-26 12/600 2026-03-28 12:14 by zllcz
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +9 丹青奶盖 2026-03-26 10/500 2026-03-28 07:45 by barnett0632
[考研] 求调剂推荐 材料 304 +15 荷包蛋hyj 2026-03-26 15/750 2026-03-28 04:13 by fmesaito
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[考研] 266求调剂 +11 阳阳哇塞 2026-03-27 12/600 2026-03-27 17:56 by yu221
[考研] 304材料求调剂 +4 钟llll 2026-03-26 4/200 2026-03-27 03:42 by wxiongid
[考研] 325求调剂 +5 李嘉图·S·路 2026-03-23 5/250 2026-03-27 00:42 by wxiongid
[考研] 333求调剂 +6 wfh030413@ 2026-03-23 6/300 2026-03-26 22:45 by 学员8dgXkO
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
[考研] 一志愿吉林大学材料与化工303分求调剂 +4 为学666 2026-03-24 4/200 2026-03-25 11:27 by BruceLiu320
[考研] 086003食品工程求调剂 +6 淼淼111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
信息提示
请填处理意见