zero counts for all genes in RNAseq data of Ferret
Entering edit mode
7 months ago
Sara ▴ 260

I have bulk RNA-seq data from Ferret and trying to get counts per gene. to do so I used hisat2 and got the genome from here:

After aligning the fastq files I used htseq and the following command:

htseq-count \
  -f bam \
  -r pos \
  -a 27 \
  -m intersection-strict \
  --stranded=no \
  -i gene_id \
  test_R1_001.fastq.sam.bam Mustela_putorius_furo.MusPutFur1.0.110.gtf > test_R1_001.fastq.sam.bam_sorted.bam.txt

with the gtf file from here:

but all the counts are 0. do you know what the problem could be?

RNA-seq DGE • 408 views
Entering edit mode
7 months ago
GenoMax 143k

got the genome from here:

Just to be certain. You did convert the 2bit file to fasta sequence using twoBitToFa and then built a HiSat2 index with the fasta? Matching GTF files from UCSC should be in

You are almost certainly going to run into trouble trying to mix and match sources of genome and annotation.

If you are getting the GTF from Ensembl then you could get the genome from the same location. That way you are assured that the identifiers will match between genome and GTF. That is one of the top causes of problems with counting. Use


Login before adding your answer.

Traffic: 1405 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6