注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Bioinformatics home

 
 
 

日志

 
 

Bioperl:把Genbank格式的序列转换为基因结构图  

2009-09-24 22:56:15|  分类: 生物信息编程 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

下面是一个BioPerl的小程序,可以把Genbank格式的序列转换为png图片-基因结构图,显示序列的长度,CDS区,exon区,STS区等。简单地讲就是把该序列的Genbank格式里的信息用图片表示。 这个程序来源于网络。

下面举个例子,以NM_172587为例。

NM_172587的Genbank格式下面的链接,另存为文件名’NM_172587.gb’:

http://www.liucheng.name/entrez/entrez.php?acc=NM_172587&db=nucleotide&format=gb&start=Submit

下面是Bioperl脚本:

#!/usr/bin/Perl -w
use Bio::Graphics;

use Bio::SeqIO;

use Bio::SeqFeature::Generic;

my $file = shift    or die "provide a sequence file as the argument";

my $io = Bio::SeqIO->new(-file=>$file)    or die "couldn't create Bio::SeqIO";

my $seq = $io->next_seq    or die "couldn't find a sequence in the file";

my @features = $seq->all_SeqFeatures;# sort features by their primary tags
my %sorted_features;

for my $f (@features) {
  my $tag = $f->primary_tag;
  push @{$sorted_features{$tag}},$f;
  }
 
  my $panel = Bio::Graphics::Panel->new(
    -length => $seq->length,
    -key_style => 'between',
    -width => 800,
    -pad_left => 10,
    -pad_right => 10
   );
 
  $panel->add_track(
          arrow =>Bio::SeqFeature::Generic->new(-start => 1,-end => $seq->length),
    -bump => 0,
    -double=>1,
    -tick => 2
    );
 
  $panel->add_track(generic =>Bio::SeqFeature::Generic->new(-start => 1,-end => $seq->length,-bgcolor => 'blue', -label => 1));# general case
  my @colors = qw(cyan orange blue purple green  chartreuse magenta yellow aqua);
  my $idx = 0;   
 
  for my $tag (sort keys %sorted_features)
  {
          my $features = $sorted_features{$tag};
    $panel->add_track(
       $features,-glyph => 'generic',
       -bgcolor => $colors[$idx++ % @colors],
       -fgcolor => 'black',
       -font2color => 'red',
       -key => "${tag}s",
       -bump => +1,
       -height => 8,
       -label => 1,
       -description => 1
     );   
  }
 
  print $panel->png;

存为文件名:graphics.pl

用法:

perl graphics.pl NM_172587.gb >1.png

结果就生成了一个名为1.png的图片,如下图

Bioperl:把Genbank格式的序列转换为基因结构图 - xiaofeng1982 - Leon

当然了,用这个脚本首先你的电脑已经装好了Bioperl。这个就不多讲了~

  评论这张
 
阅读(1515)| 评论(0)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017