How can I compute gene expression for a set of RNA reads?
It's unclear if your paired-end reads are actually in fasta format, I'll presume that they're in fastq instead.

The easiest tool to use is salmon, which nicely deals with things like multimapping. If you're trying to judge the quality of an assembly, then I recommend having a look at transrate, which uses some related methodology for assessing contig quality.