Thanks Mete and Jenifer for your information.
Last time, I did mRNA sequencing analysis and decided to go with TopHat>Htseq>DESeq for DE genes. Just because the results match up quite nicely with my qPCR validation. Although Cuufflink/Cuffdiff produced results with quite similar trend with DESeq, It seem to me that Cuffdiff tend to inflate the fold-change and is not good at statistical analysis. Thanks Jenny and Ross.
Anyway, about the miRNA I am working on now. My miRNA data is 51bp, single-ended. I am going to cut adapter using FASTX-toolkit and align reads using Novoalign ( as Mete suggests) and Bowtie.
Just have a question for now: my 3' adapter sequence is : 5-rAppAGATCGGAAGAGCACACGTCT-NH2-3. How many bases should I put in the " Enter custom clipping sequence" ? Is " AGATCGGA" is optimal? Just because I observed that many reads only contain part of 3' adapter sequence.
Also, Do I need to trim 5' adapter as well? and How?
Thank so much for your help
Thanh