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