Question: How to calculate per base coverage from the .bam file?
1
gravatar for barrypraveen
9 months ago by
barrypraveen30
barrypraveen30 wrote:

Hi all, I have some genes which are expressed in differentially compare to wild type, so i want to calculate perbase coverage from the BAM files. How to calculate that one.

rna-seq next-gen R • 533 views
ADD COMMENTlink modified 9 months ago by ATpoint36k • written 9 months ago by barrypraveen30
1

samtools depth but may I ask why you want to do that?

ADD REPLYlink written 9 months ago by ATpoint36k

I have to plot that expression. I have RPKM values for that gene but if I plot in the excel file it comes like a point. my idea is like showing the distribution of the expression. lets say x-axis is a length of the gene in the y-axis expression of the gene.

ADD REPLYlink modified 9 months ago • written 9 months ago by barrypraveen30

I do not really see how this would make sense. Gene coverage in RNA-seq is the sum of transcript/exon coverage so you cannot easily get a representative per-base coverage. Yes, gene expression is commonly reported as one value per gene, be it the RPKM metric or any other normalized expression value. What do you think would you gain from plotting a per-base coverage of a gene? Are you referring to exon usage?

ADD REPLYlink written 9 months ago by ATpoint36k

actually it is not a gene. it is a snoRNA. Just for saying i wrote as a gene. but when i try samtool depth

samtools depth -b missed_bed.bed -f PC1_Yale_PRS_sorted.bam

i got error

"Does the file "PC1_Yale_PRS_sorted.bam" really contain a list of files and do all exist?"

I really don't know???

ADD REPLYlink modified 9 months ago • written 9 months ago by barrypraveen30
1
gravatar for ATpoint
9 months ago by
ATpoint36k
Germany
ATpoint36k wrote:

Purely from the technical side the command is samtools depth -b missed_bed.bed PC1_Yale_PRS_sorted.bam without -f. Explore if -a might help as well.

ADD COMMENTlink written 9 months ago by ATpoint36k

got it....now its fine

ADD REPLYlink written 9 months ago by barrypraveen30
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 2.3.0
Traffic: 1358 users visited in the last hour