#!/usr/bin/perl

@gts = ("B6", "C3", "B6C3F1", "C3B6F1");
@genders = ("M", "F");
@reps = (1,2,3);
%dates = (1 => 20130410, 2 => 20130917, 3 => 20130917);

$bigwigDir = "/ddn/gs1/project/mousemeth/methylseq/analysis/sara/analysis_pilot_current/files_coverage";

foreach $gt (@gts) { foreach $mf (@genders) { foreach $r (@reps) {
  $animal = "$gt"."_$mf"."_$r";
  unless (-e "$animal.mm9.methseqdepth.plus.bigwig") {
    print "cp $bigwigDir/$animal.$dates{$r}.mm9.plus.methseqdepth.bigwig $animal.mm9.methseqdepth.plus.bigwig &\n";
  }
  unless (-e "$animal.mm9.methseqdepth.minus.bigwig") {
    print "cp $bigwigDir/$animal.$dates{$r}.mm9.minus.methseqdepth.bigwig $animal.mm9.methseqdepth.minus.bigwig &\n";
  }
} } }