Skip to content

0 gene counts for each region #107

Description

@Valentin-Bio-zz

Hello, I'm trying to count the gene abundance from Metagenomic libraries using htseq-count program
from blast otuput format table 6 I converted it into a gff file with blast2gff, after that I used gat program
agat_convert_sp_gff2gtf.pl to convert the gff to gtf. This is how the gtf file looks like:

Captura de Pantalla 2022-06-16 a la(s) 13 39 43

and this is how one of the htseq-count output table looks like using a bam file

Captura de Pantalla 2022-06-16 a la(s) 13 44 36

It seem something is missing on the gtf table but I cannot figure out what is lacking, this is the code that I used:

htseq-count -r pos -n 9 -t match_part -f bam sample1.bam sample_1.gtf > sample1_genes.count

So why it outputs "0" values of counts ? even the first column of the htseq-count program provide "nbisL1" names that I don't know to what they refer to.

Thanks for your time :)

Bests,

Valentín.

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions