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

Initial bolus of code #1

Merged
merged 3 commits into from
Feb 14, 2022
Merged

Initial bolus of code #1

merged 3 commits into from
Feb 14, 2022

Conversation

nh13
Copy link
Member

@nh13 nh13 commented Jan 31, 2022

No description provided.

@nh13 nh13 marked this pull request as draft January 31, 2022 16:29

override def execute(): Unit = {
val source = SamSource(input)
val writer = bam.map(path => SamWriter(path, header=source.header, sort=Some(SamOrder.Queryname)))
Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

fixme: sort is not needed if input is in the correct order.

@nh13 nh13 requested a review from tfenne February 4, 2022 18:09
@nh13 nh13 force-pushed the feature/initial-tools branch from 6b8b5d9 to 9c13392 Compare February 4, 2022 18:15
Copy link
Member

@tfenne tfenne left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Have yet to look through tests of FilterAndMerge, but wanted to stop and talk about a few things.

Comment on lines 16 to 20

// HACK
SamSource.DefaultUseAsyncIo = true
SamWriter.DefaultUseAsyncIo = true
VcfWriter.DefaultUseAsyncIo = true
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Huh, maybe we should implement a better way to override defaults in FgBioMain?

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Agreed, leaving that to later

import htsjdk.samtools.util.CoordMath

/** A genomic interval and useful methods. */
case class GenomicRange(refIndex: Int, start: Int, end: Int) extends Ordered[GenomicRange] {
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Can we push this back down into fgbio please?

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

todo

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

On re-review, I think we need to think about how to handle refIndex vs. refName more generally before doing that.

@tfenne
Copy link
Member

tfenne commented Feb 8, 2022

Oh, also meant to say in a comment somewhere, that I think the util package is a bit of an anti-pattern. Would suggest either bumping those classes up a level or renaming that package core.

Comment on lines 16 to 20

// HACK
SamSource.DefaultUseAsyncIo = true
SamWriter.DefaultUseAsyncIo = true
VcfWriter.DefaultUseAsyncIo = true
Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Agreed, leaving that to later

import htsjdk.samtools.util.CoordMath

/** A genomic interval and useful methods. */
case class GenomicRange(refIndex: Int, start: Int, end: Int) extends Ordered[GenomicRange] {
Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

todo

import scala.collection.immutable

/** An enumeration over types of SV evidence. */
object EvidenceType extends FgBioEnum[EvidenceType] {
Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

todo

nh13 and others added 3 commits February 14, 2022 16:17
* Attempting to fix breakpoint construction from alignment blocks.
* Clean up after a messy rebase.
A bunch of refactoring, most notably of how we accrue counts and assign IDs to breakpoints, so that we are guaranteed to count all evidence for breakpoint the same way and only assign it a single ID.
@tfenne tfenne force-pushed the feature/initial-tools branch from 3b720d6 to 2626d06 Compare February 14, 2022 23:17
@tfenne tfenne marked this pull request as ready for review February 14, 2022 23:17
@tfenne tfenne merged commit 974b92e into main Feb 14, 2022
@tfenne tfenne deleted the feature/initial-tools branch February 14, 2022 23:18
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