Question: How to get output with read counts and read count percentages from bam file?
gravatar for bayramliaziz
15 months ago by
bayramliaziz0 wrote:

I have a bam file with lots of genes included and I would like get read counts and their percentages as a bam file. I tried samtools idxstats but It doesn't give percentages. Is there any tool that can do what I want? If there is I would like to get you opinion. Also my bam file is derived from pair-end reads.

Thanks in advance.

bam samtools picard read counts • 559 views
ADD COMMENTlink modified 15 months ago by WouterDeCoster44k • written 15 months ago by bayramliaziz0
gravatar for WouterDeCoster
15 months ago by
WouterDeCoster44k wrote:

I would suggest using featurecounts to count the number of reads per gene, and calculate the percentages yourself, e.g. using R or Python.

It is not entirely clear whether you aligned to the genome or to the transcriptome, or how you obtained the sequencing data. Often alignment to the genome is the best approach to avoid a bias on (incomplete?) annotation. But since you state that samtools idxstats only gave you the counts it looks like you aligned to a transcriptome, in which case my answer may not be exactly what you need. It is also unclear what you aim to achieve in the long run, so we cannot comment if this is the best strategy.

ADD COMMENTlink written 15 months ago by WouterDeCoster44k

I aligned to the genome by using bwa aln after that I used samtools to get bam file. I am new to this thing so I am just trying to explore and experiment with different approaches. I'll check out the featurecounts.

ADD REPLYlink written 15 months ago by bayramliaziz0
Please log in to add an answer.


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