Entering edit mode
11.6 years ago
Tomáš Beluský
▴
90
Hi, I am implementing tool for detecting of structural variants and I am combining read pair and split read methods. The input is BAM with already mapped paired ends. When I find pair where one end is properly aligned and second is aligned with clipped sequence, I want to remap this sequence in surroundings of properly aligned end. So I need tool which can take these paired ends as input and do that for me. Or input can be clipped sequence and region where it should be mapped. So my question is if there is any tools, which can handle this. Also all other ideas are welcome. Thanks!
I am not sure what you want but Pindel (https://trac.nbic.nl/pindel/) does exactly what you are trying to do.
Hydra and Lumpy-SV are also good @ detecting structural variation. Writing your own code is always fun, but validation and benchmarking is not.
GASV (https://code.google.com/p/gasv/) also works great. Though it doesnt use the split read method but it uses combination of mate pair discordance and read depth analysis to find structural variants but unfortunately no exact breakpoints.
You're right, but this is my master thesis ;)
As I mentioned before that Pindel is one of the tool. Have you read about it? Lumpy-sv is another tool. i don't think there is some tool/module which will do something specific you want.
Yes, I heard about it, but I want to detect variations by myself. All I need is remap clipped sequences. But finally I have solution. Inspired by ClipCrop, I am cutting clipped sequence with a minimal length and try to remap it with bwa on same chromosome, where read was found. Thanks everybody for your help!
Hi Tomas, I knew the post is very old. But can I know your solution for Remap Clipped Sequence? I am trying to remap the Clipped sequence for detecting some eccDNA from our NGS data. My strategy is exactly like yours. If I can follow some of your code to make it work on my data. That will be great help for my project! Thank you!