SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   RNA Sequencing (http://seqanswers.com/forums/forumdisplay.php?f=26)
-   -   StringTie merge- can someone help me understand its output? (http://seqanswers.com/forums/showthread.php?t=88406)

BioLion 03-27-2019 05:20 AM

StringTie merge- can someone help me understand its output?
 
Hello,

I am currently using StringTie merge on a StringTie output in order to discover genes that are not already in my reference annotation. I tried looking at the documentation but it is not very clear to me how to distinguish some things in the output or how it proceeds.

The command I have been using is :
Code:

${stringtie} --merge -G reference.gtf -o genome.StringTie.merged.gtf -F 1 -l GenomeX GenomeX.StringTie.sorted.gtf
with GenomeX.StringTie.sorted.gtf the output of stringTie.

I do not exactly understand how it works: when it finds a region which maps or overlaps, does it replace the line from "GenomeX..." by the line from "reference.gtf"?
Also, what is the difference with this command?

Code:

${stringtie} --merge -o genome.StringTie.merged.gtf -F 1 -l GenomeX listGenomes.gtf
with listGenomes.gtf being:
"reference.gtf
GenomeX.StringTie.sorted.gtf"

When I looked at the output of the second command, I did not see how we can distinguish between "reference.gtf" lines and "GenomeX.StringTie.sorted.gtf", or which field I am supposed to look at. Is there a way to do so after a merge?Also, does it choose between the two files when a region is similar?


Thanks a lot for reading me,

BioLion


All times are GMT -8. The time now is 12:55 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.