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

update picard 2.18.11 -> 2.20.3 #977

Merged
merged 34 commits into from
Sep 11, 2019
Merged

update picard 2.18.11 -> 2.20.3 #977

merged 34 commits into from
Sep 11, 2019

Conversation

tomkinsc
Copy link
Member

@tomkinsc tomkinsc commented Jul 11, 2019

This PR does the following:

  • updates picard 2.18.11 -> 2.20.3

compare only reads in sam post-depletion, allowing us to ignore header differences due to varying versions of picard/samtools/htslib; remove related permutations of expected output from tests
@dpark01
Copy link
Member

dpark01 commented Jul 11, 2019

Wait, regarding this commit, maybe instead of "fixing" the test input file, should we be adding code that sanitizes input all the time now? Because apparently this wasn't an issue for old versions of Picard, but now it is, so these changes will break things that used to work?

@tomkinsc
Copy link
Member Author

tomkinsc commented Jul 11, 2019

Here's the (human-data-centric) discussion upstream about the breaking changes to htsjdk:
samtools/hts-specs#333
Allowing certain characters in the reference name seems put resulting files out of spec with the SAM standard. We can sanitize if you think that's best. Or maybe if the reference sequences we include are from GenBank including the accession in addition to the sequence name would resolve the issue and be more 'correct.'

add new test assertion, "assertEqualSamHeaders" to address this inline ToDo for TestFastqBam(): "Note for developers: if you're fixing the tests to handle non-bugs (ie our testing here is too brittle), let's just replace a lot of this in the future with code that just reads the header, sorts it, and tests for equality of sorted values in the RG line (and stricter equality in the non-header lines). This is kind of hacky."
add debug statements to install-viral-ngs.sh Docker script
required by GATK 3.8 conda package
….file to sanitize fasta file IDs

fasta ID is sanitized for picard CreateSequenceDictionary to adhere to character set restrictions in SAM/BAM RNAME spec. see: samtools/hts-specs#333

Roll back test input now that we are doing the sanitization. As of this commit, sanitization is only performed on fasta files destined for CreateSequenceDictionary
change lz4 dependency from lz4-bin from bioconda to lz4-c from conda-forge, v131->v1.9.1
@tomkinsc tomkinsc merged commit 1dd180d into master Sep 11, 2019
@tomkinsc tomkinsc deleted the ct-update-picard-2.20.3 branch September 11, 2019 18:54
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