-
Notifications
You must be signed in to change notification settings - Fork 596
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 htsjdk to 2.21.0 #6250
update htsjdk to 2.21.0 #6250
Conversation
change in UpdateVCFSequenceDictionary to prevent 0 length sequences
dict.getSequences().stream() | ||
.filter( s -> s.getSequenceLength() == SAMSequenceRecord.UNKNOWN_SEQUENCE_LENGTH) | ||
.map(SAMSequenceRecord::getSequenceName) | ||
.collect(Collectors.joining(","))); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Do you think we should use .limit(20)
or something? There might be a lot...
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
that's a good point
@@ -209,11 +211,22 @@ else if (hasReference()) { | |||
if (resultDictionary == null || resultDictionary.getSequences().isEmpty()) { | |||
throw new CommandLineException.BadArgumentValue( | |||
String.format( | |||
"The specified dictionary source has an empty or invalid sequence dictionary", | |||
"The specified dictionary source has an empty or invalid sequence dictionary: %S", |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is the upper casing deliberate %S
?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
oh.. no.. accidental
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
LGTM once tests pass.
Upgrading htsjdk to 2.21.0.
This includes a change that relaxes restrictions on loading vcfs with sequence dictionaries that are missing lengths. This effected one of the tests which I changed a bit. We could also implement blanket ban on sequence dictionaries with missing lengths in our sequence dictionary validation code, but we currently allow them for sam/bam as far as I can tell and have tests that take them into account.