I have got RNA-seq reads from Illumina which I've trimmed adapters using trimmomatic and then aligned it with reference genome using hisat2. Now I've run through string tie for transcript identification and quantification. I would like to use ballgown for differential gene expression determination. Please advise if the pipeline I 've followed so far is right. also, I would like to use Ballgown for differential gene expression determination. Can anyone advise me on a good documentation for ballgown usage?