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

wangyufeng的博客

祝愿BB 健康开心快乐每一天

 
 
 

日志

 
 

phylip2fasta:Converts an aligned fasta seq file to phylip format  

2012-03-16 16:02:53|  分类: Perl & bioperl |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
#!/usr/bin/perl  my $usage = "Usage: $0 [-h] [-v] [infile]\n" .             "  -h  help\n" .             "  -v  verbose (print name conversion in STDERR)\n" .             " infile should be a phylip or paml format (one liner), " .             "STDIN is used if no infile is given\n";  use IO::File; use Getopt::Std; getopts('hv') || die "$usage\n";  die "$usage\n" if (defined ($opt_h));  my $totNumChar = 10;  # number of characters allowed for name in phylip my $numFrontChar = 7; # When the name is too long, this amount of characters                       # are used from the beginning of the name, and the rest                       # are from the end of the name.  while(<>){      next if ($. == 1);  # skip the first line      chomp;     s/^\s+//; s/\s$//;     next if (/^$/);      my @line = split (/\s+/);      my @nameChar = split (//, $line[0]);      if (@nameChar > $totNumChar) {  if ( /^(.{$totNumChar})/ ) {      $name = $1;  }     } else {  $name = $line[0];     }      s/$name//;     s/^\s+//;      print ">$name\n$_\n"; }
  评论这张
 
阅读(1139)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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