htseq-count returns : does not contain a 'gene' attribute
Dear BIOSTAR community,
I'm trying to make count matrix with htseq-count,
htseq-count -s yes -t gene -i gene 01.sorted.sam annotation_cattle.gff > 01.txt
even with --idattr=gene , it returns error:
Error processing GFF file (line 1864255 of file annotation_cattle.gff):
Feature gene-D1Y31_gp1 does not contain a 'gene' attribute
[Exception type: ValueError, raised in features.py:329]
The 1864255 line of my gff file is:
NC_039210.1 RefSeq gene 1095 8093 . + . ID=gene-D1Y31_gp1;Dbxref=GeneID:37627025;Name=D1Y31_gp1;gbkey=Gene;gene_biotype=protein_coding;locus_tag=D1Y31_gp1
Is there any solution that maybe I could solve it?
Best regards
• 1,176 views
•
link
written
by
ziziqolo •
1 answer
It solved by -i ID, and also I'm thinking what would happen if I enter gene name by myself in that line?
• 0 views
•
link
written
by
ziziqolo •
Log in to answer this question.
More posts like this
-
query regqrding genome feature extraction (gff/gtf) file extraction
written by microorganism_001 •
-
Read count using htseq-count.
written by Saeed •
-
htseq-count returns 0 count for all
written by ziziqolo •
-
Query on htseq count
written by lavanyac790 •
-
Query on htseq count
written by lavanyac790 •
-
[Title Edited] How to update the GFF3 wrong file? [was:] featureCounts may be ignore chloroplast se…
written by marcelolaia •
-
How to print lines of two file where column 2 and 3 of File 1 fall in between column 2 and 3 of fil…
written by Obex33 •
-
'CDS' but not 'exon' in GFF
written by NielQC •
-
split a GFF3 with many Genomes to GFF3 files with individual Genomes
written by Chris •
-
htseq-count error Feature does not contain a 'gene_id' attribute
written by madzayasodara •
