missed supporting reads

[helloeoh001]

Hello,

I'm using pindel to detect tandem duplication frequently occurring in FLT3 gene in leukemia.

The size of duplicated sequence in my data is about 60bp, and the number of supporting read was about 140 when checking the bam file with IGH viewer.

However, only 90 of the them were detected by pindel, and most of the supporting reads on the left side were missed (refer to the figure below).

I need correct number of supporting reads to calculate VAF.

Please help me~

Thanks in advance

[liangkaiye]

Is it possible to share the local bam file to look into this? Kaiye@xjtu.edu.cn<mailto:Kaiye@xjtu.edu.cn> Hello, I'm using pindel to detect tandem duplication frequently occurring in FLT3 gene in leukemia. The size of duplicated sequence in my data is about 60bp, and the number of supporting read was about 140 when checking the bam file with IGH viewer. However, only 90 of the them were detected by pindel, and most of the supporting reads on the left side were missed (refer to the figure below). I need correct number of supporting reads to calculate VAF. Please help me~ Thanks in advance [image]<https://user-images.githubusercontent.com/16658129/38807046-36749a74-41b6-11e8-86af-49329d5c82ac.png> — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub<#88>, or mute the thread<https://github.com/notifications/unsubscribe-auth/AB9s-4LsxlT5FXInKMUP_inuoWaU_kWyks5tpIN6gaJpZM4TWXJ0>.

[philip-holmgren]

Hi,

I think this is related to the issue I raised (#86).

In that case the reason why the left reads were 'missed' is because they are used as the 'mapped' read to find the 'unmapped' read which would be the right read of the pair (using 'mapped' and 'unmapped' here in reference to the description of the Pindel algorithm).

As far as I know the 'mapped' read is never taken into account for Pindel's VAF calculation although it's quite possible that read also contains the ITD and is softclipped.

[helloeoh001]

Dr. Ye,

I e-mailed you with the link that the bam file can be downloaded.

If there are any troubles in downloading, let me know.

[liangkaiye]

Got the files. I am currently working on a proposal and will start to work on the code after that. Closed #88<#88>. — You are receiving this because you commented. Reply to this email directly, view it on GitHub<#88 (comment)>, or mute the thread<https://github.com/notifications/unsubscribe-auth/AB9s-whA4d9g-Fy71Eyw-z3IrXvUoCJvks5tpaF0gaJpZM4TWXJ0>.