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

Can't find first and/or second of pair in sam block of length 2 #389

Open
jonahcullen opened this issue Jan 18, 2023 · 0 comments
Open

Can't find first and/or second of pair in sam block of length 2 #389

jonahcullen opened this issue Jan 18, 2023 · 0 comments

Comments

@jonahcullen
Copy link

Hello, thank you for producing and supporting this tool, it is most excellent!

I have been struggling with a problem that seems to crop up for a small proportion of all the samples I attempting to process, namely the samblaster error below

samblaster: Can't find first and/or second of pair in sam block of length 2 for id: A00351:169:HLCY7DSXX:3:1220:15591:9862

I am running lumpyexpress as

lumpyexpress -B sample.bam -o sample.lumpy.vcf

The analysis-ready, coordinate-sorted sample.bam was generated following GATK's best practices from 16 paired FASTQs. When I look at the region where the error occurred I think I understand what it is complaining about but am not entirely sure

A00351:169:HLCY7DSXX:3:1220:15591:9862  163     chr1    1452 ... RG:Z:D06257_J
A00351:169:HLCY7DSXX:3:1220:15591:9862  163     chr1    1452 ... RG:Z:D06257_B
A00351:169:HLCY7DSXX:3:2220:16215:7592  163     chr1    1452 ... RG:Z:D06257_J
A00351:169:HLCY7DSXX:3:2220:16215:7592  163     chr1    1452 ... RG:Z:D06257_B

I was hoping to avoid pulling the split and discordant reads prior to input lumpy as I've already processed a couple hundred samples without issue. Is there some other filtering I should apply to sample.bam or should I rethink/reprocess everything to be consistent?

Thanks again,
Jonah.

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