BedOverlapName: Difference between revisions
From genomewiki
Jump to navigationJump to search
No edit summary |
No edit summary |
||
| Line 1: | Line 1: | ||
Add overlapping bedname from second bed: | |||
<pre> | <pre> | ||
#!/bin/bash | #!/bin/bash | ||
Revision as of 14:08, 2 October 2006
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
