Skip to content

Conversation

VJalili
Copy link
Member

@VJalili VJalili commented Aug 29, 2025

This PR fixes the medianCoverage.R script, which fails to report the sample ID in the output.

 > printf "Chr\tStart\tEnd\tNA12878\nchr1\t10000\t10100\t1014\nchr1\t10100\t10200\t378\nchr1\t10200\t10300\t218\n" > test_input.bed

Current output:

> Rscript src/WGD/bin/medianCoverage.R test_input.bed -H test_output.bed
> cat test_output.bed

#sample_id	Med_withoutZeros
378	NA

Note that the sample ID is missing.

Corrected output:

> Rscript src/WGD/bin/medianCoverage.R test_input.bed -H test_output.bed
> cat test_output.bed

#sample_id	Med_withZeros	Med_withoutZeros
NA12878	378	NA

Testing

Successfuly tested EvidenceQC, the output medianCov.transposed.bed looks as expected.

gs://fc-bd728187-9620-4fee-af18-967e29390e66/submissions/69e133a4-6d19-42a0-bfce-340fdf13e895/EvidenceQC/150a181f-9d12-45cf-b81a-59cdaa5e135a/call-MedianCov/MedianCov/e7667941-a396-4f15-95fb-39ba5d5a64e0/call-CalcMedCov/all_samples_medianCov.transposed.bed

@VJalili VJalili requested a review from mwalker174 August 29, 2025 18:35
@VJalili VJalili merged commit f091af0 into broadinstitute:main Sep 2, 2025
4 checks passed
@VJalili
Copy link
Member Author

VJalili commented Sep 2, 2025

Thanks, @mwalker174!

@VJalili VJalili deleted the fix-median-cov branch September 2, 2025 19:25
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

Successfully merging this pull request may close these issues.

2 participants