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

Merge overlapping read pairs into single LocusRead #59

Open
iskandr opened this issue Nov 4, 2016 · 1 comment
Open

Merge overlapping read pairs into single LocusRead #59

iskandr opened this issue Nov 4, 2016 · 1 comment

Comments

@iskandr
Copy link
Contributor

iskandr commented Nov 4, 2016

If the insert size is small enough then many reads actually overlap the sequence of their mate pairs, which provides an opportunity for merging them into a longer read length. Doing this would also prevent over count of variant allele support.

@iskandr
Copy link
Contributor Author

iskandr commented Nov 18, 2016

May be simplified when we switch to using pysam's fetch instead of pileups (#64) since we can first group all reads at a locus by their template name and then when processing each group merge multiple reads which come from the same fragment.

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