Hi,
I have 5 replicates of RNAseq data. I aligned the read files with STAR and then assembled the gtf (using StringTie) for each replicate individually giving me 5 gtf files. I then used Stringtie Merge on the 5 gtf files to get a single merge.gtf.
When I do wc -l on the individual gtf files, I get an average of 700,000+ transcripts. When I do the same for the merge.gtf, I get 1.7 million transcripts. Since these are all replicates, shouldn't I get roughly the same number of transcripts in the merged gtf and the individual gtfs?
Did I do something wrong in the merge step?
•
written
by
nattzy94 •
1 answer
The vast majority of them could be low abundance transcripts, i.e., artifacts of expression of neighbouring high abundance transcripts, otherwise known as 'transcriptional noise'. The possibility exists that they are also false-positive alignments from STAR. You could increase your thresholds when running stringtie merge - take a look at the bottom of this page:
If this does not help, please share your full STAR and StringTie commands.
Kevin
Log in to answer this question.
-
GTF file generated by stringtie merge has two source ensemble and stringtie
written by musk •
-
Clarification on StringTie's `--merge` parameters
written by DGTool •
-
when should I merge in annotation transdecoder pipeline
written by Mike •
-
StringTie merged transcripts for two different conditions
written by vitor •
-
Error when running Stringtie Merge: Error getting header ID# for gseq 1
written by nattzy94 •
-
Question about gffcompare using multiple reference annotations from different databases
written by newbie
-
Checking workflow for assembly of time series expt using Stringtie-Ballgown
written by nattzy94 •
-
Transcriptomic assembly Stringtie
written by Nelislam4 •
-
Expression data from cufflinks merged gtf
written by ddowlin
-
Different results when visualizing tdf files in IGV and when doing the differential expression ana…
written by iraia.munoa
