Hi all,
I am using featureCounts to quantify my RNA-seq signal over gene regions. I am using it as follows:
For my single-end reads:
featureCounts -t exon -g gene_id -s 0 -O -T 10 -a <.gtf> -o <.counts> <.bam>
For my paired-end reads: featureCounts -p -t exon -g gene_id -s 0 -O -T 10 -a <.gtf> -o <.counts> <.bam>
Where I am running into something confusing is with the -s option. I have run these with all options for -s, with differing proportions of 'successfully assigned reads'. -s 0 (unstranded) and -s 2 both assign about ~25% of the reads, however, -s 1 assigns only ~3% of reads to gene regions. The fact that both 0 and 2 are producing similar proportions of assigned reads/fragments is interesting, and must say something about the library... but I need help determining exactly what that is.
Thanks in advance. Let me know if any additional information is needed.
•
written
by
Noah •
1 answer
-s 2 captures sense transcription,
-s 1 captures the antisense transcription.
-s 0 captures both
so it seems you have low-levels of antisense transcription,
everything is working as expected
Log in to answer this question.
-
Featurecounts has more reads than star assigned reads ?
written by PK •
-
Counting gene vs exons in featureCounts
written by Meisam •
-
Exon-level counts with featureCounts
written by rependo •
-
Why there are no reads for gene quantification but many reads can be quantified for exon using feat…
written by mina95921m •
-
Different results from featureCounts using gene_id and gene_biotype
written by CuriusScientist
-
Very low percentage of assigned reads using featureCounts
written by Meghdad •
-
Illumina hiseq 2500 strand specificity for paired-end datasets
written by nnry06 •
-
featureCounts in predicting strand specificity
written by vrehaman •
-
featureCounts segmentation fault, Is there any solution yet ?
written by EagleEye
-
Cufflinks genome guided de novo transcript assembly- increase in ambiguous reads
written by sunil.mangalam •
Your library is reverse stranded, so if you set it to stranded (
-s 1) you would be aligning most of your reads to the wrong strand relative to the gene coding direction. With unstranded (-s 0) you align to both strands so you would capture any transcripts from that section of DNA regardless if the library was unstranded, stranded,or reverse stranded.Another important detail when counting paired end reads is that the meaning of
-phas changed:basically when read pairs are counted the counts will be halved