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

SvPileup should allow some minimal filtering #26

Open
tfenne opened this issue Sep 28, 2023 · 1 comment
Open

SvPileup should allow some minimal filtering #26

tfenne opened this issue Sep 28, 2023 · 1 comment

Comments

@tfenne
Copy link
Member

tfenne commented Sep 28, 2023

I'm envisioning we could add parameters like:

--min-split-reads (default: 0)
--min-read-pairs (default: 0)
--min-total-support (default: 1)

The only problem with doing this is that we would no longer be able to write the evidence BAM on the fly - we'd have to write a temporary BAM with all reads supporting all break-points, and then re-process that to remove reads/templates for breakpoints with insufficient support. That's not a big overhead - these BAMs are usually pretty small. And it would also give us a logical place to add re-sorting the BAM file into coordinate order for use in IGV.

@tfenne
Copy link
Member Author

tfenne commented Sep 28, 2023

Alternatively this could be a separate tool FilterPileup that takes in the raw pileup, filters the text file and then filters the BAM file emitted by SvPileup.

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

No branches or pull requests

1 participant