Wednesday, January 8, 2014

map bed file to gene(another bed)

sort-bed hsa.bed > hsa.bed.2

sort-bed miRNA/GSM1097888_CD34_smallRNA-seq.bed > miRNA/GSM1097888_CD34_smallRNA-seq.bed.2

bedmap --skip-unmapped --echo --echo-map --count --sum hsa.bed.2 miRNA/GSM1097888_CD34_smallRNA-seq.bed.2 > b.bed.2

or

bedmap --skip-unmapped --echo --max hsa.bed.2 GSM973662/miRNA.bed.2 > a.bed.2

No comments:

Post a Comment