BedOverlapName: Difference between revisions
From genomewiki
Jump to navigationJump to search
No edit summary |
No edit summary |
||
| Line 2: | Line 2: | ||
#!/bin/bash | #!/bin/bash | ||
PREFIX=$3 | PREFIX=$3 | ||
overlapSelect $1 $2 -mergeOutput stdout -dropped=stdout| cut -f 1,2,3,4,8 | gawk '{OFS="\t";print $1,$2,$3,$4" | overlapSelect $1 $2 -mergeOutput stdout -dropped=stdout| cut -f 1,2,3,4,8 | gawk '{OFS="\t";print $1,$2,$3,$4",'$PREFIX'"$5}' | ||
</pre> | </pre> | ||
example: | |||
<pre> | |||
bedOverlapName randomFeatures.bed geneAnnotation.bed EXON_ | | |||
bedOverlapName stdin conservedRegions.bed CR_ | | |||
bedOverlapName stdin UTRs UTR_ | |||
</pre> | |||
will annotate the bed names with various additional information, each seperated by a comma | |||
[[Category:User Developed Scripts]] | [[Category:User Developed Scripts]] | ||
Revision as of 13:17, 29 September 2006
#!/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
