Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

0 gene counts for each region #107

Open
Valentin-Bio-zz opened this issue Jun 16, 2022 · 0 comments
Open

0 gene counts for each region #107

Valentin-Bio-zz opened this issue Jun 16, 2022 · 0 comments

Comments

@Valentin-Bio-zz
Copy link

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.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant