Extract reads counts with featureCounts

0

Hello all,

I'm trying to extract reads counts with featureCounts.

I obtained the gtf file from UCSC genome browsers through the command:

wget https://hgdownload.soe.ucsc.edu/goldenPath/dm6/bigZips/genes/dm6.ensGene.gtf.gz

I ran the following command to extract reads counts

featureCounts -a dm6.ensGene.gtf -o Male_RNAseq.tx sorted.files/*bam

However, I got error that says

Failed to open the annotation file dm6.ensGene.gtf, or its format is incorrect, or it contains no 'exon' features.

Does anyone know how I can fix this error? I greatly appreciate any help, thank you


RNAseq


BAM


subread


featureCounts

• 192 views

updated 2 hours ago by

100k

written 10 hours ago by

0



Source link