Remap Clipped Sequence
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!
• 3,548 views
•
link
written
by
Tomáš Beluský •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Which means the default library orientation?
written by lee.jingu0626 •
-
Command Line Tools To Do Long Sequence Alignment
written by Xingyu Yang
-
In Paired-End Data, If One Read Of The Pair Is Unmapped, Is That Pair Autmoatically Improper?
written by mrals89
-
Hg19 Strand Information Of Sam Output.
written by Coryza •
-
How To Rearrange Paired End Bam File?
written by deepthithomaskannan
-
Find Right Variations In Sam/Bam File
written by Tomáš Beluský •
-
Detecting Rearranged Reads From Sam/Bam File
written by Tomáš Beluský •
-
Represent Precise Deletion In Vcf
written by Tomáš Beluský •
-
What Type Of Filtering Is Typically Applied To A Bam File?
written by Ashutosh Pandey
-
Counting The Whole Insert Size From Paired-End Reads As Coverage
written by Alastair Kerr
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!