Question: Structural Variation; Discordant Pairs?
gravatar for Rm
8.3 years ago by
Danville, PA
Rm7.9k wrote:

For structural variation (SV) study: Following hydra or SVdetect or other pipelines;

I Need some clarification on type(s) of discordant reads from a mate pair or paired end analysis to be retained as inputs for the above programs.

We remove concordant reads (both forward and reverse reads aligned within the permisible range (depends on insert size) to a chromosome. (read mapped in proper pair)

For filtering a bam file what "samtools view" flags to be used? Some I come accross are

  1. -F 2
  2. -F 1806

aditionally we can also apply quality threshold too.

My question is whether to retain mate unmapped or singletons? as discordant reads.

All inputs on this topic is highly appreciated

ADD COMMENTlink modified 8.3 years ago by Chris Whelan550 • written 8.3 years ago by Rm7.9k
gravatar for Chris Whelan
8.2 years ago by
Chris Whelan550
Portland, OR
Chris Whelan550 wrote:

For "whether to retain mate unmapped or singletons as discordant reads", I think it is going to depend on what the SV caller you are using expects. I don't believe that either Hydra or SVDetect use pairs in which only one end is mapped as part of their algorithms (someone please correct me if I'm wrong on that!). Other tools, like Pindel, I believe, use those types of pairs. So I think you need to customize your workflow to the SV caller(s) you are planning to use.

ADD COMMENTlink written 8.2 years ago by Chris Whelan550
gravatar for Junfeng
8.2 years ago by
Junfeng310 wrote:

I would like to know whether you have solved the problems. The discordant pairs are the pairs with incorrect orientations or the pairs with mapped distance beyond a fixed range (Ruibin Xi, Tae-Min Kim and Peter J. Park. Detecting structural variations in the human genome using next generation sequencing. Briefings in Functional Genomics; 2010; 9(5-6): 405-415).

ADD COMMENTlink written 8.2 years ago by Junfeng310

thanks for the reference.

ADD REPLYlink written 8.2 years ago by Rm7.9k
Please log in to add an answer.


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