Entering edit mode
4 weeks 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
you use the term coverage and depth as if they are different things but (to my knowledge) they mean the same thing
coverage refers to the percentage of the genome (or target region) that is sequenced at least once, while depth indicates the number of times a particular base is sequenced.
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
thank you for answering!
even though this parameter shows base-by-base depth from samtools depth, can I say it also reflects coverage across the gene? Since depth per base shows which regions are well or poorly covered, is it correct to make that inference?