[BioC] Amplicon and exon level read counts and GC content

Yu Chuan Tai yuchuan at stat.berkeley.edu
Wed May 30 08:00:02 CEST 2012

Hi Wei,

Thanks a lot for your quick reply. I will look into featureCounts in 

Yu Chuan

On Wed, 30 May 2012, Wei Shi wrote:

> Dear Yu Chuan,
> You may consider using the featureCounts function in Rsubread package to count the number of reads for the features you are interested in. The built-in RefSeq annotation will allow you to get exon-level read counts. But to get amplicon-level read counts or read counts in a specific genomic interval, you'll have to provide an annotation file to this function to let it count reads for you.
> This function accepts SAM files as input, but I guess you can convert your BAM files to SAM files.  For you information about this function, see ?featureCounts.
> Hope this helps.
> Cheers,
> Wei
> On May 30, 2012, at 3:17 PM, Yu Chuan Tai wrote:
>> Hi,
>> I have some questions about DNA-Seq and RNA-Seq analyses. In Amplicon sequencing, is there any package/function which can extract amplicon-level read counts and GC content  from an aligned file of BAM format? The same question to exon-level read counts and GC content. More generally, given a genomic interval, how could I extract the read count and GC content for that interval?
>> Thanks for any input!
>> Best,
>> Yu Chuan
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
> ______________________________________________________________________
> The information in this email is confidential and inte...{{dropped:5}}

More information about the Bioconductor mailing list