I'm trying to run htseq-count with Bowtie alignments to a bacterial genome. Even though almost 100% of reads aligned to the genome, when I run htseq-count all I get is 0 for every gene. My genome fasta file used for alignment as well as annotation gff were both downloaded from Ensembl and should be correct format. I checked chromosome names in both alignment and annotation as suggested in another post on the subject but that did not fix the problem (there is only one chromosome, so just labeled as "chromosome"). Any other tricks to try? I'm convinced it is just a formatting issue of some sort.
Heads up! This is a static archive of our support site. Please go to help.galaxyproject.org if you want to reach the Galaxy community. If you want to search this archive visit the Galaxy Hub search
Question: htseq-count returning only 0
17 months ago by
lpb • 0
lpb • 0 wrote:
ADD COMMENT • link •
Please log in to add an answer.
Powered by Biostar version 16.09
Traffic: 166 users visited in the last hour