-
Notifications
You must be signed in to change notification settings - Fork 18
Description
Hi,
I want to use polysolver for HLA-typing, thank you very much for developing the conda version of polysolver! (I am not root user so cannot install it using docker as suggested by Broad Institute). I install the polysolver following your guide, all steps seems working well. But when I test on test.bam, there is difference between my output and the standard output.
(polysolver) [shiyang@statcomp hla-polysolver-master]$ diff output/winners.hla.txt test/orig.winners.hla.txt
2,3c2,3
< HLA-B hla_b_39_01_01_02l hla_b_39_01_01_03
< HLA-C hla_c_07_01_05 hla_c_07_01_05
---
> HLA-B hla_b_39_01_01_02l hla_b_39_01_01_02l
> HLA-C hla_c_07_01_05 hla_c_06_02_01_01
And below is a section of my standard output that seems wrong when I run test script:
[samopen] SAM header is present: 6597 sequences.
Tue Jul 10 22:05:17 CST 2018
get first winners
Tue Jul 10 22:05:17 CST 2018
rm: cannot remove 'output/counts1.R0k6': No such file or directory
Tue Jul 10 22:05:26 CST 2018
calculating lik2
winners1 hla_a_24_02_01_01 hla_b_39_01_01_02l hla_c_07_01_05
Tue Jul 10 22:05:28 CST 2018
get second winners
rm: cannot remove 'output/counts2.R0k6': No such file or directory
winners1 hla_a_24_02_01_01 hla_b_39_01_01_02l hla_c_07_01_05
winners2 hla_a_24_02_01_01 hla_b_39_01_01_03 hla_c_07_01_05
cleanup
Do you know which cause this discrepancy? Could you be kind to give me some advice?
Thanks!
Yang