24小时热门版块排行榜    

Znn3bq.jpeg
查看: 795  |  回复: 2

dadalia

木虫 (小有名气)

[求助] perl 语言小程序求助

希望从输入文件中提取一部分信息作为输出文件,需要用perl语言写,麻烦大家给个最标准的写法,这是一个小测试,不知道如何写才是最能得分的
输入文件:

   SW   perc perc perc  query                position in query        matching       repeat           position in repeat
score   div. del. ins.  sequence             begin  end      (left)   repeat         class/family   begin  end    (left)  ID

  476   22.2 11.7  4.1  AC189398-KBrB056G23    3194   3423 (103578) + ATHATN3        DNA/hAT-Ac        269    554    (0)   2  
  551    1.5  0.0  1.5  AC189398-KBrB056G23    4834   4902 (102099) + (TA)n          Simple_repeat       1     68    (0)   3

要求的输出文件:
##gff-version 3

AC189398-KBrB056G23        RepeatMasker        dispersed_repeat        3194        3423        476        +        .        Target=ATHATN3 269 554
AC189398-KBrB056G23        RepeatMasker        dispersed_repeat        4834        4902        551        +        .        Target=(TA)n 1 68

看起来可能比较乱,附件中提供了文本格式的,谢谢
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 2.rar
  • 2013-04-13 10:41:30, 738 bytes
  • 附件 2 : 2.rar
  • 2013-04-13 10:38:38, 738 bytes

» 收录本帖的淘帖专辑推荐

perl语言专栏

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

possetcn

铁虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
dadalia: 金币+6, ★★★★★最佳答案 2013-04-15 23:46:27
将Togff3.pl放在与RepeatMasker.out同一个目录下,运行 perl Togff3.pl  即可:
CODE:
#!/usr/bin/perl -w
use strict;

#the input filename
my $in = "RepeatMasker.out";
open IN, $in or die "Cannot open the $in:$!\n";
open OUT, ">JS.gff3" or die "Cannot create the JS.gff3:$!\n";
print OUT "##gff-versioin 3\n\n";

while (){
  if (/^\s.+\d/){
      my @LN = split /\s+/;
      print OUT "$LN[5]\tRepeatMasker\tdispersed_repeat\t$LN[6]\t$LN[7]\t$LN[1]\t$LN[9]\t\.\tTarget=$LN[10] $LN[12]\t$LN[13]\n";
  }
}

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Togff3.pl
  • 2013-04-15 17:06:02, 421 bytes
2楼2013-04-15 17:05:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青春不悔

新虫 (小有名气)

楼主,我是perl语言初学者。请问你对gff3格式的perl编程应该很了解了吧?我现在正有这方面的问题 也是把另一种格式转换成gff3格式。。 表示很纠结。
3楼2013-05-26 01:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dadalia 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 河北省自然科学基金 +5 Peterchao 2026-05-18 8/400 2026-05-24 11:58 by 晓晓爱翠翠
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +4 光伏爱好者 2026-05-17 7/350 2026-05-24 10:20 by 光伏爱好者
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 论文撤稿了 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:06 by Equinoxhua
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 pmo95bazuy 2026-05-23 7/350 2026-05-24 06:35 by fpo5ljpv91
[基金申请] 揭秘青基评审内幕:几个A才能顺利中标 +3 国自然国社科中 2026-05-23 4/200 2026-05-23 15:37 by 2000zf36392
[基金申请] 青B发送上会通知了吗 +5 chemBioBro 2026-05-22 7/350 2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +14 1234567wang 2026-05-17 16/800 2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[基金申请] 评审有感 +15 popular289 2026-05-18 26/1300 2026-05-21 10:35 by 西葫芦炒鸡蛋
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
请填处理意见