grep mir ../ucsc.gtf|awk '{print $1"\t"$4"\t"$5"\t"$10"\t"$6"\t"$7}' > miRNAcoor.bed
sort -k1,1b -k2,2n miRNAcoor.bed > miRNAcoor.bed.s

perl makeBed.pl danRer10/danRer10.fa danRer10 > danRer10/danRer10.bed

bam2x getseq -i danRer10/miRNAcoor.bed -o miRNA/miRNA.fa -m seq -g ../danRer10.2bit -b 6
