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

Compression integrity checks #93

Merged
merged 10 commits into from
Dec 7, 2023
Merged

Conversation

yashpatel6
Copy link
Collaborator

Description

Adding integrity check option for compressed files.


Test Results

Tested with:

#!/bin/bash
echo "empty BAM"
pipeval validate  test_files/empty_bam.bam

printf "\n"

echo "invalid BAM"
pipeval validate  test_files/invalid.bam

printf "\n"

echo "pass BAM"
pipeval validate  test_files/pass.bam

printf "\n"

echo "BAM with no index"
pipeval validate  test_files/noindex.bam

printf "\n"

echo "Just text file"
pipeval validate  test_files/hello.txt

printf "\n"

echo "Failing checksum MD5"
pipeval validate  test_files/hello_bad_md5.txt

echo "Failing checksum SHA512"
pipeval validate  test_files/hello_bad_sha512.txt

printf "\n"

echo "Generate md5 checksum"
pipeval generate-checksum -t md5 test_files/togen.txt

echo "Generate sha512 checksum"
pipeval generate-checksum -t sha512 test_files/togen.txt

echo "Validate generated checksums"
pipeval validate test_files/togen.txt

rm test_files/togen.txt.md5
rm test_files/togen.txt.sha512

printf "\n"

echo "Valid VCF"
pipeval validate  test_files/test_vcf.vcf.gz

printf "\n"

echo "Valid CRAM"
pipeval validate test_files/valid.cram -r /hot/ref/reference/GRCh38-BI-20160721/Homo_sapiens_assembly38.fasta

printf "\n"

echo "CRAM with no index"
pipeval validate test_files/noindex.cram -r /hot/ref/reference/GRCh38-BI-20160721/Homo_sapiens_assembly38.fasta

printf "\n"

echo "CRAM with default reference"
pipeval validate test_files/default_ref.cram

printf "\n"

echo "Invalid CRAM"
pipeval validate test_files/invalid.cram -r /hot/ref/reference/GRCh38-BI-20160721/Homo_sapiens_assembly38.fasta

printf "\n"

echo "Valid SAM"
pipeval validate test_files/valid.sam

printf "\n"

echo "Valid FASTQ"
pipeval validate test.fq.gz -t

printf "\n"

echo "Truncated FASTQ"
pipeval validate test_13.fastq

printf "\n"

echo "Bad record FASTQ"
pipeval validate test_invalid.fastq.bz2 -t

Checklist

File Commits

  • This PR does NOT contain Protected Health Information (PHI). A repo may need to be deleted if such data is uploaded.
    Disclosing PHI is a major problem1 - Even a small leak can be costly2.

  • This PR does NOT contain germline genetic data3, RNA-Seq, DNA methylation, microbiome or other molecular data4.

  • This PR does NOT contain other non-plain text files, such as: compressed files, images (e.g. .png, .jpeg), .pdf, .RData, .xlsx, .doc, .ppt, or other output files.

  To automatically exclude such files using a .gitignore file, see here for example.

Code Review Best Practices

  • I have read the code review guidelines and the code review best practice on GitHub check-list.

  • I have set up or verified the main branch protection rule following the github standards before opening this pull request.

  • The name of the branch is meaningful and well formatted following the standards, using [AD_username (or 5 letters of AD if AD is too long)]-[brief_description_of_branch].

  • I have added the major changes included in this pull request to the CHANGELOG.md under the next release version or unreleased, and updated the date.

Testing

  • I have added unit tests for the new feature(s).

  • I modified the integration test(s) to include the new feature.

  • All new and previously existing tests passed locally and/or on the cluster.

  • The docker image built successfully on the cluster.

Footnotes

  1. UCLA Health reaches $7.5m settlement over 2015 breach of 4.5m patient records

  2. The average healthcare data breach costs $2.2 million, despite the majority of breaches releasing fewer than 500 records.

  3. Genetic information is considered PHI.
    Forensic assays can identify patients with as few as 21 SNPs

  4. RNA-Seq, DNA methylation, microbiome, or other molecular data can be used to predict genotypes (PHI) and reveal a patient's identity.

@yashpatel6 yashpatel6 merged commit 290fa89 into main Dec 7, 2023
@yashpatel6 yashpatel6 deleted the yashpatel-compression-integrity-check branch December 7, 2023 22:07
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