Skip to content

Commit

Permalink
attempt #2 to fix ClinvarSignificance output
Browse files Browse the repository at this point in the history
  • Loading branch information
rsasch committed Jan 30, 2024
1 parent 3457c37 commit 5ef1c51
Showing 1 changed file with 3 additions and 3 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -1019,18 +1019,18 @@ task ClinvarSignificance {
risk factor
uncertain significance' > expected_clinvar_classes.csv

comm -23 <(sort bq_clinvar_classes.csv) expected_clinvar_classes.csv > missing_clinvar_classes.csv
comm -23 <(sort bq_clinvar_classes.csv) expected_clinvar_classes.csv > missing_clinvar_classes.txt

NUMRESULTS=$( wc -l bq_clinvar_classes.csv | awk '{print $1;}' ) # we expect this to be 13+
NUMMISS=$( wc -l missing_clinvar_classes.csv | awk '{print $1;}' ) # we expect this to be 0
NUMMISS=$( wc -l missing_clinvar_classes.txt | awk '{print $1;}' ) # we expect this to be 0

echo "false" > ~{pf_file}
# If the result of the query has fewer than 13 rows, that means clinvar_classification must not have all the expected values
if [[ $NUMRESULTS -ge 13 && $NUMMISS -eq 0 ]]; then
echo "The VAT table ~{fq_vat_table} has the correct values for clinvar classification" > ~{results_file}
echo "true" > ~{pf_file}
else
tr '\n' ', ' < missing_clinvar_classes.csv
tr '\n' ', ' < missing_clinvar_classes.txt > missing_clinvar_classes.csv
echo "The VAT table ~{fq_vat_table} has missing values for clinvar classification: " > ~{results_file}
cat missing_clinvar_classes.csv >> ~{results_file}
fi
Expand Down

0 comments on commit 5ef1c51

Please sign in to comment.