8000 Haplotag command only tagging a small percentage of reads! · Issue #97 · twolinin/longphase · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

Haplotag command only tagging a small percentage of reads! #97

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.

Sign up for GitHub

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

Jump to bottom
Open
DHmeduni opened this issue Jan 14, 2025 · 2 comments
Open

Haplotag command only tagging a small percentage of reads! #97

DHmeduni opened this issue Jan 14, 2025 · 2 comments

Comments

@DHmeduni
Copy link

Hi,
The phasing vcf looks good, but only a small percentage of reads are being haplotagged...I do have a deep coverage (>100x)
The same data with whatshap is able to phase and haplotag alot more.
the haplotag command for reference.

haplotag -b /data/test2/XXXXXXXX_bed.bam -s /data/test2/phased_prefix_SV.vcf -r /data/#references/chm13v2.0/chm13v2.0.fa --sv-file /data/test2/phased_prefix_SV.vcf -o /data/test2/haplotag --log -t 32

Any help would be must appreciated.
Thanks,
DH
@ythuang0522
Copy link
Collaborator

Hi @DHmeduni

From the command it looks like you provided the wrong SNP vcf file (phased_prefix_SV.vcf), assuming phased_prefix_SV.vcf only stores the phased SVs. We need both SNP and SV vcfs provided for haplotagging, because the number of SVs is way much less than SNPs. Can you check this first?

Yao-Ting

@DHmeduni
Copy link
Author

I totally oversaw that!...WOW!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet < 3B5E /div>
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
2 participants
@ythuang0522 @DHmeduni
0