Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Does svtyper calculate read pairs which have long overlap twice? #100

Open
xyw1 opened this issue Feb 15, 2019 · 0 comments
Open

Does svtyper calculate read pairs which have long overlap twice? #100

xyw1 opened this issue Feb 15, 2019 · 0 comments

Comments

@xyw1
Copy link

xyw1 commented Feb 15, 2019

Hi Dave,
I read you script, and noticed that RO is calculated by ref_seq+ref_span. and ref_seq is calculated by the following code.

#  svtyper/svtyper/singlesample.py line 255 to 259 
is_ref_seq_A = sam_fragment.is_ref_seq(read, None, chromA, posA, ciA, min_aligned)
is_ref_seq_B = sam_fragment.is_ref_seq(read, None, chromB, posB, ciB, min_aligned)
if (is_ref_seq_A or is_ref_seq_B):
    p_reference = prob_mapq(read)
    ref_seq += p_reference

I'm wondering that if this read has a paired end read that overlaps with it a lot (shown as the illustration below), will svtyper take read1 and read2 as two independent observations?

image

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant