Question: Obtaining coverage for rna-seq data
gravatar for rotemkat
6 weeks ago by
rotemkat10 wrote:

I have many bam files (containing rna-seq data), and i need to obtain the total coverage in each file. I tried to use "samtool depth", but since the output is so huge - and i need to sum up all the values there, i was wondering of there's a simpler/faster way to get this value.

ADD COMMENTlink written 6 weeks ago by rotemkat10

It does not make sense to calculate total coverage for each file since regions of your genome may be deeply covered in some places while others may have no reads. If you just need theoretical coverage then you could simply do: (total number of bases sequenced/size of the genome in bases).

You may want to try out mosdepth (LINK) and use one of the coverage options described there in (e.g. specific regions, chromosomes, BED intervals etc).

ADD REPLYlink modified 6 weeks ago • written 6 weeks ago by GenoMax95k
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: 1700 users visited in the last hour