I used the below-mentioned command using ncbi-blast+ in order to find the conserved miRNA. I have created the miRbase database and ran the below-mentioned command. I needed read_count and description of the aligned sequence in the output file. I will really appreciate if you could kindly provide your suggestion in this regard.
blastn -db <db_source> -query <query_source> -out <outfile> -outfmt "6 qseqid sseqid slen qstart qend length mismatch gapopen gaps sseq" -word_size 4 -perc_identity 100