I have used StringTie on mice RNA-seq data. StringTie was run with parameters given in Nature paper. Reference gtf annotation was downloaded from the Gencode Genes.
Most of the transcripts in the reference annotation have very low expression values (> 95% transcripts have 0 fpkm; for one sample all reference transcripts have 0 fpkm). However, half of newly merged transcripts have non-zero expression values. My questions are:
- Is it normal for all the reference transcripts to have expression values equal to zero?
- After running StringTie with default parameters I got many (maybe too many) newly merged transcripts (140K vs 25K in the reference annotation). This looks kinda strange and I would prefer to use transcripts from the reference, but as I mentioned above - they have 0 expression values.