Chinaunix首页 | 论坛 | 博客
  • 博客访问: 83773
  • 博文数量: 36
  • 博客积分: 0
  • 博客等级: 民兵
  • 技术积分: 0
  • 用 户 组: 普通用户
  • 注册时间: 2016-07-28 10:06
个人简介

安徽屯溪,哈尔滨的雪地,扬州的瘦西湖,想必知道我是谁了吧!!对,小金思密达

文章分类

全部博文(36)

文章存档

2015年(7)

2014年(29)

我的朋友

分类: PERL

2015-11-26 20:41:33

#!/usr/bin/perl
$example = < >A|c21852_g1f1
XXXXXX
>B|10011958m
XXXXXX
>C|311240
XXXXXX
>D|AT1G01880
XXXXXX
............
E
;

open IN,"
open OUT,">result";
while (){chomp;
if(/^>/){ $key =$_;}else{$hash{$key}.=$_}
}
foreach $line(keys %hash){ chomp$line;
    $line=~/>(\w+)/;
    if ($seq{$1}){ $seq{$1}.=$hash{$line}  }else{ $seq{$1}.=$hash{$line}
}
}

foreach (keys %seq){
print OUT">$_\n$seq{$_}\n";
}
阅读(752) | 评论(0) | 转发(0) |
0

上一篇:去除重复行(依据特征保留一行)

下一篇:没有了

给主人留下些什么吧!~~