BedRegion

From genomewiki
Jump to navigationJump to search
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.
#!/usr/bin/gawk -f
# max 11/05
BEGIN { 
        if (ARGC>2) {
          print;
          print "Will display the minimum start and maximum end position of a bed file
          print "in UCSC-style format (chr4:minstartpos-maxendpos)"
          print ""
          print "SYNTAX:"
          print " bedRegion" 
          print "EXAMPLE:"
          print " cat bla.bed | bedRegion"
         exit 1

        }
         FS=" ";
         min=10000000000;
         max=0;
        chrom = "noChrom";
}
        
/track type=wiggle_0/ { wiggle=1;
    next;}
/track/ { wiggle=0; 
            next;}
/browser/ {next;}
/^#.*/ {next;}
/^$/ {next;}


//      { if (wiggle!=1) {
          if ($2 < min) {
           min = $2
          }
        if ($3 > max) {
                max = $3 
        }
        chrom = $1
    }
}

END { 
    if (chrom!="noChrom") 
        print chrom ":" min "-" max; 
    else {
        print "no chrom found, might be a wiggle track or empty";
        exit (1);
    }

}