Large Deletions With Bwa And Single End Reads
0
2
Entering edit mode
12.8 years ago
Travis ★ 2.8k

Hi,

I am trying to detect deletions approximately 66bp in length by mapping 90 base reads to the human genome with BWA (and using GATK). By fine tuning the alignment parameters I have managed to increase my maximum detected indel length to around 30 bases but so far I have been unable to improve on that. I have increased the maximum number of gap extensions to 200 and reduced the gap extension penalty to zero so in my head at least, I thought the larger indels would be detectable. Can anyone offer any insight?

Thanks in advance.

indel next-gen sequencing alignment bwa • 5.9k views
ADD COMMENT
1
Entering edit mode

You might consider artificially splitting your read into two fragments in various lengths and mapping them separately? Tophat uses that strategy for transcript splicing.

ADD REPLY
0
Entering edit mode

Hi @Travis, are you trying to find deletions in the reference that are insertions in your sequencing sample, or deletions in your sample that are insertions in the reference, and split your reads into halves?

ADD REPLY
0
Entering edit mode

At the minute I am trying to find two deletions in the 60-66bp range. The deletion is in the reads, not the reference.

ADD REPLY
0
Entering edit mode

I have detected deletions in the 45-64b range using bwa/gatk (100bp PE Illumina), but haven't had them validated so.. I don't want to say much more than that. I was under the impression from recent conference talks that identification of deletions over 30bp (1/3rd read length) wasn't likely to happen with this particular toolchain.

ADD REPLY
0
Entering edit mode

Cheers Daniel. I am attempting to tweak another parameter or two so will see what happens on that front and update. I am determined to do it somehow!

ADD REPLY
0
Entering edit mode

Good advice, but I am specifically attempting to address the program without customization of anything other than existing program parameters.

ADD REPLY
0
Entering edit mode

You could also try other aligners such as blat, gsnap, and bfast.

ADD REPLY
0
Entering edit mode

Bfast will be my next port of call if I can't get the desired resuls from bwa. I spotted Nils commenting somewhere about the fact that his tests showed better indel detecting abilities than with bwa. The only drawback is speed.

ADD REPLY
0
Entering edit mode

Bfast will be my next port of call if I can't get the desired results from bwa. I spotted Nils commenting somewhere about the fact that his tests showed better indel detecting abilities than with bwa. The only drawback is speed.

ADD REPLY
0
Entering edit mode

Does anyone know what is limiting BWA's ability to find a 63bp deletion when it can find a 30bp deletion without a problem? Gap extension penalties are turned off.

ADD REPLY

Login before adding your answer.

Traffic: 762 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6