-
Notifications
You must be signed in to change notification settings - Fork 3
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
NimaGen 29-plex #134
NimaGen 29-plex #134
Conversation
markers = sorted(interval.data, key=lambda m: m.sortkey) | ||
markers = sorted(interval.data, key=lambda m: m.sources[0].sortkey) |
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.
Fixing a bug that was introduced in #133.
def rsid_union(self, *others): | ||
rsids = set(self.rsids) | ||
for other in others: | ||
rsids |= set(other.rsids) | ||
for rsid in rsids: | ||
for marker in chain([self], others): | ||
if rsid not in marker.rsids: | ||
marker.rsids.append(rsid) |
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.
Some markers are defined only using positions, no RSIDs. Some of these are exact matches for markers with some or all RSIDs defined. This function merges RSID lists when markers are merged during the build.
This PR adds definitions for 29 markers from NimaGen's (currently beta) microhap kit.