1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
#!/usr/bin/perl
#-------------------------------------------------------------------------------
#
# Split a WIG file by chromosome
#
#
# Pablo Cingolani
#-------------------------------------------------------------------------------
$chrPrev = "";
# Parse STDIN
while( $l = <STDIN> ) {
$chr = "";
if( $l =~ /^variableStep\s+chrom=(\S+)/ ) {
$chr = $1;
# Open new chromosome file?
if( $chr ne $chrPrev ) {
# Close old file
close OUT if $chrPrev ne "";
# Open new file
$file = "$chr.wig";
print STDERR "Creating new file '$file'\n";
open OUT, ">$file";
$chrPrev = $chr;
}
}
print OUT $l;
}
close OUT;
|
