Question: zeroes in ht-seq output
0
gravatar for komal.yasmin
3.4 years ago by
komal.yasmin10
komal.yasmin10 wrote:

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 • 1.1k views
ADD COMMENTlink modified 3.2 years ago by Biostar ♦♦ 20 • written 3.4 years ago by komal.yasmin10
2

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

ADD REPLYlink written 3.4 years ago by Devon Ryan92k

how do I confirm if they have different chromosome names?

ADD REPLYlink written 3.4 years ago by komal.yasmin10
cut -f 1 foo.gtf | sort | uniq -c
samtools view -H foo.bam
ADD REPLYlink written 3.4 years ago by Devon Ryan92k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1684 users visited in the last hour