We have to find the expression level of a gene on different conditions. We have assembled both RNA-seq sample data using velvet-Oases. Now I need to find the gene from a full length CDNA library using the assembled contigs. Then I need to map the reads with the gene to find the expression level. I think I need to use BLAST for finding gene from cDNA library and then some alignment tool for finding expression level. I am just wondering if I am in right direction. If not, can you please suggest what should I do. Also what are the best tools I could use for this. I am just new to this bioinformatics, so need some help from experts.