Entering edit mode
7.3 years ago
mkh
▴
60
Hi all
I am trying to do genome assembly by stringtie software. The command that I am using is as follow:
" stringtie -p 2 -G dmel.gtf -o out.gtf STAR.bam "
The input bam file is from STAR.
The output file is look like this:
3R StringTie transcript 1 656 1000 - . gene_id "STRG.1"; transcript_id "STRG.1.1"; cov "3.112805"; FPKM "382.908844"; TPM "572.137207";
3R StringTie exon 1 656 1000 - . gene_id "STRG.1"; transcript_id "STRG.1.1"; exon_number "1"; cov "3.112805";
3R StringTie transcript 734 1776 1000 - . gene_id "STRG.2"; transcript_id "STRG.2.1"; cov "6.242569"; FPKM "767.903931"; TPM "1147.391479";
3R StringTie exon 734 1776 1000 - . gene_id "STRG.2"; transcript_id "STRG.2.1"; exon_number "1"; cov "6.242569";
So as it shown above there is no "transcript_id" and "gene_id" corresponding to those in the input gtf file. Stringtie end with no Error so i can not understand what is wrong with my inputs.
Thanks
See this: C: Help regarding RNA-seq STRINGTIE pipeline and output