Error when running Stringtie Merge: Error getting header ID# for gseq 1
0
0
Entering edit mode
19 months ago
nattzy94 ▴ 50

Hi,

I am trying to merge 3 GTFs that have assembled using stringtie --merge. I ran the following command:

/home/e0470749/stringtie/stringtie --merge -G /home/e0470749/reference_genome_ensembl/Homo_sapiens.GRCh38.99.gtf -o /gpfs/eplab/Nathaniel/HPAF2_celllines_star-stringtie/stringtie-merge/HPAFII_merged.gtf -p 4 /gpfs/eplab/Nathaniel/HPAF2_celllines_star-stringtie/stringtie/assembly_GTF_list.txt

However, I get the following error:

Error getting header ID# for gseq 1 (file: /gpfs/eplab/Nathaniel/HPAF2_celllines_star-stringtie/stringtie/H2_Veh24hr-1/H2_Veh24hr-1.gtf)

These are the first few lines of H2_Veh24hr-1.gtf:

# /home/e0470749/stringtie/stringtie H2_Veh24hr-1_Aligned.sortedByCoord.out.bam -p 12 -G /home/e0470749/reference_genome_ensembl/Homo_sapiens.GRCh38.99.gtf --rf -o /gpfs/eplab/Nathaniel/HPAF2_celllines_star-stringtie/stringtie/H2_Veh24hr-1/H2_Veh24hr-1.gtf -v -eB # StringTie version 2.1.2 1 StringTie transcript 1167104 1167198 1000 + . gene_id "ENSG00000207730"; transcript_id "ENST00000384997"; ref_gene_name "MIR200B"; cov "1.000000"; FPKM "0.164038"; TPM "0.379526"; 1 StringTie exon 1167104 1167198 1000 + . gene_id "ENSG00000207730"; transcript_id "ENST00000384997"; exon_number "1"; ref_gene_name "MIR200B"; cov "1.000000"; 1 StringTie transcript 1167863 1167952 1000 + . gene_id "ENSG00000207607"; transcript_id "ENST00000384875"; ref_gene_name "MIR200A"; cov "0.866667"; FPKM "0.142166"; TPM "0.328922";

My assembly_GTF_list.txt looks like this:

/gpfs/eplab/Nathaniel/HPAF2_celllines_star-stringtie/stringtie/H2_Veh24hr-1/H2_Veh24hr-1.gtf /gpfs/eplab/Nathaniel/HPAF2_celllines_star-stringtie/stringtie/H2_Veh24hr-2/H2_Veh24hr-2.gtf /gpfs/eplab/Nathaniel/HPAF2_celllines_star-stringtie/stringtie/H2_Veh24hr-3/H2_Veh24hr-3.gtf

I have run stringtie merge a few times before and it has always worked. Can someone help me to solve this error?

assembly stringtie • 421 views
ADD COMMENT

Login before adding your answer.

Traffic: 2056 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6