Question: What are the rules Pindel uses to determine which pairs of reads it will use?
1
gravatar for mohesan
3.0 years ago by
mohesan10
Canada
mohesan10 wrote:

I am trying to determine how Pindel decides on which (pair)reads it uses to analyse for structural variances. Pindel takes in a bam file that has been aligned using BWA. Pindel then applies the pattern-growth method,to find indels, on reads that could not be mapped but their pair was mapped successfully as well as on reads that BWA matched both pairs successfully but found an indel on them.  According to the BWA aligned bam files that I am working with, there are reads that support an indel but Pindel does not return these reads when I run it. My question is: What are the rules Pindel uses to determine which pairs of reads it will use? Any help at all will be appreciated! Thanks!

ADD COMMENTlink modified 3.0 years ago • written 3.0 years ago by mohesan10

For example should there be a certain MS score(i doubt this is the case)?

ADD REPLYlink written 3.0 years ago by mohesan10
0
gravatar for liangkaiye
3.0 years ago by
liangkaiye250
United States
liangkaiye250 wrote:

when you say Pindel does not use these reads, what kind of the reads you are referring to? Is it discordant read pair or the other type of difficult-to-map reads?

 

 

ADD COMMENTlink written 3.0 years ago by liangkaiye250

I am not sure. Is that a possible reason for them to remove reads?

ADD REPLYlink written 3.0 years ago by mohesan10

Majority of the reads that Pindel did not use were mapped in proper pair

ADD REPLYlink written 3.0 years ago by mohesan10
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: 867 users visited in the last hour