Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

.DepthOfCoverage_Target.txt has odd output #90

Open
vinjana opened this issue Jan 31, 2024 · 0 comments
Open

.DepthOfCoverage_Target.txt has odd output #90

vinjana opened this issue Jan 31, 2024 · 0 comments

Comments

@vinjana
Copy link
Member

vinjana commented Jan 31, 2024

The file ${sampleName}_${pid}_targetExtract.rmdup.bam.DepthOfCoverage_Target.txt is produced by the coverageQc D program and has somewhat odd output.

  1. The normal chromosomes are expressed as overall coverage, i.e. all (non-N) positions in the genome are considered for both enumerator and denominator of the coverage value. That does not fit well to the "target"-specificity.
  2. The "all" row in the output is restricted to the target regions, but because coverageQc is called on the complete file, instead of on the target-region specific file, the enumerator uses all bases (in- and off-target), while the denominator is the target regions size. If the target region size is not provided by the user to exclude Ns, then both enumerator and denominator are not as expected.

In a new version, make sure that these issues are fixed.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant