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

possorted_bam.bam versus phased_possorted_bam.bam #28

Open
burrir opened this issue Feb 10, 2022 · 2 comments
Open

possorted_bam.bam versus phased_possorted_bam.bam #28

burrir opened this issue Feb 10, 2022 · 2 comments

Comments

@burrir
Copy link

burrir commented Feb 10, 2022

Hi,
I would be curious to hear whether you expect the possorted_bam.bam output by "LongRanger align" (rather than the phased_possorted_bam.bam output by "LongRanter wgs") should work adequately in linkedSV?

We are looking to replace LongRanger wgs, as it is tremendously time consuming, and to our understanding, LongRanger align provides a phased output as well. We yet fail to understand the difference between the bam files output by the two different pipelines.

Any insights would be very appreciated!

Best wishes,
Reto

@fangli80
Copy link
Collaborator

Hello Reto,
Yes. I understand that the Longranger pipeline is pretty slow.
I have not tested on the "possorted_bam.bam". Does the bam file has a "BX" tag for barcode and an "HP" tag for haplotype?

Thanks,
Li

@burrir
Copy link
Author

burrir commented Feb 14, 2022

Hi Li,

thanks a lot for getting back on this question! The merely possorted bam files contain the BX field. This is what e.g. hapCUT2 uses to infer phase information (which is why I figured that there is phase info, but it's the info to infer phase and not phase info itself). However, as I now saw when checking, these files do not feature the HP tag. Thus the difference in naming and why using the possorted_bam will not work for linkedSV.

Thanks again & best wishes!
Reto

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

2 participants