| 5 | 1/1 | 返回列表 |
| 查看: 1580 | 回復(fù): 8 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
湖人湖人鐵桿木蟲 (小有名氣)
|
[求助]
請(qǐng)幫忙編寫一個(gè)perl程序處理文本(提取GI號(hào)) 已有1人參與
|
|
|
現(xiàn)在我碰到這樣一個(gè)問題,我有一個(gè)文本文件1.txt,是Blastx的結(jié)果,內(nèi)容如下: 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左右。 現(xiàn)在想提取前面核酸序列的GI號(hào):gi|293567070,gi|293567065。。。輸出為txt文本,其中重復(fù)的只輸出一個(gè), 另外,再輸出一個(gè)文本,提取后面蛋白序列的GI號(hào):gi|219980050,gi|388506924。。。重復(fù)的也是只輸出一個(gè)。 各位大蝦,請(qǐng)問怎樣編寫一個(gè)perl程序處理上面的問題。恐x謝了! |
至尊木蟲 (著名寫手)
驃騎將軍

鐵桿木蟲 (正式寫手)
木蟲 (正式寫手)
|
#!/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吧,一下就搞定了,不用寫程序。 |
鐵桿木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料求調(diào)劑 +9 | 呢呢妮妮 2026-04-01 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 08生物與醫(yī)藥專碩初試346找調(diào)劑 +5 | dianeeee 2026-04-01 | 6/300 |
|
|
[考研] 366求調(diào)劑一志愿東北大學(xué) +3 | 運(yùn)氣來得若有似?/a> 2026-04-02 | 3/150 |
|
|
[考研] 304求調(diào)劑 +12 | 素年祭語 2026-03-31 | 15/750 |
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-04-01 | 3/150 |
|
|
[考研] 材料專碩322分 +10 | 哈哈哈吼吼吼哈 2026-04-01 | 10/500 |
|
|
[考研] 273求調(diào)劑 +19 | 李芷新1 2026-03-31 | 19/950 |
|
|
[考研] 292求調(diào)劑 +17 | 木蟲er12138 2026-04-01 | 17/850 |
|
|
[考研] 307分求調(diào)劑 +14 | (o~o) 2026-03-31 | 15/750 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 5/250 |
|
|
[考研] 一志愿中國(guó)科學(xué)院大學(xué)265求調(diào)劑 +8 | 恬淡ye 2026-03-31 | 9/450 |
|
|
[考研] 291求調(diào)劑 +3 | 迷蒙木木 2026-04-01 | 4/200 |
|
|
[考研] 求0861交通運(yùn)輸專碩or材料專碩調(diào)劑 +4 | 勒布朗@ 2026-03-31 | 4/200 |
|
|
[考研] 考研材料工程351分調(diào)劑 +5 | 整個(gè)好的 2026-03-31 | 5/250 |
|
|
[考研] 340求調(diào)劑 +4 | 希望如此i 2026-03-31 | 4/200 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +13 | 辰沐5211314 2026-03-26 | 13/650 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|