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

featureCounts process fails when setting --fc_count_type to gene #440

Closed
pcantalupo opened this issue Jul 14, 2020 · 0 comments
Closed

featureCounts process fails when setting --fc_count_type to gene #440

pcantalupo opened this issue Jul 14, 2020 · 0 comments
Assignees
Labels
bug Something isn't working

Comments

@pcantalupo
Copy link
Contributor

I'm trying to run pipeline with a BK virus GTF file (contents of the GTF file are here) which does not contain exon as a count-type. Therefore, I set --fc_count_type to gene.

When I run:
nextflow run nf-core/rnaseq -profile docker --aligner hisat2 -r 1.4.2 --reads 'data/SRR11074360_{1,2}_10k.fastq' --reverseStranded --max_cpus 3 --max_memory 12GB -c config_bkv.txt --fc_count_type gene --fc_group_features gene_id --fc_group_features_type gbkey

The first featureCounts command executes properly:
featureCounts -a GCF_000837865.1_ViralProj14074_genomic.gtf -g gene_id -t gene -o SRR11074360_1_10k.sorted_gene.featureCounts.txt --extraAttributes gene_name -p -s 2 SRR11074360_1_10k.sorted.bam

But the second one featureCounts -a GCF_000837865.1_ViralProj14074_genomic.gtf -g gbkey -o SRR11074360_1_10k.sorted_biotype.featureCounts.txt -p -s 2 SRR11074360_1_10k.sorted.bam fails with the error: ERROR: no features were loaded in format GTF. The annotation format can be specified by the '-F' option, and the required feature type can be specified by the '-t' option... I think this is because the -t gene option is not added to the FC command line on this line: https://github.com/nf-core/rnaseq/blob/master/main.nf#L1401. Therefore FC rightly declares that it cannot find any "features".

@pcantalupo pcantalupo self-assigned this Jul 14, 2020
@ewels ewels added the bug Something isn't working label Jul 22, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

3 participants