Hi Yan,
You may use the HTseq count wrapper in the
http://galaxy.nbic.nl/
.
It does a good job and I could employ edgeR on that count matrix.
Good luck.
Best wishes,
Anto