Question: How do I summarize coverage across many genomes to inform filtering cutoffs?
0
gravatar for selplat21
5 months ago by
selplat2110
selplat2110 wrote:

Hello,

I used samtools depth to get a textfile of coverage across my genome for ~600 samples. In order to inform coverage cutoffs in the next step I'd like to make a summary/histogram of some sort that summarizes coverage for all of these genomes together.

For example, eventually I'd like to say, filter out reads below x coverage and above y coverage. This would be based on the distribution of coverage for all samples to filter out likely duplicates etc.

Any help is appreciated!

ADD COMMENTlink modified 5 months ago by onestop_data260 • written 5 months ago by selplat2110
0
gravatar for onestop_data
5 months ago by
onestop_data260
onestop_data260 wrote:

Please check this out Filter Bam File Based On Coverage and also take a look int the new samtools named coverage (samtools coverage). Please sure you have the latest samtools.

ADD COMMENTlink written 5 months ago by onestop_data260

Do you mean samtools stats --coverage? I only see this as a subcommand in the new samtools. This is very helpful thank you!

ADD REPLYlink written 5 months ago by selplat2110

No. samtools coverage It was introduced here (samtools 1.10)

ADD REPLYlink written 5 months ago by onestop_data260
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1226 users visited in the last hour