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

Add a tool to make IGV-compat. BEDPE files from AggregateSvPileup #37

Merged
merged 8 commits into from
Apr 10, 2024

Conversation

clintval
Copy link
Member

Hopefully this addition is welcome! Converting the output of AggregateSvPileup to BEDPE facilitates genome-browsing in IGV and helps with understanding how breakpoint events relate to the alignments that cause them. There exist more advanced features of BEDPE in IGV (like record-specific coloring) but I chose not to implement that at this time to keep the tool simple and forwards compatible. Here is what test data looks like in IGV:

Screenshot 2024-04-10 at 12 38 26 PM

@clintval clintval requested a review from nh13 April 10, 2024 21:12
Comment on lines -8 to +19
class _All extends ClpGroup {
override val name: String = "All tools"
override val description: String = "All tools."
class _BreakpointAndSv extends ClpGroup {
override val name: String = "Breakpoint and SV Tools"
override val description: String = "Primary tools for calling and transforming breakpoints and SVs."
}

final val All = classOf[_All]
class _Utilities extends ClpGroup {
override val name: String = "Utility Tools"
override val description: String = "Helper tools for working with breakpoint or SV data."
}

final val BreakpointAndSv = classOf[_BreakpointAndSv]
final val Utilities = classOf[_Utilities]
Copy link
Member Author

Choose a reason for hiding this comment

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

This is a backwards compatible change and make the CLI look like:

Screenshot 2024-04-10 at 12 38 05 PM

Comment on lines +72 to +75
end1 = pileup.left_max_pos + 1,
chrom2 = pileup.right_contig,
start2 = pileup.right_min_pos,
end2 = pileup.right_max_pos + 1,
Copy link
Member Author

@clintval clintval Apr 10, 2024

Choose a reason for hiding this comment

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

Max positions are 0-based so to make them half-open we need to add 1.

EDIT: This was wrong! See #40

Copy link
Member

@nh13 nh13 left a comment

Choose a reason for hiding this comment

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

nitpicks, but lgtm

@clintval clintval assigned clintval and unassigned nh13 Apr 10, 2024
@clintval clintval merged commit 30baecb into main Apr 10, 2024
@clintval clintval deleted the cv_bedpe branch April 10, 2024 21:35
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