fgsv
fgsv copied to clipboard
SvPileup should allow some minimal filtering
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.