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
hapdup fails with multiple primary alignments #25
Comments
Never seen this before, but looks like the error might be inside PEPPER/Margin packages. Do you mind sending |
I don't believe there to be anything particularly unusual about the sample, but it does have some ultra long read data in it. The log file is below. I don't have the full terminal output but here is where it was at before the error:
|
Thanks. According to this, it might have something to do with the duplicated read ids: kishwarshafin/pepper#67 Could you try command line: |
Apologies - I've realised this was an error on our side where duplicated reads had got into the assembly pipeline. Thanks for your help. |
I've got one sample in a series which is failing in hapdup with the following error.
Any suggestions?
Thanks
`
Expected three tokens in header line, got 2
This usually means you have multiple primary alignments with the same read ID.
You can identify whether this is the case with this command:
[2022-12-13 17:15:57] ERROR: Missing output: hapdup/margin/MARGIN_PHASED.haplotagged.bam
Traceback (most recent call last):
File "/data/test_data/GIT/hapdup/hapdup.py", line 24, in
sys.exit(main())
File "/data/test_data/GIT/hapdup/hapdup/main.py", line 206, in main
file_check(haplotagged_bam)
File "/data/test_data/GIT/hapdup/hapdup/main.py", line 114, in file_check
raise Exception("Missing output")
Exception: Missing output
`
The text was updated successfully, but these errors were encountered: