Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Bwa Soft Clipped Indels Realignment
0
gravatar for stolarek.ir
6.2 years ago by
stolarek.ir • 640
Poland
stolarek.ir • 640 wrote:

Hi guys.

I generated a sam file by aligning single end reads with bwa-mem.

i have identified softclipped fragments which suggest indels at their breakpoints. Now I wanted to select those fragments and realign around those regions. And here is my question:

Should I take left and right breakpoint softclipped fragments putting them together in one file as one would normally do for pair-end data? And them using those 2 files to realign against reference file composed of the regions around the breakpoints?

I never tried indel realignment. And my idea is utilizing soft clipped fragments, which to complicate matters more come from single end reads.

Can any one suggest anything?

• 2.2k views
ADD COMMENT • link •
modified 6.2 years ago • written 6.2 years ago by stolarek.ir • 640
0
gravatar for Ashutosh Pandey
6.2 years ago by
Ashutosh Pandey ♦ 11k
Philadelphia
Ashutosh Pandey ♦ 11k wrote:

Have you ever tried using SVMerge (http://svmerge.sourceforge.net/). It exactly does what you are looking for. Little complicated to use but it does the local realignment around the softclipped reads to get you the exact breakpoints.

ADD COMMENT • link written 6.2 years ago by Ashutosh Pandey ♦ 11k

Hmm actually it depends on the calls made by other tools in the first step. Only Pindel is using proper split read approach in which I'm interested (but don't know if for single end ones). Other tools are no good for me, as my pipeline has much greater sensitivity and base pair resolution in first step of analyzing SAM file.

I now it sounds very geeky, that I'm writing almost everything from scratch, but well, that's the best way to learn something.

I think that it shouldn't matter now if I will use soft clipped parts as pair end or just dump into one file. Proper naming will resolve this problem, and adding +- 1k bp either side of the breakpoint in reference should give nice specificity

Anyways, thanks for the SVMerge

ADD REPLY • link written 6.2 years ago by stolarek.ir • 640
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2401 users visited in the last hour