Sequencing Depth (Read Depth) Calculations
1
0
Entering edit mode
14 days ago
R.L. • 0

Hi Community,

I have to calculate the read depth of HG002 UCSC Nanopore Sequencing data for my benchmarking.

I downloaded 3 fastq files provided in GIAB and mapped them with minimap2.

minimap2 -ax map-ont --secondary=no --MD

ChatGPT gives me following cmd to calculate:

sequencing_depth=$(samtools depth -a "$sorted_bam_file" | awk '{sum += $3} END {print sum/NR}')

echo "Sequencing Depth: $sequencing_depth"

This cmd works and gives me 40.5235, I was wondering whether this cmd is right or wrong, and does anyone have a sense on the correctness of this 40x depth?

Thanks a lot

depth • 292 views
ADD COMMENT
0
Entering edit mode

it looks looks ok to me.

ADD REPLY
0
Entering edit mode

Thank u very much!

ADD REPLY
0
Entering edit mode
14 days ago
GenoMax 142k

You could also use pandepth (LINK) or mosdepth (LINK) as well. These will allow you to be flexible about windowed coverage etc.

ADD COMMENT
0
Entering edit mode

Understood! Gonna try those two!

ADD REPLY

Login before adding your answer.

Traffic: 2803 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