Question: How to filter paired-end reads to detect inversions according to insert size
0
gravatar for wang.yiguan
13 months ago by
wang.yiguan0 wrote:

Hi,

I have paired-end NGS data of a fruitfly population, and I am trying to detect inversions according to paired-end insert size under the premise that a much larger insert size will be observed in presence of inversions. ( A break point between pair-end reads will increase the insert size when aligned to reference genome)

But I find the situation is more much complicated that I expected. The reads can be mapped in different ways, eg. supplementary alignment or chimeric reads... I also noticed the sam flag(second column in sam files) provides such information, but I am not clear how to filter reads according to these flags.

My question is: in order to infer inversions based on insert size, how should I filter reads?

Thanks in advance!

alignment next-gen genome • 322 views
ADD COMMENTlink written 13 months ago by wang.yiguan0
2

But I find the situation is more much complicated that I expected

Exactly, this is not trivial at all and I strongly encourage to use dedicated software for it, such as lumpy (or any other structural variant caller). Naive approaches might work but will take time to develop, must be tested etc... Use standard software!

ADD REPLYlink written 13 months ago by ATpoint40k
Please log in to add an answer.

Help
Access

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