## get mirbase18 annotation and convert to a bed like format
wget ftp://mirbase.org/pub/mirbase/CURRENT/genomes/hsa.gff
grep -v "#" hsa.gff | cut -f 1,4,5,9 >mirbase18.temp1
grep -v "#" hsa.gff | cut -f 6,7 >mirbase18.temp2
paste mirbase18.temp1 mirbase18.temp2 | sed s/^/chr/ >mirbase18.bed
Extra credit: Is there a clever shell trick to piping two separate stdins into paste, thus avoiding the temp files?
[...] note: Using the <(command) trick, my post from March 16 can now be [...]
ReplyDelete