Number of reads

Question

[isantanna]

Hey everyone,

I have one problem in my fpkm normalization. So I would like to know the number of reads per genes as we have in cufflinks output. Anyone knows how to do this? Is it possible do this with samtools?

Thank you

Answer 1

[Carlo Yague]

Try featureCounts from the alignment files. It does exactly what you want.