Question: featureCounts for multiple features
0
gravatar for hwu12
12 months ago by
hwu1210
hwu1210 wrote:

Hi, I am using the linux version of featureCounts, and I would like to count the number of reads in CDS, 5'UTR and 3'UTR. I was wondering if I should just count three times like:

featureCounts -t CDS -g gene_id -a $GTF -o counts.txt $BAM

featureCounts -t five_prime_UTR -g gene_id -a $GTF -o counts.txt $BAM

featureCounts -t three_prime_UTR -g gene_id -a $GTF -o counts.txt $BAM

Or is there a way to do all three together? Thanks!!!

 

 

 

 

featurecounts • 670 views
ADD COMMENTlink modified 12 months ago by Wei Shi3.1k • written 12 months ago by hwu1210
Answer: featureCounts for multiple features
3
gravatar for Wei Shi
12 months ago by
Wei Shi3.1k
Australia
Wei Shi3.1k wrote:

Regarding questions related to featureCounts that is released as part of the SourceForge Subread package, please send your question to Google Subread group (https://groups.google.com/forum/#!forum/subread).

You cannot count reads to three different types of features at one time with your current GTF annotation file. However you can get around this by changing the names of the three feature types to the same name, for example you may rename all of them to 'CDS_or_5UTR_or_3UTR'. Then you can count reads to all the three types in one command:

featureCounts -t CDS_or_5UTR_or_3UTR -g gene_id -a $GTF -o counts.txt $BAM

But this might take more effort than simply summing up the counts you got from the three individual counting.

ADD COMMENTlink modified 12 months ago by Gordon Smyth37k • written 12 months ago by Wei Shi3.1k

Thanks so much, Wei and Gordon!!!

ADD REPLYlink modified 12 months ago • written 12 months ago by hwu1210
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 261 users visited in the last hour