You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
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".
The text was updated successfully, but these errors were encountered:
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
togene
.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".The text was updated successfully, but these errors were encountered: