BedOverlapName

From genomewiki
Revision as of 14:08, 2 October 2006 by Max (talk | contribs)
Jump to navigationJump to search

Add overlapping bedname from second bed:

#!/bin/bash
PREFIX=$3
overlapSelect $1 $2 -mergeOutput stdout -dropped=stdout|  cut -f 1,2,3,4,8 | gawk '{OFS="\t";print $1,$2,$3,$4",'$PREFIX'"$5}'

example:

  bedOverlapName randomFeatures.bed geneAnnotation.bed EXON_  |
  bedOverlapName stdin conservedRegions.bed CR_ |
  bedOverlapName stdin UTRs UTR_

will annotate the bed names with various additional information, each seperated by a comma