I was trying to use the hla-ligand-atlas branch for HLA type II calling however I was running into errors when running, after step 4 the process exits. Here is the log file I got for my samples, but I got the same issue with the example files.
[E::idx_find_and_load] Could not retrieve index file for 'mapped/test/840113-107/r1.bam'
0:00:01.30 imports ready
0:00:01.30 start bam loading
aln read 0/0 ...Traceback (most recent call last):
File "py/enumtyper.py", line 46, in
b1 = bam_to_hdf(args.bam1, min_readlength=30)
File "/home/szabolcs_tarapcsak_servier_com/LocalSoftware/ot2-snake/py/otutils.py", line 281, in bam_to_hdf
read_hits = aln_to_hits(aln, i_seq, n_deduction=aln.query_sequence.count('N'))
File "/home/szabolcs_tarapcsak_servier_com/LocalSoftware/ot2-snake/py/otutils.py", line 225, in aln_to_hits
refname, start, cigar, nm = subtag.split(',')
ValueError: too many values to unpack (expected 4)
Szia Andras,
I was trying to use the hla-ligand-atlas branch for HLA type II calling however I was running into errors when running, after step 4 the process exits. Here is the log file I got for my samples, but I got the same issue with the example files.
[E::idx_find_and_load] Could not retrieve index file for 'mapped/test/840113-107/r1.bam'
Thank you for your help! Best, Szabi