SvPileup should handle circular contigs better
Closed this issue · 1 comments
tfenne commented
When running SvPileup on a BAM file aligned carefully to a circular genome (plasmid, mito) it generates:
- A really well supported breakpoint using split-reads than span the origin
- A huge number of "abnormal pair" supported breakpoints
For intra-contig evidence, if the contig is labeled as circular (TP:circular
in the relevant SQ
header), sv-pileup should treat the transition around the origin as contiguous and re-calculate the "inner distance" between segments for split reads and between reads for pairs. Doing this should remove the vast majority of origin-related breakpoints being emitted.
jdidion commented
Do you have some example data I can use for test case(s)?