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

Puriney's Notes

Puriney=purine+Y, my Wonderland

 
 
 

日志

 
 

[bio-tips] genomeCoverageBed: BAM2Wig  

2012-09-26 10:51:42|  分类: Bio |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

Last post  introduced the basic usage of genomeCoverageBed. Here is some update to deal with BAM2WIG processing. 

As the default output of genomeCoverageBed is formatted as following: 
chr start depth

However, WIG file format is different. The chromosome information is recorded in the header line:

variableStep chrom=chr1 span=1
Thus needed information exact work should be done. 
Here I share with you my homemade Perl script that manages to help you out. 
#!/usr/bin/perl -w
my ($in, $out) = @ARGV;
die $! && <DATA> unless (@ARGV ==2);
 
open IN, $in;
open OUT, ">$out";
 
print OUT "track type=wiggle_0 name="NAME" color=250,100,100 visibility=full autoScale=off viewLimits=0.0:99.0 yLineOnOff=off priority=10\n";
my $tmp = "NA";
while ( <IN> ) {
    chomp;
    my ($chr, $start, $v) = (split /\t/, $_)[0,1,2];
    my $tag = 0;
    $tag = 1 if ($tmp eq $chr);
         if ($tag == 1 ){
             print OUT "$start\t$v\n";
         }
         elsif ($tag == 0){
              $tmp = $chr;
              print OUT "variableStep chrom=$chr span=1\n";
              print OUT "$start\t$v\n";
         }
}   
== End && reference
 
== afters
thx god and thx girlfriend
  评论这张
 
阅读(788)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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