#!/usr/bin/perl

$ls = `ls -1 *.bam`;
foreach $bamfile (split/\n/, $ls) {
  unless (-e "$bamfile.bai") { system "samtools index $bamfile"; }
  @ar = split/\./, $bamfile; $id = $ar[0];
  if ($bamfile =~ /plus/) { $str = "plus"; } else { $str = "minus"; }
  print "track type=bam name=\"$id $str\" db=hg19 bamColorMode=strand pairEndsByName=yes visibility=full bigDataUrl=http://snpinfo.niehs.nih.gov/ucscview/fessler/GTEx/devSara/$bamfile\n";
}