This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to get coverage ratio?

Dear colleagues:

I have got whole exome sequencing data of dozens of samples. I want to know how we can get the exome data capture quality index of coverage ratio.

Any help that can help me to understand the question would be appreciated.

Best!

sequence next-gen-sequencing
• 3,059 views • link
updated by Ram Ram • written by 897598644 897598644

1 answer

Sounds like a job for Picard CalculateHsMetrics.

Run it on each bam file, then check the column PCT_SELECTED in the output metrics file for each sample.

• 0 views • link
updated by Ram Ram • written by Danielk Danielk

Are there any examples of command lines to show how to achieve it?

• 0 views • link
updated by Ram Ram • written by 897598644 897598644

It's something like

java -jar CalculateHsMetrics.jar INPUT=in.bam METRICS=metrics.txt

but do check the details of the syntax on their website, it changed recently so the above might not be correct (but correct-ish).

• 0 views • link
updated by Ram Ram • written by Danielk Danielk

What would be the metrics.txt file look like?

• 1 views • link
updated by Ram Ram • written by MAPK MAPK

Log in to answer this question.

More posts like this