Say I want to parse out all split reads. These are reads with a clipped position and a secondary mapping.
How do I identify these reads in samtools? Is there a flag to find reads that are split with a secondary mapping.
Once I have these reads, how I can I tell which chromosome they are mapped to?
The blue read should have a clipped signal, I want to parse out these reads that have the alternate mapping position on the same chromosome.