-
Notifications
You must be signed in to change notification settings - Fork 1
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
Conversation
|
||
override def execute(): Unit = { | ||
val source = SamSource(input) | ||
val writer = bam.map(path => SamWriter(path, header=source.header, sort=Some(SamOrder.Queryname))) |
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.
fixme: sort is not needed if input is in the correct order.
6b8b5d9
to
9c13392
Compare
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.
Have yet to look through tests of FilterAndMerge, but wanted to stop and talk about a few things.
|
||
// HACK | ||
SamSource.DefaultUseAsyncIo = true | ||
SamWriter.DefaultUseAsyncIo = true | ||
VcfWriter.DefaultUseAsyncIo = true |
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.
Huh, maybe we should implement a better way to override defaults in FgBioMain
?
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.
Agreed, leaving that to later
src/main/scala/com/fulcrumgenomics/fgsv/util/PutativeBreakpoint.scala
Outdated
Show resolved
Hide resolved
src/main/scala/com/fulcrumgenomics/fgsv/util/PutativeBreakpoint.scala
Outdated
Show resolved
Hide resolved
src/main/scala/com/fulcrumgenomics/fgsv/util/PutativeBreakpoint.scala
Outdated
Show resolved
Hide resolved
import htsjdk.samtools.util.CoordMath | ||
|
||
/** A genomic interval and useful methods. */ | ||
case class GenomicRange(refIndex: Int, start: Int, end: Int) extends Ordered[GenomicRange] { |
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.
Can we push this back down into fgbio please?
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.
todo
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.
On re-review, I think we need to think about how to handle refIndex vs. refName more generally before doing that.
src/main/scala/com/fulcrumgenomics/fgsv/util/GenomicRange.scala
Outdated
Show resolved
Hide resolved
Oh, also meant to say in a comment somewhere, that I think the |
|
||
// HACK | ||
SamSource.DefaultUseAsyncIo = true | ||
SamWriter.DefaultUseAsyncIo = true | ||
VcfWriter.DefaultUseAsyncIo = true |
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.
Agreed, leaving that to later
src/main/scala/com/fulcrumgenomics/fgsv/util/PutativeBreakpoint.scala
Outdated
Show resolved
Hide resolved
src/main/scala/com/fulcrumgenomics/fgsv/util/PutativeBreakpoint.scala
Outdated
Show resolved
Hide resolved
src/main/scala/com/fulcrumgenomics/fgsv/util/GenomicRange.scala
Outdated
Show resolved
Hide resolved
import htsjdk.samtools.util.CoordMath | ||
|
||
/** A genomic interval and useful methods. */ | ||
case class GenomicRange(refIndex: Int, start: Int, end: Int) extends Ordered[GenomicRange] { |
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.
todo
import scala.collection.immutable | ||
|
||
/** An enumeration over types of SV evidence. */ | ||
object EvidenceType extends FgBioEnum[EvidenceType] { |
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.
todo
* 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.
3b720d6
to
2626d06
Compare
No description provided.