MikeAxtell / ShortStack

ShortStack: Comprehensive annotation and quantification of small RNA genes
MIT License
88 stars 29 forks source link

differential gene expression #87

Closed mustafaraza1987 closed 3 years ago

mustafaraza1987 commented 5 years ago

Hello Dear Sir, Can you please explain about the analysis of differential gene expression of miRNA (replicate data) with short stack. I already mine miRNA from different replicates, but now issues is that how to create count matrix from count data as well as rpm , Each replicates provide the clusters different in number. Suppose replicate_1 =500 ,replicate_2=200, replicate_3=100. How to create count matrix in such a way .

groverj3 commented 5 years ago

I'm not Mike, but I think it sounds like you've run ShortStack on all your libraries independently. If you do that, it only uses that particular replicate's reads to generate loci. In order to have loci that are defined with all your samples, such that every sample has the same number of loci defined, with a value in the counts.txt file you just need to run them all together.

So, either provide all the .fastq files to ShortStack at once, or map them all separately and use samtools merge to generate a single .bam file to feed into ShortStack to define loci.