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.
Question: htseq-count returning only 0
16 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: 160 users visited in the last hour