Obtaining coverage for rna-seq data
0
0
Entering edit mode
18 months ago
rotemkat ▴ 10

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.

samtools bedtools RNA-Seq variant calling • 418 views
ADD COMMENT
0
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 1576 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6