Coverage depth using samtools
0
0
Entering edit mode
1 day ago
slzr_ • 0

Is it correct to use the parameter below to evaluate coverage base per base or is it only going to provide depth?

  samtools depth -a -b

samtools ngs depth sequencing coverage • 167 views
ADD COMMENT
0
Entering edit mode

yes but the option -b is missing a file if you want to compute all bases overlapping a bed file.

https://www.htslib.org/doc/samtools-depth.html

-a

    Output all positions (including those with zero depth) 
-a -a, -aa

    Output absolutely all positions, including unused reference sequences. Note that when used in conjunction with a BED file the -a option may sometimes operate as if -aa was specified if the reference sequence has coverage outside of the region specified in the BED file. 
-b FILE

    Compute depth at list of positions or regions in specified BED FILE. [] 
ADD REPLY
0
Entering edit mode

you use the term coverage and depth as if they are different things but (to my knowledge) they mean the same thing

ADD REPLY

Login before adding your answer.

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