2
$\begingroup$

I am following the tutorial here to learn how to process raw RNA-seq data and get gene read counts. I've had no issues until I got to step 3.9, "Assignment of RNA-Seq Reads to Genes." Here, the authors run the command

fc0 <- featureCounts(sam.list, annot.ext = 
"Ara-port11_GFF3_genes_transposons.201606.gtf", allowMultiOverlap = T, isPairedEnd = F, 
nthreads = 2, strandSpecific = 0, isGTFAnnotationFile = T)

I ran the same command except I set isGTFAnnotationFile = F because I couldn't find the gtf file, so I used the gff file TAIR10_GFF3_genes_transposons.gff instead, downloaded from here (they mentioned that using gtf only speeds up the process but that you can still do this using gff). I'm not sure if this is what is causing the issue, but the error I get is:

Error: Line 1 contains a format error. The expected annotation format is SAF.
Error in featureCounts(sam.list, annot.ext = "TAIR10_GFF3_genes_transposons.gff",  : 
  No counts were generated.
Improve this question
$\endgroup$
1
  • $\begingroup$ Thanks for answering the question @AnIgnorantWanderer please consider "Accepting" your answer below because it helps with site statistics. $\endgroup$
    – M__ ♦
    Sep 13 at 16:28

1 Answer 1

Reset to default
1
$\begingroup$

Answer from @an-ignorant-wanderer, converted from comment:

I ended up using the annotation file from https://support.illumina.com/sequencing/sequencing_software/igenome.html and it worked.

[Please update/improve this answer if possible]

Improve this answer
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.