zeroes in ht-seq output
0
0
Entering edit mode
7.9 years ago
komal.yasmin ▴ 10

I am using ht seq to analyse RNA- seq data for Drosophila. I have used GFF3 file as well as GTF file in the script. I do not get an error report but the output file generated only contains gene IDs against zeroes for feature counts . Can someone please tell me what am I doing wrong? I am using the following script:

htseq-count --order=name --type=gene --idattr=ID --stranded=no --mode=union $ALIGNMENT_FILE $GFF_FILE > $EXPRESSION (Idattr is ID for Gff file and gene_id for gtf file)

RNA-Seq • 2.2k views
ADD COMMENT
2
Entering edit mode

Odds are good that your GTF and BAM file have different chromosome names.

ADD REPLY
0
Entering edit mode

how do I confirm if they have different chromosome names?

ADD REPLY
0
Entering edit mode
cut -f 1 foo.gtf | sort | uniq -c
samtools view -H foo.bam
ADD REPLY

Login before adding your answer.

Traffic: 1977 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