makePindelList.pl

#!/usr/bin/perl
#

&help if (@ARGV != 2);

open(IN, "<$ARGV[1]") or die;
while ($line = <IN>) {
  chomp($line);
  next if ($line =~ /^Locus/);
  @element = split(/\t/,$line);
  @ele = split(/:/,$element[0]);
  $depth{$ele[0]}[$ele[1]] = $element[1];
}
close IN;

open(IN, "<$ARGV[0]") or die;
while ($line = <IN>) {
  chomp($line);
  next if ($line !~ /BP_range/);
  @element = split(/\s+/,$line);

  $sv_type = $element[1];
  $sv_size = $element[2];
  $chr1 = $element[7];
  $chr2 = $element[7];
  $pos1 = $element[9];
  $pos2 = $element[10];

  next if ($chr1 eq "Mt");
  next if ($chr1 eq "Pt");
  next if ($sv_size < 10);

  $freq{$sv_type} ++;

  $depth_val = ($depth{$chr1}[$pos1] + $depth{$chr2}[$pos2]) / 2;

  $key = "$chr1\t$pos1\t$chr2\t$pos2\t$sv_type";
  $inf{$key} = "$sv_size\t$element[26]\t$element[15]\t$depth_val";

  $dat{chr} = $chr1;
  $dat{pos} = $pos1;
  $dat{key} = $key;
  push(@keysort,{%dat});
}
close IN;

@keysort = sort{$a->{chr}<=>$b->{chr} || $a->{pos}<=>$b->{pos}} @keysort;

for($i=0; $i<@keysort; $i++) {
  $key = $keysort[$i]{key};
  print STDOUT "$key\t$inf{$key}\n";
}

foreach $key (sort keys %freq) {
  print STDERR "$key\t$freq{$key}\n";
}

exit(0);


### SUB : help
sub help {
  print STDERR "USAGE: makePindelList.pl pindel.output depth\n";
  exit(-1);
}