Error with HTseq RNAseq read count

Hi,

I am getting error while running HTseq. This is the command and the error:

htseq-count -q -f bam -s yes Ac1_mapped/ac1_mappedAligned.bam /global/home/users/catalinacastro/star/genome/genomic_v2.gtf > count.txt

Error occurred when processing GFF file (line 637338 of file /global/home/users/catalinacastro/star/genome/genomic_v2.gtf):
  not enough values to unpack (expected 9, got 1)
 [Exception type: ValueError, raised in __init__.py:221]

Thanks, Catalina

rnaseq htseq

• 739 views • link

updated 2.7 years ago by GenoMax

• written 2.7 years ago by Catalina •

Is that a blank line at end of your GTF file?

• 0 views • link

written 2.7 years ago by GenoMax

0 answers

No answers yet.

More posts like this

How to fix "ValueError: not enough values to unpack (expected 3, got 1)" when excuting batch analys…
written 6.0 years ago by mychung3265 •

CNVkit 0.9.7 Detected file format: bed Splitting large targets Wrote results/120430_HG19_ExomeV3_UTR_EZ_HX1_primary_annotated.target.bed with 407920 regions Detected file format: bed Traceback (most recent call last): File "/home/lab/miniconda3/envs/CNVKIT/lib/python3.7/site-packages/skgenome/tabio/util.py", …
htseq gene atrribute error
written 6.5 years ago by Balatheskulter •

$ htseq-count -q --additional-attr 534accepted_hits.bam 535accepted_hits.bam Rsolani.gff Error occured when processing GFF file (line 10 of file Rsolani.gff): Feature exon-AG1IA_00001-2 does not contain a 'gene_id' …
Error with HTseq RNAseq read count
written 7.0 years ago by Bioinfonext

Hi, I am getting error while running HTseq, could you please suggest what is the issue; I am using this command $ htseq-count Oryza_indica.ASM465v1.44.chr.gff3 Leaf_T1_F_R10_S1_L001.sam …
HTseq code error repeatedly
written 7.5 years ago by bnayer26 •

Hi, I am really new to bioinformatics so please help me figure out this error. I have tried looking at other threads with similar questions …
Error on HT-Seq-count while RNA quantification
written 8.5 years ago by k.kathirvel93

Hi all, I have started human RNA-seq data analysis recently. After the mapping of draft genome against the GRCh38 & GTF by using TopHat and …
Gff file error in HTSeq-count
written 9.2 years ago by y.gladbach •

Hi there, I need your help:). I guess my gff file is not correctly formatted, but not sure. Help is much appriciated. So please let …
HT-seq error and empty result
written 9.6 years ago by biologo •

hello, friends: i was using the htseq to calclate the reads, but it was sth wrong, like that: Error occured when processing SAM input (line …
HTSeq error when reading BAM file
written 9.9 years ago by Mike

Hi, I am running HTseq for read count with STAR output BAM file, but I got following errors: STAR --genomeDir /databank/igenomes/Homo_sapiens/STAR --readFilesIn R1.fastq --runThreadN 8 …
htseq-count error: "does not contain 'gene_id'
written 10.2 years ago by lynne.m.a.harris •

I'm attempting to use htseq-count to get read counts for some *Plasmodium *RNA-seq data. The command line looks as follows: htseq-count --format=bam --order=pos 3D7_1-1.bam PlasmoDB-28_Pfalciparum3D7.gff …
HTSeq analysis for Illumina tissue samples?
written 11.0 years ago by pbio

Hi, I have downloaded a few tissue bam files from illumina body map from the link ftp://ftp.ensembl.org/pub/release-81/bamcov/homo_sapiens/genebuild/ Now to get the read counts from each …